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CHAPTER  1 


Introduction  And  Executive  Summary 

This  document  is  the  final  report  on  novel  techniques  for 
reducing  the  ECM  vulnerability  of  tactical  communications. 
Tactical  communications  have  been  plagued  by  having  had  so 
many  constraints  placed  on  them  that  there  has,. in  the  past, 
been  little  attention  paid  to  the  vulnerability  of  field  sys¬ 
tems  to  electronic  countermeasures.  Army  tactical  communica¬ 
tions  are  applied  in  a  wide  variety  of  circumstances  and  hostile 
environments.  The  systems  employed  include  field  radio  HF, 

VH?  and  UHF  single  channel  and  multichannel  and  SHF  satellite 
communications  -both  fixed  and  mobile.  In  addition,  there  are 
line  of  sight  microwave  radio  systems  and  UHF  tropospheric 
scatter  over  the  horizon.  The  deployment  of  various  tactical 
field  radios  is  prodigious  and  each  system  must  frequently 
operate  in  a  grid  network  tied  to  a  backbone  transmission  sys¬ 
tem  and  switching  centers.  The  application  on  the  tactical 
communication  is  widely  varied  and  includes  voice,  data,  mess¬ 
age  store  and  forward.  Any  or  all  of  the  information  may  be 
analog  or  digital  and  may  or  may  not  be  encrypted. 

Vulnerability  to  jamming  of  such  tactical  communication 
systems  has  been  recognized  for  some  time,  but  because  of  con¬ 
flicting  restraints  cost,  and  the  apparent  unavailability  of 
a  sufficiently  small,  inexpensive,  and  reliable  technology  to 
implement  ECCM  techniques,  these  techniques  have  lagged  behind. 

ECCM  techniques  come  in  many  different  forms  and  they  in- 
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dude  various  modes  of  spread  spectrum  (?N  direct  sequence, 


frequency  and  time  hopping,  chirp,  etc.),  signal  cancellers, 
diversity,  coding,  hall  steering  antennas,  power  control,  and 
transmission  coordination  among  the  favorites.  In  this  study, 
where  we  concentrate  on  multichannel  digital  communications  we 
address  some  ideas  which  would  appear  to  have  the  greatest 
possible  impact  on  ECCM  effectiveness  given,  the  constraints 
imposed  by  the  tactical  scenario  and  even  more  specifically  by 
the  nature  of  tactical  multichannel  communications.  These 
constraints  form  a  checklist  of  issues  in  tactical  multichannel 
communications  which  guided  our  thinking  and  they  are: 

•  Small,  stationary,  isolated  transmitters 

•  Relatively  low  cost 

•  Limited  bandwidth  availability  and  preassigned  frequency 
allocation 

•  The  quality  of  transmission  may  be  allowed  to  depend  on 
jamming  conditions  and  on  the  priority  class  of  the  user 

•  The  availability  of  channels  and  the  concomitant  grade  of 

service  (G.Q.S.)  may  be  allowed  to  depend  on  jamming  conditions 

and  on  the  priority  class  of  the  user 

/ 

•  Highly  sophisticated  jammers  are  unlikely  Cit  could  be 
easier  to  destroy  the  antenna  tower) 

•  There  may  exist  a  geometric  distribution  of  both 
friendly  and  unfriendly  interference. 

This  checklist,  while  formidable  in  the  limitations  that 
it  imposes  on  what  is  achievable,  does  nevertheless  narrow 


the  focus  on  those  ECCM  techniques  which  stand  a  chance  for 


success.  The  bright  side  of  the  picture  is  the  anticipation 
that  VHSI  and  VLSI  technique  which  will  be  available  by  1985 
will  allow  the  practical  implementation  of  many  of  the  novel 
ideas  which  today  would  appear  too  sophisticated.  With  this 
in  mind  we  have  chosen  for  investigation  six  ideas  which  appear 
to  be  promising.  They  are: 

•  Channel  monitoring  and  adaptive  channel  capacity 
allocation  (Chapters  2  and  3) 

•  Means  for  protecting  the  synchronization  of  multi¬ 
channel  digital  communications  (Chapter  4) 

•  The  use  of  low  rate  block  and  convolutional  codes  in 
lieu  of  traditional  spread  spectrum  (Chapter  5) 

•  Adaptive  null-steering  antenna  arrays  (Chapter  6) 

•  Bandwidth  reduction  with  multiple  alphabet  signalling 
(Chapter  7 ) 

•  Frequency  Hopping  systems  (Appendix  A) . 

The  following  is  a  brief  synopsis  of  each  Chapter  ox  this 
report. 


Chapter  2  -  Bit  Error  Rate  Monitoring  Using 
Extreme  Value  Theory 

In  order  to  be  able  to  react  to  a  jamming  condition,  it 
is  necessary  to  monitor  the  channel.  In  a  digital  system  the 
bit  error  rate  (3ER)  is  the  most  logical  and  directly  signif¬ 
icant  variable  to  observe.  Unfortunately,  a  working  system 
requires  a  very  low  3ER  (<  10* so  that  long  measurement 


times  are  often  necessary  to  obtain  an  estimate  of  3ER  with 
any  reasonable  confidence.  Mow  since  the  signal  and  noise 
structure  is  known,  it  is  possible  to  obtain  estimates  using 
extreme  value  theory  (EVT)  more  directly  and  more  quickly  (with 
fewer  measurements  than  would  appear  to  be  necessary) .  The 
monitoring  technique  described  here  is  relatively  simple  to 
implement  and  is  directly  amenable  for  use  in  conjuction  with 
the  dynamic  channel  allocation  scheme  described  in  Chapter  3. 

Chanter  3  -  Dynamic  Allocation  Of  Channel  Capacity 

Given  that  we  can  monitor  the  condition  of  the  channel  and 
determine  whether  or  not  a  jammer  is  present,  we  could,  ordi¬ 
narily  go  into  a  spread  spectrum  mode  to  provide  processing 
gain  (PG)  and  corresponding  jam  protection  for  communications. 
However,  since  mast  tactical  radios  are  bandwidth  limited  and 
tend  to  use  the  available  bandwidth  to  provide  6,  12,  or  IX 
or  more  voice  and  data  channels,  it  is  not  possible  to  merely 
spread  the  emissions  without  some  sacrifice.  One  way  to  ensure 
performance  and  protection  for  an  important  class  of  users  is 
to  use  the  bandwidth  of  possibly  less  important  users  during 
an  attack.  This  need  not  mean  that  the  less  important  users 
are  denied  access  to  communications  but  rather  that  their  grade 
of  service  (G.O.S.)  may  be  lowered  by  increasing  the  probability 
that  they  will  be  denied  access  to  the  service  (depending  on  the 
total  demand).  For  example,  while  in  a  non-stressed  'situation 
all  users  might  have  access  to  ail  channels,  under  stress 
priority  class  #1  might  be  assigned  3 0 - 1 00 %  of  the  channel 


capacity  while  a  loweT  priority  class  might  have  access  to  only 
say  40*  o£  the  channel  resources.  The  question  is  then  how 
does  one  make  the  assignment  to  keep  some  measure  of  overall 
performance  at  a  satisfactory  level.  In  Section  2  we  show  a 
method  for  making  the  allocation  so  that  high  priority  users 
maintain  their  GOS  under  jamming  while  lower  priority  users 
still  have  reasonable,  albeit  reduced  access  to  communications. 
Specific  examples  illustrate  the  method.  We  believe  the 
method  can  be  made  adaptive  in  conjunction  with  the  channel 
monitoring  scheme  proposed  in  Chapter  2. 

Chapter  4  -  Susceptibility  Analysis  Of  A  Frame  Synchronication 
Technique  For  Multichannel  Frequency-Hopped  BNCFSK 
Systems  In  Partial-Band  Noise  Ot  CW-Tone  Jamming 

Perhaps  the  most  vulnerable  aspect  of  a  multichannel  time 
division  multiplexed  digital  signal  is  the  mechanism  and  signal 
format  used  for  frame  synchronization.  If  frame  synchroniza¬ 
tion  is  lost  or  otherwise  interfered  with  then  the  bit  stream 
becomes  useless  as  an  information  signal.  In  this  Chapter  we 
investigate  the  effect  of  various  partial  band  jamming  signals 
on  the  frame  synchronization.  The  analysis  is  presented  in  a 
general  form  in  order  to  permit  parametric  studies  for  optimiz¬ 
ing  the  frame  synchronization  anti- jam  technique.  Several 
numerical  examples  are  worked  out  of  the  given  formulas  in 
order  to  show  the  effect  of  noise  jamming  without  frequency 
hopping,  noise  jamming  and  partial  hard  jamming  with  frequency 
hopping  and  various  signal  to  jam  ratios  on  the  probability  of 
correct  frame  synchronization.  It  would  appear  that  perhaps 
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greater  protection  ought  to  be  given  to  the  frame  synch  signal 
than  has  hitherto  been  the  case.  Several  coding  schemes  are 
possible  candidates. 


Chapter  5  -  Low  Rate  Error  Correcting  Codes 


As  An  Alternative  To  Conventional  Spread  Spectrum 


Since  direct  sequence  spread  spectrum  derives  its  processing 
gain  by  integrating  over  many  DS  chips  per  data  bit,  the  modulated 
DS  signal  can  be  viewed  as  a  law  rate  code.  The  DS  modulation 
however  is  not  exploited  as  a  code  since  none  of  its  structural 
properties  (except  autocorrelation)  is  exploited.  Indeed  DS 
pseudo-noise  can  be  viewed  as  an  "inner"  channel  operation  which 
is  undone  by  the  receiver  correlation  and  quite  independent  or  any 
data  that  is  sent.  As  such  it  is  an  "applique"  and  we  should  be 
able  to  do  better,  at  the  same  information  rate;  by  coding.  In 
this  Chapter  we  propose  a  particular  class  of  low  rate  linear 
block  codes  because  1)  they  have  the  best  distance  properties 
for  their  length  and  2)  they  are  amenable  to  effective  simple 
majority  logic  decoding.  We  show  the  construction  of  both  encoder 
and  decoder  and  present  an  analysis  of  the  performance  of  these 
codes.  Because  of  the  code  structure,  we  can  not  only  correct 
errors  and  thereby  improve  the  bit  error  rate  (3ER)  for  a  given 
spread  spectrum  over  conventional  methods  (in  the  range  of 
interesting  H^/Nq)  but  we  can  also  detect  errors  which  allow  us 
either  to  retransmit  or  to  delete  blocks  with  many  more  errors 
than  we  can  correct;  There  are  many  situations  especially  in 
the  transmission  of  digiticed  speech  when  it  is  preferable  to 
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delete  error  blocks  than  accept  them.  Thus  the  several  perfor¬ 
mance  measures  are  analysed.  They  are  bit  error  rate  (3ER) , 
undetected  bit  error  rate  (USER),  and  data  bit  deletion  rate 
(DBDR).  The  analysis  prompted  numerical  investigations  by 
computer  and  asymptotic  analysis  to  determine  the  precise  range 
where  coding  gain  (over  spread  spectrum)  is  achieved.  In  the 
range  of  2^/No  of  10  d3  or  greater  for  hard  decision  decoding, 
significant  coding  gain  is  achieved.  The  BER  can  be  reduced  by 
orders  of  magnitude  below  that  achievable  by  straight  spread 
spectrum.  The  U3ER  and  DBDR  are  even  more  impressive.  These 
codes  also  lend  themselves  to  soft  decision  decoding  when  coding 
gains  of  10  d3  or  more  are  theoretically  possible  down  to  3  dB. 

Further  in  this  chapter,  we  propose  a  new  class  of  very  low 

rate  convolutional  codes  which  we  call  "Complete  Convolutional 

Codes.”  These  codes  are  characterized  by  the  property  than  for 

k  - 1 

constraint  length  K  the  generators  consist  of  all  2  binary 
polynomials  (all  tap  connections)  and  hence  form  ta  vector  space. 
These  codes  have  demons tratable  good  performance  and  coding  gain 
since  we  show  that  the  free  distance  of  the  codes  can  be  explicitly 
derived  and  an  upper  bound  on  3ER  can  be  computed  from  this  free 
distance . 

Significant  advantages  can  be  had  by  using  coding  as  an 
alternative  or  in  conjunction  with  spread  spectrum.  These  advantage 
are  .summarized  in  the  chapter. 


Chapter  6  -  Adaptive  Null  Steering  Antenna  Arravs 

Perhaps  the  most  promising  of  ECCM  techniques  available  in 
a  bandwidth  limited  environment  is  that  of  null  steering  antennas. 

It  may  be  possible  to  achieve  40  or  50  d3  of  jamming  margin  with 
hardly  any  substantial  bandwidth  expenditure.  Furthermore, 
tactical  multichannel  radio  communications  usually  exist  on 
stationary  sites  and  do  not  have  to  operate  on  rapid  angular 
dynamics  (as  in  a  plane]  so  that  perhaps  costs  can  be  manageable. 

The  problem  in  adaptive  arrays  however,  is  that  they  may  be 
vulnerable  to  "blinking”  jammers  matched  in  period  to  the  natural 
modes  of  the  array  loops.  In  this  section,  a  basic  analysis, 
starting  with  first  principles  is  made  for  adaptive  antenna 
arrays.  Both  the  linear  mean  square  (LMS]  or  Widrow  algorithm 
and  the  maximum  Signal  to  noise  (SNR]  or  Applebaum  algorithm  are 
considered.  The  important  result  is  that  we  can  determine  the 
exact  deterministic  transient  solution  and  periodic  solution  of 
the  differential  equations  of  the  array.  We  believe  that  this 
has  never  been  done  previously  and  that  this  will  have  an  enormous 
impact  on  the  understanding  of  the  vulnerabilities  and  design  of 
adaptive  arrays.  In  particular,  it  should  shed  light  on  the  effects 
of  "blinking"  jammers  (a  subject  which  has  recently  received  much 
attention  by  experimenters].  Furthermore,  the  ability  to  solve 
the  equations  explicitly  for  deterministic  signals  allows  us 
fortuitously  to  obtain  the  analytic  behavior  of  the  adaptive  array 
under  stochastic  or  noisy  jamming  signals.  The  pursuit  of  this 
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approach  should  give  us  a  fairly  complete  picture  of  the  effects 
of  various  jammers  on  adaptive  arrays  and  the  means  to  make  them 
more  effective. 


Chapter  7  -  Modulation  and  Spreading  Techniques 

This  chapter  presents  the  results  of  an  extensive  effort  to 
discover  the  optimum  multialphabet  modulation  scheme  to  be  used 
with  spread  spectrum  that  will  result  in  the  ’’best”  combination 
modulation/ spreading  with  respect  to  a  jammer  under  the  constraint 
of  limited  total  bandwidth.  Various  modulation  schemes  were 
considered  including  BPSK,  QPSK,  16-ary  QASK ,  ?SK  and  the  effects 
that  a  number  of  jammers  would  have  on  them.  It  appears  that 
QPSK  for  direct  sequence  on  how  coherent  FSK  with  coding  in  a 
frequency  hopped  system  is  best.  The  curves  in  this  chapter  can 
be  used  to  demonstrate  the  trade-off  possible  by  using  adaptive 
and  dynamic  channel  allocation  as  discussed  in  chapter  3. 

Chanter  8  -  Conclusions  and  Recommendations 

This  very  brief  chapter  is  included  to  provide  an  overview 
of  which  has  been  accomplished,  how  and  whether  some  of  the  ideas 
that  evolved  during  the  course  of  this  work  ought  to  be  implemented, 
and  some  recommendations  for  further  work  both  in  analysis  and  in 
development . 


Appendix  A  -  Spread  Spectrum  Background 
Frequency  Hopping 


In  this  appendix  we  provide  an  up  to  date  and  thorough 
review  of  frequency  hopping.  This  appendix  together  with  those 
of  a  companion  report  by  S  consultants  are  a  self-contained 
exposition  of  spread  spectrum,  its  important  features,  problems 
and  limitations.  Several  of  the  topics  in  this  appendix  have  no 
appeared  in  the  literature  before.  The  treatment  here  includes 
the  properties  of  frequency  hopping,  correlation,  effects  of 
partial  band  jamming,  the  use  of  coding,  acquisition,  digital 
frequency  synthesiser  and  a  comparitive  analysis  of  frequency 
hopping  and  direct  sequences  modulation. 


CHAPTER  2 


Sit  Error  Rate  Monitoring  Using  Extreme-Value  Theory 

The  general  technique  and  rationale  for  using  extreme-value 
theory  (EVT)  as  the  basis  for  a  bit  error  rate  CBER]  performance 
monitor  is  described  in  [2.1].  An  appendix  is  included  at  the 
end  of  this  section  which  contains  extracts  from  [2.1]  and  hence 
it  provides  the  basis  for  the  remainder  of  this  section.  The 
specific  manner  in  which  it  will  be  used  in  the  current  study 
will  be  described  below. 

As  can  be  seen  from  either  [2.1]  or  the  Appendix,  each  of 
the  two  SER  estimates  (i.e.one  which  estimates  the  number  of 
♦  l's  in  error  and  one  which  estimates  the  number  of  -l's  in 
erTor]  depends  upon  two  parameters.  One  parameter  Clabeled  u] 
is  a  function  of  both  the  location  of  the  initial  distribution 
as  well  as  the  spread  ot  the  standard  deviation  of  the  distri¬ 
bution  of  the  underlying  samples,  while  the  second  parameter 
(labeled  a]  is  inversely  proportional  to  the  standard  deviation 
of  the  initial  distribution.  When  estimating  the  number  of 
-l's  in  error  one  is  interested  in  the  right  hand  tail  of 
the  density  and  hence,  in  the  notation  of  [2.1]-[].5],  in  the 
parameters  un  and  an.  Analogously,  the  parameters  and 
are  appropriate  for  estimating  the  left-hand  tail  of  a  density 
and  hence  are  used  to  estimate  the  number  of  ♦l’s  in  error. 

If  it  is  assumed  that  we  are  dealing  with  a  coherent 
QPSK  system  and  that  we  are  estimating  the  BER  in  one  of  the 
two  biphase  channels,  then  for  equally  likely  input  data,  the 


number  or  +l's  in  error  should,  on  the  average,  equal  the  number 


of  *  1 ' s  in  error.  Also,  from  the  defining  equations  for  u^, 
u^,  3^  in  the  Appendix,  it  can  be  seen  that  for  this1  situation, 


C2.ll 


and 

=n  *  al  C2‘ 

Using  these  two  equations,  it  is  reasonably  straightforward  to 

determine  when  a  jammer  is  present.  Consider,  for  example,  a  tone 
jammer,  the  effect  of  which,  is  to  cause  a  shift  in  the  probability 
density  of  the  final  test  statistic.  Figure  2.1  illustrates  a 
typical  situation.  Figure  2.1a  shows  the  density  functions  when 
the  jammer  is  not  present.  Figure  2.1b  shows  the  effect  of  the 
jammer  assuming  the  shift  Cdue  to  the  jammer)  is  to  the  right. 

It  is  clear  that  many  more  decisions  will  be  made  than  -1 
decisions,  hence  the  existence  of  an  outside  disturbance  !>•?• 
a  jammer)  has  been  confirmed.  Alternately,  u1>>|unJ,  thus  vio¬ 
lating  Eq.  .(2.1) • 

If  instead  of  being  a  tone  j  ammer  the  interference  had  been 
a  noise  jammer,  the  situation  would  look  as  shown  in  Figure  2.1c. 
Mow  the  number  of  ♦l  decisions  is  still  equal  to  the  number  of 
-1  decisions,  but  because  the  variance  of  either  density  is  much 
larger  than  what  it  initially  was,  the  values  of  or  an  will  be 
much  smaller,  thereby  indicating  the  presence  of  an  outside  dis¬ 
turbance.  This  latter  statement  assumes  that  reasonable  esti¬ 
mates  of  a,  and  a  ,  say  d.  and  d  .  respectively,  can  be 
i,  n  i  n 

obtained  when  the  jammer  is  present  and 
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before  any  counter-measures  have  been  taken  to  combat  the  jammer. 
This  will  only  be  the  case  if  a  learning  period  is  provided 
during  which  time  the  receiver  knows  whether  a  given  sample 
is  a  si. 

If  this  learning  period  is  not  provided,  one  can  still 
get  an  indication  of  the  jammer’s  presence  in  the  following 
manner:  As  described  in[l] ,  the  parameters  c^,  a  ,  u^,  and  un 
are  all  estimated  using  the  extremes  (either  maxima  or  minima] 
of  sets  of  n  samples.  When  the  learning  period  is  absent,  some 
of  those  extremes  will  be  incorrect,  and  in  particular,  when 
the  noise  jammer  is  present,  most  of  these  samples  will  be 
incorrect.  However,  they  will  be  incorrect  in  a  very  specific 
manner,  namely  they  will  all  be  very  close  to  the  threshold 
used  to  distinguish  between  il's  (zero  in  this  case).  This 
will  result  in  a  very  small  standard  deviation  and  hence  a 
very  large  value  of  S.  Hence  even  though  an  incorrect  value 
of  a  is  being  obtained,  it  will  nonetheless  be  different 
enough  from  the  value  of  d  obtained  when  the  jammer  was.  absent 
to  clearly  indicate  the  jammer's  presence. 

Having  determined  the  fact  that  the  system  is  being  jammed, 
it  remains  to  indicate  what  response  to  the  jammer  the  system 
should  initiate.  Ideally,  one  would  like  to  use  just  enough 
counter-measure  (e.g.  spectrum  spreading]  to  bring  the  system 
back  to  where  it  was  (in  a  BER  sense]  before  it  was  jammed. 

This,  however,  requires  precise  knowledge  of  the  BER  while 
the  system  is  being  jammed  and  before  any  counter-measure  is 
initiated.  Because  EVT  is  only  accurate  on  the  tail  of  a 
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density,  and  because  the  system  at  this  point  is  no  longer 
operating  on  the  tail,  accurate  BER  estimation  using  EV7  cannot 
take  place.  Therefore,  a  compromise  is  to  spread  the  spectrum 
some  predetermined  amount,  say  SOI  of  the  total  possible  spread, 
and  test  to  see  if  this  puts  the  system  back  on  the  tails  Ce*g- 
test  to  see  if  the  estimated  parameters  0^,  S1 ,  and  fin 
return  to  their  pre-jamming  values).  If  it  does,  operate  the 
system  as  it  is,  if  it  does  not,. spread  the  spectrum  the  remain¬ 
ing  SOI.  Assuming  the  system  has  at  its  disposal  enough  pro¬ 
cessing  gain  to  nullify  the  jammer,  the  system  will  at  one  of 
the  two  steps  above  return  to  operating  at  a  low  BER,  and  accur¬ 
ate  BER  estimation  can  then  once  again  take  place. 

To  verify  the  above  procedure  essentially  requires  a  com¬ 
puter  simulation.  However,  some  indication  of  how  it  will 
perform  can  be  obtained  by  assuming  a  specific  system  and  using 
the  true  parameter  values  whenever  they  are  needed.  Therefore 
assume  we  are  looking  at  the  inphase  channel  of  a  QPSK  system 
Cor  equivalently  a  simple  3PSK  system)  consisting  of  a  coherent 
demodulator  followed  by  an  integrate  and  dump  filter.  The 
signal  is  received  in  the  presence  of  additive  white  gaussian 
noise  (AWGN)  of  two-sided  noise  spectral  density  as  well 
as  the  Cpotsatial)  presence  of  a  gaussian  noise  jammer  occupying 
a  bandwidth  equal  to  that  of  the  main  lobe  of  the  signal 
spectrum. 

There  are  two  considerations  that  must  be  addressed  in 
evaluating  the  error-rate  monitoring  procedure  for  this  situa¬ 
tion.  Since,  as  was  indicated  above,  for  the  noise  jammer  the 
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key  parameters  are  and  a^,  it  has  to  be  ascertained  that 
these  values  do  indeed  undergo  a  significant  enough  change  in 
value  when  a  jammer  is  present  (from  what  they  were  when  the 
jammer  was  absent},  and  that  given  a  siraole  enough  change  in 
value,  this  change  can  be  estimated  from  samples  of  the  detected 
waveform. 

Consider  the  first  question.  Because  both  the  noise  and 

interference  are  gaussian,  the  output  statistic  of  the  receiver 

has  a  symmetrical  probability  density.  Under  this  condition,  it 

is  straightforward  to  show  that  a,  *a  .  Therefore  onlv  values 

1  n 

of  a^  need  be  considered.  For  any  density,  is  defined  by 

ax  *  n£Cux)  (2.3} 

where  u^  is  the  solution  of  the  following  equation  (see  the 
Appendix} : 

FTu1)  -  i  (2.4) 


In  (2.3)  and  (2.4),  F(*)  and  f(*)  are  the  cumulative  distribu¬ 
tion  function  and  probability  density  function,  respectively,  of 
the  output  statistic  and  the  significance  of  n  is  explained  in 
[2.1]-[2.5]  and  will  be  seen  below.  If  the  signal  component 
of  the  output  statistic  for  the  problem  at  hand  is  normalized 
to  unity,  it  can  be  shown  in  a  straightforward  manner  that  the 
following  equation  holds: 
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where  E  is  the  signal  energy,  S  the  signal  power,  J  the  jammer 


power,  2£c  the  bandwidth  in  Hertt  of  the  mainlobe  of  the  signal 
spectrum,  T  the  duration  in  seconds  of  one  data  symbol,  and 


Examule 


* (x)  £  _1_  f 
/iir 


*  '  dy 


(2.6) 


Assume  now  that  the  system  is  operating  with  an  E/No  of 

10  d3  and  without  any  spread  spectrum  processing  gain.  With 

no  jammer  present,  can  be  shown  to  equal  14.6  when  n*10J. 

If  a  jammer  with  30  dB  more  power  than  the  signal  (i • e '•  ^  ■  30  dB) 

2 

in  a  bandwidth  of  2fc  *  if  suddenly  starts  to  interfere  with 
signal,  the  value  of  a^  reduces  drastically  to  1.7  x  10 

If  the  system  uses  this  reduction  in  a^  as  an  indication 
of  the  jammer's  presence  and  hence  spreads  its  spectrum  by  a 
factor  of,  say,  500  (i.e.  2fcT  *  1000),  a^  will  then  increase 
to  4.4.  This  clearly  indicates  that  while  significant  improve¬ 
ment  in  system  performance  has  been  obtained  by  spreading  the 
spectrum,  the  spread  was  not  sufficient  to  bring  the  system 
back  to  the  point  it  was  before  the  jamming  signal  was  received. 

In  any  event,  it  is  seen  that  a^  does  indeed  undergo  a  signif¬ 
icant  enough  change  in  value. 

Relative  to  the  second  question,  suppose  we  have  a  learning 
period  available  and  that  a^  is  to  be  estimated  by  defined 

as  follows: 
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where 


ana 


N' 


2  [xt 

i*l  min 


(2, 


is  the  estimate  of  the  standard  deviation  or  the  minimum  of  a 

samples  of  the  output  statistic  and  is  obtained  by  computing 

the  sample  standard  deviation  of  N  samples,  each  one  being  the 

smallest  from  a  group  of  n  independent  samples  of  the  output 

of  the  receiver.  The  rationale  for  the  estimate  of  [2.6)  is 

described  in  [2.2]  and  [2.3].  Each  of  the  X.  in  [2.7) 

xmin 

is  the  minimum  of  n  indeoendent  samoles,  and  X  .  is  the  samole 

*  mm 

mean  of  the  N  minima. 

Obtaining  the  exact  probability  density  of  3,  does  not 
appear  to  be  tractable.  Therefore,  an  approximation  which 
becomes  more  and  more  valid  is  N  increases  will  be  used.  In 
[2.2]  and  [2.3],  it  is  shown  that  (3^)  "  is  asymptotically  nor¬ 
mal  with  mean  equal  to  l/a^  and  variance  equal  to  l.l/N^".  There¬ 
fore,  using  this  asymptotic  result,  we  obtain  a  confidence  inter¬ 
val  on  5,  as  follows: 


P  jacj 


) 


where  a  and  b  are  constants  and  $  [•)  is  defined  in  (2.6). 

If,  for  example,  we  are  interested  in  the  probability  of 
2^  being  within  a  factor  of  wo  of  its  true  value  when  N'*23, 
then 


(2.3) 
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Therefore,  notwithstanding  the  fact  that  we  are  using  an 
asymptotic  result  for  finite  N  (e.g.  N  ■  20),  it  appears  that 
one  can  indeed  obtain  an  accurate  enough  estimate  of  to  deter 
mine  when  the  system  is  being  jammed. 

Conclusions 

From  the  above  results,  it  is  clear  that  EVT  can  serve  both 
to  indicate  the  presence  of  a  jammer  as  well  as  to  estimate  the 

average  probability  of  erroT  of  the  system  once  a  sufficient 

/ 

degree  of  processing  gain  is  given  the  system  to  combat  the 
jammer.  Furthermore,  the  actual  implementation  of  such  an 
error  rate  monitor  is  relatively  straightforward,  requiring  only 
that  extreme  samples  taken  from  sets  of  samples  of  the  final  re¬ 
ceiver  test  statistic  be  accumulated  and  then  used  to  compute 
the  estimates  of  the  u’s  and  the  a’s  (e.g.  Eq.  (2.7)).  From 
these  estimates,  the  desired  estimates  of  probability  of  error 
as  given  by  rqs.  (A.l)  and  (A. 2)  of  the  Appendix  are  simply 
obtained. 
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Appendix  to  Chapter  2 
Review  of  Extreme-Value  Theory 

Given,  a  collection  of  n  iid  random  variables,  there  are 
certain  conditions  on  the  distribution  function  of  the  random 
variables  such  that  the  distribution  of  the  largest  random 
variable  from  the  set  of  n  can  be  expressed  in  an  asymptotic 
form  (see,  e.g.  (2 . 5 ] - [2 . S ] ) .  Specifically,  Gumbei  in  [2.3] 
shows  there  are  three  distinct  asymptotic  forms.  Of  the  three, 
the  so-called  "exponential  type"  appears  to  be  of  most  interest. 
This  limiting  form  (for  maximum  values}  applies  for  any  initial 
distribution  which  goes  to  tero  on  its  right-hand  tail  in  a 
manner  similar  to  an  exponential.  In  particular,  it  is  shown 
in  the  above  references  that,  with  F(x)  and  f(x)  the  initial 
distribution  and  density  functions,  respectively,  if 


lim  f  fx)  _  _  _  lim  f  ‘  (x 


■ray 


x- 


r x ) 
uT 


then  the  distribution  function  of  the  maximum  value  is  given  by 


where 


and 


F(un)  -  1  -  i 


'  aiCV 


:o 


The  analogous  condition  for  the  left-hand  tail  (i.e.. 


minimum  value)  is  as  follows.  If 


then 


lim  ffx)  .  lia  |lQp 

X*  -  F (X)  X--«  f [X) 


1-exp 


-exp[a1 (x-ux)] 

m  m 


where 

F(Ui)  -  i 

and  ■  nfCu^) 

Since,  in  general,  the  underlying  initial  distribution  is 
unknown,  the  parameters  u^,  an,  u^,  and  have  to  be  estimated 
from  sample  values.  This  usually  proceeds  as  follows.  A  total 
of  K*nN  samples  from  a  given  distribution  are  taken  and  divided 
into  N  groups  of  n  samples  per  group.  If  it  is  desired  to  esti¬ 
mate  u^  and  an>  in  each  of  the  N  groups  the  largest  of  the  n 
samples  is  chosen.  From  these  N  largest  values,  it  can  be 
shown  that  a  variety  of  procedures  exist  from  which  un  and  can 
be  estimated.  A  method  for  obtaining  maximum- likelihood  estimates 
of  un  and  a^  is  described  in  (2.3)  along  with  various  other 
estimation  techniques. 

It  is  desired  to  estimate  u.,  and  a^,  a  similar  procedure 
using  minimum  instead  of  maximum  values  can  be  performed.  Finallv, 
appropriate  estimates  for  the  two  types  of  error  in  a  binary 


digital  communication  system  are  given  by 


pfa  •  CA-l) 

and 

Pfd  *  5«3cpCa1Cxt- u^]  CA.2) 

where  xt  is  the  threshold  that  the  system  is  operating  with.  (e.g. 
tero  for  a  3PSX  system)  and  where  the  two  probatilities  have 
been  labeled  ?^a  and  ?£d  in  analogy  with  classical  radar  notation, 
to  represent  the  false  alarm  and  false  detection  probability 
respectively. 
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CHAPTER  3 


Dynamic  Allocation  of  Channel  Capacity 

In  the  Interim  Report  for  the  period  September  30,  1978 
to  January  31,  1979,  a  specific  technique  for  the  dynamic  allo¬ 
cation  of  channel  capacity  was  presented.  This  technique  made 
use  of  an  equivalent  traffic  model  in  which  users  from  different 
priority  classes  shared  a  certain  percentage  of  the  total 
resources  available,  and  were  then  credited  with  a  certain 
’’effective''  traffic  which  was  in  general  greater  than  their 
actual  traffic.  This  increase  in  traffic  was  of  course  due  to 
the  fact  that  the  channels  assigned  to  a  given  user  class,  say 
the  irl1  class,  were  in  general  available  to  other  useT  classes 
and  hence  affected  the  grade  of  service  (GOS)  (i.e.  the  proba¬ 
bility  that  an  incoming  call  will  be  blocked)  cf  a  user  in  the 
ith  class. 

However,  further  analysis  of  that  model  demonstrated  that 
it  had  a  significant  flaw  to  it.  This  can  be  demonstrated  as 
follows:  Assume  a  channel  allocation  between  two  user  classes 
is  as  shown  in  either  Table  3.1  or  Table  3.2  below. 


Table  3.1 
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able  3.2 


N.  channel  ? 

user  *  >v  1  2  3  .  .  .  .n^  n^+1  ni't‘“ . C 

1  X  X  X  ...  .0  X  X . 0 

2  o  o  o  ....  o  x  x . x 


In  both,  allocation  schemes,  there  are  a  total  of  C  channels. 

In  the  first  scheme,  user  class  #1  has  n^  channels  fully  dedicated. 

to  it,  while  in  the  second  scheme,  user  class  *1  has  access  to 

n^+l  channels,  n^  on  a  fully  dedicated  basis,  and  one  channel 
th 

Cthe  n^+1  * )  on  a  shared  basis  with  user  class  *2.  For  a  given 
offered  traffic  of  user  class  #1,  say  e^  erlangs,  it  is  physically 
obvious  that  the  GOS^  of  user  class  #1  with  the  second  allocation 
scheme  has  to  be  better  than  the  GOS^  with  the  first  allocation 
scheme.  Yet  because  the  GOS^  of  user  class  #1  was  computed  in 
the  Interim  Report  based  upon  the  "effective”  traffic,  which  was 
a  linear  combination  of  the  actual  traffic  of  user  classes 
*1  and  *2,  it  is  clear  that  for  some  value  of  e7  (and  all  higher 
values),  the  actual  (i.e.  offered)  traffic  of  user  class  *2, 
the  GOS^  of  the  second  allocation  scheme  will  appear  to  be 
poorer  than  the  GOS1  of  the  first  allocation  scheme. 

In  other  words,  the  flaw  with  the  traffic  model  presented 
in  the  Interim  Report  was  the  mdnner  in  which  the  shared  traffic 
was  combined  to  yield  the  "effective"  traffic,  this  latter 
quantity  being  the  basis  for  the  GOS  calculations.  Therefore, 
for  this  present  report,  it  was  necessary  to  use  a  different 
traffic  model,  and  such  a  model  will  be  presented  below.  While 
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the  original  motivation  for  changing  to  this  new  model  was 
to  avoid  the  flaw  just  described  in  the  old  model,  it  turned 
out  that  this  new  model  had  the  further  advantage  of  being  easier 
to  use.  In  particular,  whereas  the  first  model  required  the 
computation  of  the  pseudo-inverse  of  an  NxN  matrix,  N  being  the 
number  of  priority  classes,  no  such  calculation  is  necessary 
with  the  current  technique.  In  what  follows,  the  model  will 
first  be  described,  and  then  numerical  results  illustrating  its 
behavior  will  be  presented. 

Assume  that  there  are  three  basic  priority  classes  with  n^ 
users  in  the  i  1  class,  i-1, 2, 3,  and  let 

n  £  I  nt  (3.1) 

i-1  1 

be  the  total  number  of  users.  Let  e^  be  the  number  of  erlangs 
of  traffic  of  the  ir^  priority  class,  and  let 

3 

«  £  2  e-  (3.2) 

i-1 

While  there  are  many  types  of  priority  assignments  that 
could  be  made,  the  following  one  will  arbitrarily  be  considered 
in  the  remainder  of  this  section:  User  class  *1  will  have  highest 
priority,  followed  by  *2  and  then  *3.  All  users  from  all  classes 
will  have  access  to  every  channel.  However,  if  a  user  from  class 
*1  is  blocked  when  attempting  to  access  one  of  the  channels, 
that  user  will  be  allowed  to  "bump"  one  of  the  *3  users  (pro¬ 
viding  at  least  one  of  the  channels  is  in  fact  occupied  by  a 
*3  user)  and  hence  gain  immediate  access  to  the  channel. 

Given  the  above  model,  the  grade-of-service  (GCS)  for  a  user 
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in.  each  or  the  three  priority  classes  will  be  computed.  Assume 
there  are  a  total  of  c  channels,  where  in  general  c<n  and  in 
most  instances  c<<n,  and  denote  by  ?(c,e)  the  probability  of 
blocking  in  a  system  with  c  channels  and  e  erlangs  of  traffic 
as  computed,  from  the  so-called  Erlang- B  model. 

Consider  user  class  *1.  A  call  will  be  blocked  in  this 
case  only  if  all  c  channels  are  occupied  by  either  other  ?1 
users  or  *2  users.  To  compute  this  probability,  consider  the 
probability  of  one  channel,  say  the  channel,  being  occupied 
by  a  user  from  class  *i.  Denoting  this  latter  probability 
by  P^,  we  have 

P^  »  P [channel  j  occupied  by  a  user  form  *i  |  channel  j  is 
occupied]  •  P [channel  j  is  occupied]. 

The  conditional  probability  on  the  rhs  will  be  taken  to  be  n-/n. 

The  second  term  on  the  rhs  is  computed  as  follows  : 

c 

P [channel  j  is  occupied]  *  Z  ? [channel  j  is  occupied! 

k»o 

exactly  k  channels  occupied]  •  P [exactly  k  channels  occupied] 


Therefore 


?i 
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Using 


P  l  i » j  >  ^ ^ ) 


3),  there  is  straightforward  to  show  that 
?  {i  channels  used  by  #1,  j  channels  used  by  #2, 

k  channels  used  by  #3 
i  channels  not  being  used} 


TTrfbr  VVV'1' Wp3> 


(3.4) 


where  i+j+k+1 
is  given  by 


wnere 


c.  Finally,  GOS^,  the  GOS  of  a  user  in  class  #1 


G0S1  -  P(c,e1>e2)Pc 


A  ^  ^  C 1 

pc  ■  1  1  T^T r  ?(i»j  ,o,o) 

i»o  j*o  ‘ J  * 

i*j*c 


(3.5) 


(3.o) 


and  where  PCc^e^+e,)  is  the  Erlang  B  blocking  formula  for  c 
channels  and  e2+e2 .erlangs  of  ordered  traffic. 

For  a  user  in  class  #2,  the  GOS  is  simply 

GOS2  »  P(c,e1^e2*e3)  •  P(c,e)  (3.7) 


Finally,  for  a  user  in  class  #3,  the  definition  of  GOS  has  to  be 
generalited  somewhat.  Since  a  #3  user  can  be  "bumped"  from  a 
channel  after  he  has  started  sending  a  message  by  a  #1  user, 
the  GOS  of  a  *3  use  will  be  defined  as  the  probability  of  either 
being  blocked  or  of  being  "bumped."  Since  these  are  mutually 
exclusive  events,  we  have 

GCS3  *  ?  (block)-  *  ?(bump)  (3.8; 
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The  last  term  on  the  rhs  of  (3.3)  is  the  probability  that  a 
user  from  class  *1  wants  to  access  a  channel  during  the  tine  that 
a  user  from  class  #3  is  transmitting  and  that  no  free  channel 
is  available.  To  upper  bound  this  probability,  assume  a  *3 
message  uses  one  channel  during  the  interval  (o,t),  where  t 
is  a  random  variable  with  an  exponential  probability  density. 

The  probability  that  at  least  one  user  from  class  #1  wants  to 
access  the  channel  during  that  t  second  interval  is  given  by 


P{at  least  one  *1  arrives  in  (o,t)jt>  *  1-e"  l1 


(3.9) 


where  a,  is  the  average  number  of  arrivals/sec  for  user  class  *1. 
Eq.  (3.9)  followed  because  of  the  assumed  Poisson  arrival  density 
for  any  of  the  users.  To  obtain  the  unconditional  probability 
that  a  *1  user  arrives  while  the  *3  is  transmitting  it  is  necessary 
to  average  (3.9)  over  the  density  of  t.  Therefore 

P(at  least  one  *1  arrives  while  the  *3  user  is  on} 


Cl-e  *  ]± 


Alc,l  .-t/u.  ,  u 

Jr  *  dt  -  1  -i — 


u  ,vl 


1  ♦  a  ^  y 


(3.19) 


where  1/u  is  the  average  duration  of  a  message,  assumed  the 
same  for  any  user  in  any  priority  class.  Using  (3.10),  we  have 
P(bump)  «  ?{(the  *3  user  is  not  initially  blocked), 

(a  *1  user  arrives  while  the  #3  user  is  on) , 

(all  channels  are  occupied  when  the  *1  user  arrives)  , 
(of  all  the  *3  users  using  channels  when  the  *l  user 
arrives,  this  particular  ?3  is  the  one  bumped); 
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<  ?{user  #5  not  initially  blocked,  a  *1  user  arrives  while  *3 
user  is  on} 

el 

-  [1-PCc, e)]^~ 

*  el 

Therefore,  combining  (3.8)  and  (3.11)  yields 

G0S3  <  P(c,e)  ♦  [1*P (c ,e)  ]  y-f- 


(3.11 


With  the  above  equations,  the  performance  of  the  system  can 
be  computed.  Notice  that  the  procedure  itself  is  independent  of 
whether  or  not  the  system  is  being  jammed.  The  presence  or 
absence  of  a  jammer  will  determine  how  many  channels  are  availa¬ 
ble  and  hence  what  the  GOS  will  be.  Therefore  there  is  a  direct 
coupling  between  the  performance  monitoring  procedure  described 
in  Chapter  2  and  the  capacity  allocation  scheme  currently  being 
discussed.  A  simple  flow  chart  is  shown  in  Figure  3.1.  The 
boxes  labeled  "spread  3W  by  SOI"  control  the  number  of  channels 
that  are  available.  For  example,  if  the  only  means  at  the  dis¬ 
posal  of  the  system  at  a  given  point  in  time  to,  say  halve 
the  net  information  bandwidth  (and  hence  allow  this  now  unused 
bandwidth  to  be  consumed  by  extra  spreading  of  the  spectrum)  is 
to  reduce  the  number  of  available  channels,  then  these  boxes 
would  initiate  a  command  to  reduce  the  number  of  available 
channels  by  SOI. 

Notice  that  both  the  priority  arrangements  and  the  specific 
channel  assignment  can  be  performed  by  the  same  "switch". 

Ail  that  is  necessary  is  for  the  switch  to  keep  information 
about  which  priority  class  users  are  occupying  which  channels. 
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Figure  3.1 

Flow  Chart  for  Jam-sensing  System  Control 


In  that  way  the  switch  will  be  able  to  bump  a  #3  priority  class 
user  when  necessary  (i.e.  when  all  channels  are  occupied  and 
a  #1  priority  class  user  requests  a  channel)  as  opposed  to 
possibly  bumping  a  user  from  either  classes  #1  or  #2  by  mistake. 
Finally,  if  at  a  given  point  in  time  when  it  is  necessary 
to  invoke  this  procedure  there  are  c-<co  *3  class  users  occupying 
channels,  the  choice  of  precisely  which  #3  to  bump  can  be  made 
in  a  variety  of  ways  (e.g.  randomly,  the  last  to  access  a 
channel,  etc. 

Equations  (3.5),  (3.7),  and  (3.11)  were  evaluated  numericallv 
for  several  combinations  of  traffic  parameters  and  the  results 
are  presented  in  Table  3.3.  The  assumption  implicit  in  the 
traffic  values  chosen  m  Table  3.3  is  that  the  highest  priority 
users  (i.e.,  users  in  class  #1)  only  transmit  a  small  amount  of 
traffic  whereas  the  lowest  priority  users  (i.e.,  class  #3)  have 
by  far  the  greatest  amount  of  traffic  to  transmit. 

Jpon  examining  the  results  of  Table  3.3,  it  can  be  seen  that, 
at  least  for  the  traific  values  chosen,  the  GOS  for  user  class 
*1  a2mos t  perfect  (i.e.  less  than  10  ^ )  except  when  there  are 
only  three  channels  and  even  in  that  case,  it  is  less  than  10*3. 

On  the  other  hand,  the  GOS  for  user  class  #3  is  uniformly  poor 
except  when  24  channels  are  used  in  combination  with  E1«l, 

02*2,  and  Ej«10  erlangs.  Both  of  these  results  are  direct  con¬ 
sequences  of  the  priority  assignment  (i.e.  the  ability  of  *1 
user  to  "bump”  a  *3  user). 

Also,  as  seen  from  Eq.  C3.‘),  the  GOS  of  a  user  in  class  *2 
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is  just  that  which  would  he  predicted  by  using  the  Erlang  3 
formula  in  which  the  total  number  of  channels  and  the  total  ' 
traffic  (i.e.  E*5^+E, *E-}  are  used.  It  is  seen  that  if  say, 

E-  is  IQ  erlangs  and  E,  is  2  erlangs,  GGS7  will  seriously  de¬ 
grade  .if  the  total  number  of  channels  C  is  reduced  from  24-  to  12. 
Therefore,  for  the  particular  parameter  values  used  in  Table  3.5, 
this  priority  scheme  only  gives  true  protection  in  terms  of  GCS 
to  a  user  in.  class  *1  when  the  overall  system  is  heing  jammed 
severely  enough  to  cause  a  significant  reduction  (e.g.  50%)  in 
the  number  of  channels  that  are  available  for  use. 

However,  all  the  users  in  all  priority  classes  that  gain 
access  to  the  channels  receive  the  added  protection  against  the 
jammer  that  the  increased  processing  gain  allows.  This  can  be 
illustrated  very  simply  by  an  example.  Suppose  there  are  initially 
24  channels  available  and  the  offered  traffic  of  each  of  the 
three  user  classes  is  as  follows: 


El  •  0.S 


e3  «  10 

Assume  the  system  is  operating  normally  (i.e.  it  is  not  being 
jammed)  and  no  priority  assignment  is  necessary.  From  the  Erlang 
3  formula,  the  GOS  of  any  user  will  be  1-.  3  x  1 3  *  J  .  Note  that 
this  applies  to  any  user  in  any  priority  class. 

If  a  jamming  signal  is  now  detected,  one  can  invoke  the 
priority  system.  If  this  is  done  without  changing  the  total  number 
of  channels  QL.t.  C  still  equals  24),  the  protection  against  the 
jammer  will  oe  whatever  it  would  have  been  without  the  priority 
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Table  3.3 


Grace  c£  Service  Dependence  on  Channel  Allocation 


GOS. 

o 


c 

r 

“I 

*2 

“3 

G0S1 

GOS2 

(upper  bound) 

24 

.  5 

2 

10 

.  2  x 

10'33 

.13  x  10*2 

.  33 

12 

.5 

n 

im 

10 

.  71  x 

10-12 

.22 

.48 

6 

.  5 

2 

10 

.75  x 

10'5 

.57 

.71 

3 

.  5 

2 

10 

.  34  x 

10*2 

.”8 

.85 

24 

.1 

10 

<10' 

■  37 

.87  x  10‘3 

.92  x  10'1 

12 

.1 

2 

10 

.  16  x 

10'14 

.20 

.28 

6 

.1 

2  • 

10 

.  5  x 

10-6 
-  •> 

.  56 

.60 

3 

.1 

2 

10 

.11  X 

10  “ 

.77 

.79 

24 

.  3 

2 

20 

.11  X 

10‘34 

.12 

.41 

12 

.  5 

2 

20 

.29  x 

10'14 

-6 

.51 

.67 

6 

.  5 

2 

20 

.31  x 

10  0 

.73 

.83 
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assignment.  However,  the  GOS  for  a  user  in  class  *1  will  de¬ 
crease  to  less  than  1Q*JJ  C i.e.  it  will  effectively  have  a 
dedicated  channel)  while  the  GOS  for  a  user  in  class  *3  will  become 

greater  than  0.33.  A  class  *Z  user  of  course  will  continue  to 

♦ 

have  a  GOS  of  1.3  x  IQ*3. 

If  the  jamming  signal  is  reasonably  strong,  this  will  not 
be  a  viable  procedure.  That  is,  some  additional  processing  gain 
will  be  needed  and  therefore  suppose  that  the  system  gees  to 
the  priority  scheme  and  in  addition  reduces  the  number  of  channels 
from  21  to  12.  This  will  yield  an  extra  3d3  in  processing  gain 
agains  the  jammer  and  for  all  practical  purposes  will  not  effect 
the  GOS  of  a  user  in  class  *1  Cit  will  now  be  less  than  L0~'“)  . 
However,  it  will  substantially  degrade  the  GOS  for  a  user  in 
class  #2.  Instead  of  being  1.3  x  10'-3,  it  will  now  be  0.22. 

Also,  GOS-  will  increase  from  0.33  to  0.13. 

V 

The  final  question  then  is  what  has  this  extra  3  d3  in 
processing  gain  bought  us?  The  answer  to  this  question  will  be 
deferred  until  Chapter  7  in  which  the  performance  of  various 
modulation/ spreading  techniques  under  different  jamming  conditions 
is  discussed. 

In  light  of  the  above  discussion,  it  can  be  seen  that  the  com¬ 
bination  of  the  priority  system  with  a  variable  number  of  channels 
which  is  determined  by  the  condition  of  the  channels  Ci.e. 
whether  or  not  they  are  being  jammed)  can  provide  the  necessary 
channel  accessibility  to  users  having  the  highest  prioritv  clus 
increased  jammer  protection  to  all  users  who  gain  access  to  a 
channel.  Finally,  the  specific  priority  scheme  chosen  here  was 
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Only  meant  to  be  indicative  of  the  type  of  strategy  one  might 
employ,  and  clearly  numerous  similar  schemes  (e.g.  a  #1  user 
bumps  both  * 2  and  #3  users,  while  a  #2  user  bumps  #3  users) 
would  be  considered. 

Conclusions 

In  summary,  it  has  been  demonstrated  that  a  simple  scheme 
to  dynamically  allocate  channel  capacity  on  a  priority  basis 
can,  when  used  in  conjunction  with  a  procedure  to  monitor  the 
condition  of  the  channels,  guarantee  access  to  the  channels  to 
the  high  priority  users,  plus  guarantee  improved  performance 
to  any  of  the  users  who  are  able  to  access  a  channel.  This 
improved  performance  arises  because  the  process  of  reducing 
the  total  number  of  channels  in  use  (while  still  providing  a 
good  GOS  to  high  priority  users)  allows  the  remaining  channels 
to  spread  their  bandwidths  to  the  greatest  extent  possible  and 
hence  combat  a  jamming  signal  much  more  successfully  then  they 
otherwise  would  have  been  able  to  do. 
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CHAPTER  4 


Susceptibility  Analysis  of  a  Frame  Synchronization 
Technique  for'  Multichannel  Frequency ^Hopped  3NCF5K 
Systems  in  Partial-Band  Noise  or  CW-Tone  Jamming 

4. 1  Introduction 

Jamming  can  affect  a  communication  system  in  several  ways. 

It  can  attack  the  information-bearing  signal  directly,  or  it 
can  reduce  system  effectiveness  through  attacks  on  phase  syn¬ 
chronization,  spreading-code  synchronization,  etc.  In  this  chap¬ 
ter  we  consider  the  vulnerability  to  jamming  of  a  frame  synchro¬ 
nization  technique  for  a  multichannel  frequency -hopped  binary 
noncoherent  frequency -shift -key  system.  We  investigate  the  cases 
of  partial-band  noise  or  CW-tone  jamming.  The  derivations  of 
the  key  results  will  be  presented  in  as  general  a  form  as  possible, 
in  order  to  permit  extensive  parametric  studies  and  optimizations. 

The  formulas  are  applicable  to  many  special  cases. 

We  let  Nc  represent  the  number  of  time-channels ,  with  the 
last  one  being  the  frame-sync  channel.  The  receiver  looks  for 
the  time-channel  which  exhibits  a  particular  frame-sync  pattern 
of  symbols  over  an  interval  of  several  frames.  For  the  first 
N£-l  time- channels ,  because  of  the  randomness  of  the  data  within 
them,  there  is  normally  a  very  small  probability  that  one  of 
them  would  be  mistaken  for  the  sync  channel.  Jamming  may  in¬ 
crease  this  probability,  and  also  may  increase  the  probability 
that  the  correct  channel  may  be  missed  during  the  search  for 
frame  synchronization. 
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It  is  of  interest  to  study  the  system  vulnerability  for 
non-spread  and  also  frequency- hopped  systems.  For  the  frequency 
hopped  cases >  we  let  W3  represent  the  total  number  of  frequency 
slots  available  for  use,  with  two  selected  (for  this  binary 
analysis)  at  any  one  time  as  "mark”  and  "space"  frequencies. 

In  a  partial-band  jamming  strategy,  the  jamming  waveform  pro¬ 
duces  energy  in  of  the  slots,  with  allowed  to  be, 
conceptually,  any  number  from  one  to  Ns.  We  assume  that  the 
system  hops  once  per  time-channel,  in  a  pseudo-random  manner 
which  is  unknown  to  the  jammer.  It  is  further  assumed  that 
the  receiver  is  synchronized  to  the  spreading  code. 

In  the  following  section  we  derive  certain  key  general 
results.  These  are  then  used  to  obtain  specific  results,  which 
are  then  presented  in  the  next  section. 

4.2  Derivation  of  General  Results 

For  the  results  obtained  here,  it  is  not  necessary  that 
we  consider  specific  jamming  cases.  These  follow  in  the  next 
section.  Let  us  assume  that  the  receiver  has  no  prior  knowledge 
of  the  correct  sync  channel,  and  that  it  begins  a  search  for  a 
frame-sync  pattern  in  the  i^*1  time  slot.  The  frame-sync  pattern 
is  one  of  alternating  "ones"  and  "zeros".  The  first  time-slot 
searched  can  therefore  be  any  one  of  the  N  ,  with  equal  proba- 
’oility. 

The  testing  for  a  particular  tine-slot  involves  checking 
the  n  consecutive  binary  symbols  in  that  slot  over  n  con- 

s  5 

secutive  frames.  If  at  least  b  of  the  n.  match  the  svnc 
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pattern,  then  "Frame  Sync"  is  declared  by  the  receiver.  Other¬ 
wise,  the  testing  ends  for  that  slot  and  is  begun  for  the  next 
slot.  This  continues  until  a  match  is  obtained.  The  results 
obtained  here  provide  the  means  for  optimal  choices  of  both  ns 
and  b. 

The  following  fundamental  probabilities  are  required: 

(a)  P^(FS)  «  probability  of  a  false  "Frame-Sync"  decision 

for  the  ith  time-slot 

(b)  P^(CD)  *  probability  of  a  correct  dismissal  decision 

for  the  ith  time-slot 


-  1  -  PtCFS) 

(c)  Pi CCSD  *  probability  of  a  correct  "Frame-Sync”  decision 

for  the  ick  time-slot  (i  «  N) 

c 

Cd)  P^(FD)  *  probability  of  a  false  dismissal  decision  for 
the  iti:  time-slot  (i  »  Nc) 

•  1  *  Pi  CCS) 

We  now  observe  that  the  conditional  probability  of  a 
correct  "Frame  Sync"  decision,  given  that  the  search  starts  on 
the  i  *  time-slot,  is 


? [correct  "Frame  Sync"  |  search  starts  on  the  ithslot] 
■  CPtCCD3 ]Nc'x  •  [Pi CCS) ] 

+  [PiCCDDjV1  •  [P.(FD)]  •  [PiCCD)]Nc*1  .  [P.(CS)] 


c*.i: 


[Pi(CD)]Nc*1  •  [PiCCS)3  •  [?i (FD) ] -  •  [Pi(CD} ]-  CNe 


j-0 


•1) 
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m 


In  (4.1),  the  j  *  term  is  the  probability  of  a  correct  "frame 


Sync"  on  the  (j  1)  ‘  test,  given  that  search  started  on  the 
slot. 

Therefore , 


th 


P  [correct  "Frame  Sync"] 


N, 


*c  i  -V*  i 

r  Z  f-  •  [P-CCDD 1  C  •  [?.(CS) ]  •  [P.(FD)J3 
j-0  i»l  1  1  1 


•  [PiCCD)! 


jCNc-l) 


3ut,  since 


Nc  V1 
Zc  X  c 


i-1 


M  -1  N  -2 
X  c  +  X  c 


---  *  X  +  1 


1  -  X 

nr 


We  obtain 

P [correct  "Frame  Sync"] 


«  Z 


l  eyes)]  •  [Pi(FD]j  *  [P±  CCD)  ] 


j(Nc-l) 


5T 


J-0  4,c 


T^W 


N 


•  U  -  [P^  (CD) ]  c> 

Additional  simplification  results  from 


OB 

Z  XJ  »  1  X  *  X2  - 

j-0 


*  r^T 


(4 


(4 


(4 


C* 
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with 


X  -  [?i (FD) ] 


[Pi  CCD)] 


N  -1 


(4.6 


This  leads  to 

? [correct  "Frame  Sync"] 

Nc 

[P-CCSD ]  .  U  -  [p.(cd)]  c> 

m  x  -  -  -  -  *  N1 

Nc[P.CFS)]  •  {1  -  [P.(FD)]  •  [P • (CD) ]  C*  } 


(4.7 


A  statistic  of  some  interest  is  E[j],  defined  and  e\'aluated 
below. 


E[j] 


I  j 

j-o  ir 


[Pi  CCS)]  •  [P.(FD)]j  • 

- [TJTTSr 


[Pi (CD)] 


j(Nc-l) 


Nr 

•  (1  -  [Pt(CD)]  C}  (4.8 

Using 


I  j  •  Xj  -  X  +  2X2  ♦  3X3  *  --- 

3-o 


(4.9 


yields 

N  N  - 1 

[Pi  CCS)]  .  a  -  [P  -  (CD)  ]  c}  •  [P .  (FD)  ]  •  [P.  (CD)]  c 

E[j]  ■  — - - - i - -  (4.1 

N  -1  * 

N  •  [Pi(FS)]  •  [1  -  [Pi (PD)]  •  [Pi  (CD)]  c  ) 


Perhaps  of  more  interest  than  E[j]  is  E[jicorTect  "Frame  Sync"] 
which  is  found  by  dividing  (4.10)  by  (4.7). 


43 


This  gives 


E[j  ]  correct  "Frame  Sync'*] 

N-1 

CP,  CTO)  1  •  [P s  CCD)  ]  c 

- 

1  -  [PiCTO)lCPiCCD)]  C 

In  using  C4*- *  10)  or  Cl-11).  it  should  be  recalled  that  (j  *1) 
represents  the  number  of  attempts  before  "Frame  Sync"  is  declared. 

Another  function  of  interest  is  the  probability  that  j 
is  equal  to  or  less  than  some  value  j  .  This  is  derived  below. 
PU£30I 

4  HVU 

jQ  CP.CCS3 1  •  [PiCFB)]3  •  [Pi CCD)1  c  •  {1  -  [P, CCD) ]  c> 

«  £ u  •L  •**  -  ■  ■  ...  r  ^ 

j»0  Nc  •  [PiCFS)] 


This  is  simplified  by  using 


1  +  X  + 


yielding 


i 


x 


o 


pu  1  JJ 


CPjCCS)]  -(1-[P.  (CD)]Nc?»a-([P.  CFD)]  •  CPjCCD)]^  V‘o  J 

V  [PiCFS)]  •{!-[?. CTO)]  •  [PiCCD]^*1} 
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Equation  (4.14)  could  be  used  in  solving  for  the  value  j  such 
that  P  [  j  <_  j  o )  is  a  given  value,  such  as  0.9,  C.95,  etc.  Thus, 
changes  in  j  ,  as  signal  and/or  jamming  conditions  change,  rep¬ 
resents  a  meaningful  criterion  for  comparison  (as  do  (4.7)  and 
(4.11)). 

We  next  derive  the  general  results,  corresponding  to  the 
above,  for  the  false  "Frame  Sync”  case.  Consider,  first,  a 
few  special  cases  for  the  probability  of  a  false  "Frame  Sync" 
decision  on  the  (i+l)t^1  attempt,  given  that  search  stated  on 
the  ith  time-slot.  For  example, 

P [False  "Frame  Sync”  on  1st  Attempt  |  Search  Starts  in  Slot) 

N  -i 

-  1  -  [Pi(CD)3  C  (4.1r 

P [False  "Frame  Sync"  on  2nd  Attempt  1  Search  Starts  on  ir^  Slot] 

N-i  N  -1 

»  [Pj  CCD)]  c  -[PiCFDn-a-tP.CCD)]  c  (4.16 

P [False  "Frame  Sync"  on  3rd  Attempt  [  Search  Starts  in  i1^  Slot] 

N-i  2  N  -1  N  -1 

-  [Pi  (CD)]  c  • [Pi(FD)]Z-[Pi(DC)]  c  - {1- [P- (DCD ]  c  }  (4.1" 

The  general  result  is 

? [False  "Frame  Sync  on  (i+l)^  attempt  J  Search  Starts  on  i'^Slot] 

VI 

«  1-  [?,.  (CD)]  c  for  1-0 

1 

N  *1  1  (2.-1)  (N  -1)  N.-l 

-  ^ (CD)]  C  • [?i (?D) ]  -[Pi  (CD)]  C  [1-[P.(CD)3  w  }, 
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Cfc.  i  *, 


Thus  we  have 


? [False  "Frame  Sync"  on  1st  Attempt] 
»c  M  -i 

-  s  rr*  •  a-CP^ccD)]'0  } 

1-1  *c 

Nc 

U-CPiCCD)]  c} 

*  1  :ic  -ct-LPtCCD3I> 


*  1  - 


(1-  C P t C CD ) ] *  c} 

M  -L^CrSjj 
c  i 


(4.13.. 


Equation  (4.19)  follows  from  the  result 


c  N-i  ^r”1  ‘V2 

VC  _xc  +xc 


2  X 
i«l 


---  +  X  +  1 


For  the  other  cases,  (1*1, 2, 3,---) ,  we  have 


P[False  "Frame  Sync"  on  (,L*lya  test] 


(4.20, 


X  -i 


-  E  *  (PiCCD)]  c  •(Pi(FD)]l-[Pi(CD)] 
i-1  lNc  1  1  1 


(4-1) (Nc-1) 


N.-l 

.{l-[?i(CD)l  c  > 


[P,  (PD)  ] 1 '  [P.-  (CD)]  Cl_1)  CNc*1}  -{1-[P.  (CD)]’  C  X{l-[?.  CCD)]'  C> 

*  x  *  _ 


Nc* {!-[?. (CD)} 


{[P4 (FD)]*[P, CCDD 


{ 1- [P . (CD) ]  u  }•{!-[?, CCD}] 


Nc- [P±  CCD) ] 


[Pj,  CFS)  ] 


Now  we  evaluate 


P [False  "Frame  Sync"] 


Z  P [False  "Frame  Sync"  on  (i+l]  test] 
1*0 


(Ne-l)-Ne[PiCCD)]*[Pi(CD)] 

Nc  -CPiCFS)] 

N  -1  N 

[Pi(FD)]-{l-Pi(CD]]  c  J-U-tPjCCD)]  C) 

- - - 

N  • [Pi CFS] ] • { 1 - [Pi  CPU) ] • [Pi  CCD] ]  C  } 


where  the  first  term  of  (4. 21]  is  [4.19]  in  different  form,  and 
the  second  term  uses  [4.20]  for  the  summation  over  the  (1^0] 
values . 

In  a  similar  manner  for  defining  E[j],  equation  (4.8], 
we  now  define  and  evaluate  E[i]  for  this  false  "Frame  Sync" 
case. 

E[i]  *  2  i*P[False  "Frame  Sync"  on  (i+l)^  test] 

1*0 


N  -1  N 

( 1* [?i (CD) ]  c  ) • (1- [?i (CD] ]  C 

N  -1 

N,[P,-  CCD]]  c  .[P4CFS]] 


I  !•{[?. (FD] ] * [P.  (CD]] 
1-1  1  1 


Again  using  C4.9),  with  the  ia°  rara  dra?Pi:1?  ctit »  we  obtain 


N  -1 


a-[P.(CD)]  C  }  •  (1-  rp.  (CD)  ]  c>-[Pi(FT})] 

sen  — L— - N  t r7* 

Nc-[PiCFS)]*{l-[PiCFD)]-(PiCCD)]  C  / 


To  find  E[ij False  "Frame  Sync"],  divide  (4.23)  by  (4.21). 

Mow  we  find  the  probability,  for  this  false  "Frame  Sync' 
case,  which  is  analogous  to  (4.14). 


o  th 

P[i<z  ]  *  Z  P [False  "Frame  Sync"  on  (i*l) “  test] 

“  0  2*0 


(4.:: 


M, 


(Nc-1)-Mc-[P; (CD)]>[PiCCD)]  C 
Nc  *  CPi(FS)  1 


N  -1 


M 


c, 


{1>[P. (CD)]  C  >«{!-[?! CCD) I  4, 


1CT 

Mc(Pi(CD)]  C  •CPi(FS)] 


Z°([PiCFD)] • CP^CCD) j 
2*1  1 


(4.24. 


,w 


To  simplify  (4.24),  we  can  use 


in  ?  ? 

z°  XZ  -  X  *  X“  ♦ 
2*1 


(4. 


Xll-X  °) 

1-x ' 


The  result  is 
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M 


Nc 

CN  -l)-Ne[P. (CD)]*[P.(CD)]  C 
PI1-£0J  *  - - L?iC?Sj] - 

N  - 1  N  N  -1  Z 

{1- [Pi  CCD)]  c  ).{1-[P.(CD)]  c> - [P±  CFD) ] • [P±  CCD) ]  C  )  °] 

* -  - 

N  -[Pi  CFS)].  a- [Pi  (FD)].  [Pi  (CD)]  C  } 


(4.2  t 


The  key  equations  above  represent  the  desired  set  of  gen* 
eral  results,  which  provide  the  basis  for  evaluation  of  the 
frame  synchronization  technique  with  specific  jamming  waveforms. 
In  the  next  section  we  consider  the  application  of  the  above 
equations  to  special  cases  of  interest. 


4 . 3  General  Results  Applied  to  Specific  Jamming  Cases 

The  general  results  presented  in  the  previous  section 
provide  the  framework  for  study  of  the  effects  of  jamming  on 
the  frame-synchronization  technique  described.  We  begin  by 
noting  that,  for  most  jamming  cases  of  interest  for  the  NCFSK 
system  being  considered,  the  channel  is  a  binary  symmetric 
channel.  That  is,  the  two  types  of  errors  occur  with  equal 
probability,  which  we  represent  by  the  symbol  p. 

For  the  first  Nc*l  time-slots,  the  receiver  will  produce 
"zero”  and  "one"  at  its  output  with  equal  probability.  This 


i 

i 


is  due,  of  course,  to  the  completely  random  nature  of  the  data 
in  these  slots  and  the  symmetric  nature  of  the  channel.  As 
a  consequence,  we  have 


r 


0-1  n  I  2  I  a«“J 
P£  CCD)  -  2  (  ?)C4)  (y)  S 

1  j-n.-b+l  1  “  " 


and 


1 

n_ 


b-l 

2 

j«n  -b+1 


P^FS)  -  l-P^CD) 


ns'b  n_  .  J  ,  ”s  '  "s  n  ,  J  . 

s  (4*3(7)  Cy)  +  ^  (4*3(7)  Cy) 

j  *0  J  “  j*bJ 


N  -  j  n 
s  J  s 


V* 


_Z_ 

a. 


n.-b 

a 

E 

j-0 


n 

C;S) 


Thus  P^(CD)  and  P^(FS),  for  substitution  into  the  key  equations 
of  the  previous  section,  are  independent  of  jamming,  since  they 
are  not  functions  of  the  bit  error  probability  p.  Instead, 
they  are  functions  only  of  ng,  the  number  of  symbols  tested 
for  the  sync  pattern,  and  b,  the  threshold  for  determining  if 
"Frame  Sync"  is  to  be  declared. 

On  the  other  hand,  P^(CS)  and  P^(FD),  which  apply  only 
for  time-slot  Nc  where  the  non- random  sync  pattern  is  located, 
are  functions  of  p,  and  therefore  depend  on  the  jamming 
waveform  assumed.  These  are  found  from 


Vb  a.  j  a.-j  n3  a.  -1  1 

PjCCS)  •  :  (  *H1-P)  Cp)  =  *  £  (/H1-?)  t?) 

j  *0  2  j  »b  J 


n,  -  j 

3 


30 


and 

P^FB)  *  l*?i  CCS) 


b-1 

*  z 

j  *n5-b+l 


Cj  ) Cl-p) 


: 


(4 


The  remaining  task  is  to  list,  for  various  jamming  wave* 
forms,  the  appropriate  expressions  for  p.  They  are  presented 
below. 

(1)  Random  Noise  Only,  Without  Jamming,  No  PH 

D  *  V 

*ao 


(2)  Random  Noise,  with  one  Frequency  Slot  Subjected 
to  Noise  Jamming,  No  FH 

?  *  p*i 

5 

(3)  Random  Noise,  With  Both  Frequency  Slots  Subjected 
to  Noise  Jamming,  No  FH 


iexp 


P 


P 


a2 


1*_ 


C4 


(4 


(4 
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(6)  Random  Noise,  With  Noise  Jamming  in  X  of  N 

d  3 


For  the  frequency-hopping  cases,  equations  (4.36)  and  (4.3"), 
the  following  substitution  is  used. 


'.-fc 

with  J  then  representing  the  total  jamming  power  over  the 
spread  bandwidth- 

la  addition  to  the  above  seven  jamming  cases,  there  are 
two  others  for  which  we  present  results.  Each  of  these 
represents  an  example  for  the  binary  nonsymmetric  channel 
which  requires  two  error-probabilities  to  characterize  it.  If 
we  let  PQ  be  probability-of-error  for  a  "zero"  transmitted  and 
p^  be  the  probability-of-error  for  a  "one"  transmitted,  it 
fallows  directly  that 
PCzero"  received! 


P(Or) 


l’Vi 


P  ["one"  received] 


P(lr)  * 


1+Po-?l 


We  must  find  the  new  expressions,  for  these  two  cases,  for  ?^(?S), 
Pi(CS),  etc.  The  results  below  apply  to  the  case  ns  even  and 
b  at  least  as  large  as  one-half  of  n  .  Similar  results  can 
be  found  for  ns  odd.  After  considerable  effort  we  obtain 
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As  with,  the  first  seven  cases,  only  the  binary-symbol  error- 
probability  remains  to  be  specified  in  order  to  define  a 
particular  case.  The  two  of  interest  are  specified  below. 
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with  I  (x)  being  the  modified  Bessel  funciton  of  the  first 
kind. 


1 

7 


l-Q  U|)1/2,  — 

.  ' 


1/2 


177 


(4 .  - 


e  1/2 

,S- 


1/2 


^  5  T77  «  > 


In  the  next  section  we  summarize  the  above  results  and  discuss 
the  steps  involved  in  using  them. 


4. 4  Discussion  of  Results 

The  above  set  of  comprehensive  results  provides  the 
analytical  tools  for  extensive  studies  of  the  effects  of 
jamming  on  frame  synchronization  for  the  multichannel  tactical 
system.  Both  spread  and  unspread  cases  are  included,  along 
with  a  variety  of  jamming  waveforms  and  strategies.  These  same 
tools  can  be  applied  toward  system  or  jamming  optimization, 
requiring  probably  many  sets  of  performance  curves  in  order 
to  locate  "optimum"  values,  etc. 

Since  there  are  many  cases  included  in  the  above,  it 
seems  highly  desirable  to  provide  the  following  table  which 
summarizes  the  nine  cases  and  gives  the  equation  numbers 
necessary  to  apply  the  formulas. 
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Summary  of  Results 


I.  General  equations  applicable  to  all  cases 

,4.7,  4.10,  4.11,  4.14,  4.21,  4.23,  4.26 

II.  Equations  for  cases  1-7 

4.27,  4. 23,  4.29,  4.30 

III.  Equations  for  cases  8,9 

4.40,  4.41,  4.42,  4.43 

IV.  Equations  for  individual  cases 


Case 

1 

- 

4.31 

Case 

2 

- 

4.32 

Case 

3 

- 

4.33 

Case 

4 

4.34 

Case 

S 

4.35 

Case 

6 

4.31, 

4.32, 

4  •  oo  > 

4. 36 

Case 

7 

4.31, 

4.34, 

4.35, 

4 . 37 

Case 

8 

4.44, 

4.45 

Case 

9 

4.46, 

4. 47 

Even  though  it  would  be  possible,  within  the  time  con-' 
straints  imposed  for  this  task,  to  present  a  few  results  in 
graphical  form,  it  is  believed  that  to  do  so  would  be  mis¬ 
leading  to  the  reader.  To  reach  meaningful  conclusions  con¬ 
cerning  the  effects  of  parameter  variations  and  changes  in 
EOt/ECCM  strategies  will  require  that  curves  be  found  for 
a  large  number  of  potential  situations.  Only  then  can  design 
and/or  operating  decisions  be  made  with  a  high  degree  of  con¬ 
fidence.  For  example,  results  obtained  by  the  use  of  the 
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analytical  tools  developed  here  may  lead  to  novel  techniques 
involving  adaptivity  of  the  number  of  time  channels,  frame- 
sync  pattern,  pattern  length,  threshold  setting,  jamming 
waveform,  jamming  strategy,  spreading  technique,  etc.  It  is 
highly  recommended  that  these  evaluations  be  carried  out  in 
an  extensive  manner  during  a  follow-up  effort,  with  final 
graphical  results  and  conclusions  presented  in  handbook 
fashion  for  ease  of  application. 

However,  in  order  to  illustrate  the  application  of  the 
above  tools,  the  following  sample  numerical  results  are  listed. 

Parameter  values  selected 

6 
6 

12 

16dB 

Noise  Jamming  in  both  frequency  slots,  no  fre¬ 
quency-hopping,  S/J  ■  ll.SdB 
Result:  P (correct  "Frame  Sync")  *  0.844 

Noise  Jamming,  frequency  hopping  with  Ns*1000, 
Partial-band  jamming,  P  (correct  "Frame  Sync")  ■ 

0. 844. 

Results :  Ks  S/J 

1  4.5  dB 

2  3  dB 

100  -9  dB 


b  • 

Nc  * 
S/N  - 

Example  1 

Example  2 
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5.x ample  5  Noise  Jamming  in  both  frequency  slots,  no  fre¬ 
quency-hopping,  S/J  *  ZdS 
Result:  ?  (correct  "Frame  Sync")  *  0.671 
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CHAPTER  3 


Low  Rate  Error  Correcting  Codes 
In  Spread  Spectrum  Systems 

Section  1.  Introduction 

Direct  sequence  spectrum  spreading  techniques  for  use  over 
low  signal- to-noise  ratio  channels  are  in  common  use  and  well 
understood.  We  briefly  summarise  the  salient  features. 

A  binary  digital  data  sequence  producing  digits  at  the  rate 
of  Rffl  digits  per  second  is  added  modulo- two  to  a  predetermined 
pseudorandom  (PN)  sequence  of  binary  digits  occurring  at  a 
rate  of  Rc>>Rm  digits  per  second.  The  resulting  digit  stream 
is  then  used  to  phase  modulate  an  RF  carrier  and  this  signal  is 
transmitted  over  the  given  channel.  At  the  receiver,  and  after 
synchronization,  the  in-phase  modulo- two  addition  of  a  replica 
of  the  PN  sequence  to  the  received  demodulated  waveform  will 
remove  the  PN  sequence,  resulting  in  a  data  sequence  approxi¬ 
mating  that  of  the  transmitter. 

The  quality  of  this  approximation,  as  expressed  by  the 
bit  error  rate  in  the  receiver  output,  is  normally  a  function 
of  the  so  called  processing  gain,  defined  by 


and  the  signal- to-noise  ratio  in  the  channel.  In  fact,  in  an 
additive  white  Gaussian  noise  channel  of  symbol  energy  to  noise 
power  spectral  density  Es/Nq,  a  binary  coherent  PSK  system  will 
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yield  a  bit  error  race  equal  to  that  of  an  unspread  channel 
with  a  signal  to  noise  ratio  given  by 


P.G. 


(5.1) 


Although  in  principle  any  sequence  of  binary  digits  will  serve 
as  a  spectrum  spreading  code,  the  most  commonly  used  configura¬ 
tion,  at  least  in  systems  without  stringent  security  requirements, 
is  the  so  called  maximal  length  or  pseudo- random  sequence,  whose 
basic  properties  we  describe  next. 

For  any  positive  integer  k>l,  the  maximal  length  sequence 

3c  k  - 1  k  - 1 

consists  or  2-1  binary  digits,  of  which  2  are  ones  and  2  -1 

are  zeros.  Every  possible  binary  k- tuple  can  be  found  once 

and  only  once  in  the  sequence.  The  autocorrelation  function 

v 

R(r)  has  the  value  2  -1  for  t»0  and  -1  for  all  other  integer 
values  of  r.  From  the  communications  point  of  view  these  prop¬ 
erties  assure  us  of  a  low  DC  component  in  the  binary  waveform, 
a  "random"  appearance  of  the  digits  and  optimal  synchronization 
performance  at  the  receiver. 

Beyond  these  certainly  desirable  features,  however,  the 
maximal  length  sequence  has  little  to  recommend  itself.  In 
particular,  since  the  pattern  of  zeros  and  ones  is  always  the 
same  and  hence  independent  of  the  data  sequence,  no  "coding 
gain"  beyond  the  bandwidth- -energy- -performance  tradeoff  implicit 
in  (5-1)  is  achieved. 

The  question  therefore  arises,  whether  it  is  possible  to 
influence  that  tradeoff  by  a  more  intelligent  choice  of  the 
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sequence  of  digits  transmitted  over  the  channel.  In  Section  2 
to  follow  we  examine  this  question  from  the  standpoint  of  error 
control  coding.  In  particular,  we  show  that  a  slight  change 
in  the  way  maximal  length  sequences  are  normally  used  allows 
us  to  exploit  their  characteristics  as  error  correcting  block 
codes.  We  then  derive  their  basic  properties,  stipulate  an 
error  correction  algorithm  and  obtain  performance  values  for 
the  codes  in  terms  of  bit  error  rates  and  data  deletion  rates 
over  the  hard  quantized  Gaussian  Channel. 

We  conclude  with  the  derivation  of  several  asymptotic 
bounds  on  performance  and  a  comparison  of  bit  error  rates  achieved 
under  an  optimum  correlation  decoding  scheme. 

In  Section  3  we  turn  our  attention  to  the  idea  of  using 
low  rate  convolutional  codes  as  a  means  of  spectrum  spreading. 

We  first  define  a  class  of  codes  with  good  structural  properties 
and  derive  their  free  distance  and  minimum  distance.  We  then 
show  that  £free  for  these  codes  is  optimum  in  the  limit  as  the 
constraint  length  becomes  large  and  compare  it  with  upper  bounds 
on  frree  for  finite  constraint  lengths. 

Next  we  give  generator  and  parity  check  matrices  and  obtain 
an  expression  for  the  generator  function  of  the  encoder  state 
transitions.  We  then  use  these  results  to  derive  upper  bounds 
on  the  residual  bit  error  rates  when  these  codes  are  used  on 
the  binary  symmetric  and  the  additive  white  Gaussian  channels 
and  decoded  via  the  maximum  likelihood  algorithm. 

Finally,  we  compare  the  results  with  more  conventional 
direct  sequence  spread  spectrum  schemes. 
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Sect* on  2.  Maximal  Length  Sequence  Codes 

2.1  Code  Definitions  and  Properties 

The  generation  of  the  pseudo-random  sequences  discussed  in 
the  introduction  is  usually  accomplished  by  means  of  k-stage 
shift  register,  in  which  the  contents  of  certain  stages  are 
added  modulo- two  and  applied  as  an  input  to  the  first  stage. 

When  such  a  register  is  loaded  with  a  particular  nontero  combi- 
nation  or  k  binary  digits  and  shifted  2*  -1  times,  the  sequence 
corresponding  to  the  starting  binary  k- tuple.  Loading  the  regis¬ 
ter  with  any  other  combination  of  the  digits  simply  results  in 

a  cyclically  shifted  output  sequence.  There  are,  of  course, 
v 

exactly  2*-l  cyclic  shifts  of  a  sequence,  corresponding  to  the 
v 

2* -1  nonzero  starting  k- tuples .  If  we  include  as  a  possiblity 

the  sequence  of  2  -1  zeros,  we  have  a  correspondence  between 

k  k 

2  distinct  sequences  of  binary  digits  and  all  2  binary  k- tuples. 

Without  great  difficulty  we  may  show  that  these  sequences  form 

a  cyclic  binary  linear  (n,k)  error  correcting  code,  where  n*2K-l. 

The  weight  spectrum  of  the  code  is  also  easily  determined  and 
-  k  k-1 

consists  or  2  -1  sequences  of  weight  V  ~  and  one  sequence  of 

weight  zero.  Finally,  the  distance  between  any  two  code  words 
k-1 

is  d»2  and  this  implies  that  the  code  can  correct  up  to 
,k-7 

e-2  -1  errors. 

Invoking  Plotkin's  bound  (I)  to  the  effect  that  the  mini¬ 
mum  distance  of  any  linear  binary  (a,k)  block  cede  is  bounded 
by 
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d  < 


n2 


k-1 


We  easily  show  that  these  codes  meet  the  Flotkin  bound  with 
equality  and  are  therefore  optimum. 

In  Figure  3 . 1  we  show  the  general  encoding  circuit  for  these 
codes  and  Table  1  lists  the  corresponding  connection  pattern, 
for  3<k£l0. 

Since  the  codes  are  linear  and  cyclic,  they  may  be  described 
by  a  parity  check  matrix  H  or  a  parity  check  polynomial  h(x). 

The  latter  is  also  defined  in  Table  5.1  and  the  former  is  shown 
in  Figure  5.2.  Note  the  convention  of  forming  a  codeword 
with  the  information  digits,  denoted  by  x^^  ,x^ , . .  .x^. ,  followed  by 
the  parity  check  digits,  denoted  by  »r2*  *  •  rn-k’  in  the  1151141 
way  for  cyclic  codes. 

As  an  example,  consider  the  case  k*3.  In  Figure  5.3  we 
show  the  encoding  circuit  and  Table  5.2  lists  the  resulting 
codewords.  In  each  codeword  the  first  three  digits  are  the 
information  digits  and  the  last  four  are  the  parity  checks. 

Consider  now  the  use  of  the  codes  described  above  as 

signalling  elements  in  a  spectrum  spreading  application.  Since 

with  the  exception  of  the  all  zero  sequence  the  codewords  are 

maximal  length  sequences,  the  balance  of  zeros  and  ones  and  the 

"randomness"  of  the  digits  are  roughly  preserved,  if  we  assess 

k 

tnese  properties  over  sequences  much  longer  than  2-1.  The  auto¬ 
correlation  function  is  no  longer  important  since  every  cyclic 
shift  of  a  sequence  is  also  a  legitimate  sequence.  This  fact 
of  course,  implies  the  need  for  some  other  means  of  synchroni¬ 


se 


OilCS 


Figure  5-3  Encoding  Circuit  for  k=3  Maxiaal  Lengti  Sequence  Code 


h-00  -  a.0^3.r^a2x2. .  .*akxk 


aoala2 '  *  ,ak 


k*3 

1+x+x3 

1101 

k»4 

i 

1+x+x 

11QQ1 

k-S 

-  ->  5 

l*x’+x3 

1Q10Q1 

k-6 

1>X2+XS 

110Q001 

k-7 

l+x^+x7 

10010001 

k-8 

l+x^+x^+x^+x3 

101110001 

k»9- 

..  1+x^-**9-  •  • 

•  L0G01QQQQ1 

k»10 

l+x3+x10 

1001000000 

Table  S.l 

Encoder  Connection  Patterns  and  Parity  Check 
Polynomials  for  Maximal  Length.  Sequence  Codes 


0  Q  Q  \  Q  Q  Q  Q 

0  0  1  f  1  1  Q  1 

0  1  0  |  Q  1  1  1 

a  i  i  [  i  a  i  q 

l  o  o  i  1  i  i  o 

1  0  1  1  0  Q  1  1 

i  i  o  s  i  o  o  i 

i  x  1 1  a  1  Q  o 

Table  5.2 

Codewords  of  C”,33  Maximal  Length  Sequence  Code 
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ration  of  the  sequences  such  as  a  preamble  or  a  cover  sequence. 
The  all  zero  sequence  could,  if  so  desired,  be  eliminated  from 
the  code  by  inserting  a  one  in  the  data  stream  after  k-1  con¬ 
secutive  zeros  prior  to  encoding  and  removing  it  again  at  the 
receiver  after  decoding.  This  would  require  some  elastic 
buffering,  result  in  a  decrease  in  data  rate  and  give  rise  to 
the  possiblity  of  error  propagation  in  the  decoded  output. 

Aside  from  these  difficulties,  however,  the  codes  achieve 
the  effect  of  spreading  the  energy  in  the  data  sequence,  the 
bandwidth  expansion  factor  being  n/k. 

2.2  A  Hard  Decision  Decoding  Algorithm 

We  proceed  to  study  the  performance  of  the  maximal  length 
sequence  codes  defined  in  Section  2.1  under  a  specific  decoding 
algorithm  which,  although  conceptually  quite  simple,  may  entail 
some  significant  increase  in  complexity  at  the  receiver,  com¬ 
pared  to  the  standard  direct  sequence  spread  spectrum  demodu¬ 
lators  . 

The  decoding  algorithm  is  defined  as  follows: 

Let  u  be  a  transmitted  codeword  and  v  be  the  corresponding 
received  sequence,  differing  from  u  in  i  places.  If  the  syn¬ 
drome  is  rero,  accept  v.  If  the  syndrome  is  not  zero,  change 

the  digits  of  v  successively  one  at  a  time,  two  at  a  time,  and 
k-  2 

so  on  up  to  e*2  -1  at  a  time,  until  the  resulting  sequence 

has  zero  syndrome,  i.e.,  is  a  codeword.  Accept  this  codeword 
as  the  correct  decision.  If  none  of  the  sequences  have  zero 
syndrome,  declare  an  incorrectabie  error  present.  It  is  clear 
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that  this  algorithm  works  correctly  if  the  number  of  errors 

k-  2 

made  by  the  channel  noise  does  not  exceed  e»2*  "-1.  It  is  also 

obvious  that  an  uncorrectable  and  undetectable  error  pattern 

k- 1 

will  result  in  a  decoded  output  error  of  multiplicity  2  ,  since 

the  decoder  in  this  situation  produces  an  invalid  nonrero  code- 

k- 1 

word  which  must  differ  from  the  correct  one  in  2‘  places. 

Finally,  if  the  error  pattern  is  detectable  but  not  correct¬ 
able,  the  decoder  will  pass  the  received  sequence  unchanged  and 
there- fore  the  multiplicity  of  the  error  pattern  at  the  decoder 
output  will  be  that  of  the  decoder  input. 

Thus,  if  i  is  the  multiplicity  of  the  decoder  input  error 
pattern,  the  decoder  output  error  pattern  will  be  of  multiplicity 
0,d*2  or  i.  It  remains  to  count  the  number  of  decoder  input 
error  patterns  giving  rise  to  a  decoder  output  error  pattern 
multiplicity  in  each  of  the  three  categories  listed  above. 

Let  these  quantities  be  denoted  by  and  >T  ,  with  the 

convention  that  .^-Q  for  i»0.  and  ^  are »  course, 

also  the  number  of  correctable,  uncorrectable  and  undetectable, 
and  uncorrectable  but  detectable  decoder  input  error  patterns, 
respectively. 

For  i<e  we  have  since  all  input  error 

patterns  in  this  case  will  be  corrected. 

For  i*e+l,  the  decoder  is  unable  to  produce  a  codeword 
since  it  is  restricted  to  making  at  most  e  changes.  Thus  all 
input  error  patterns  of  multiplicity  e+1  are  detectable  and 
we  have  No*N^*0, 

Next,  we  look  at  the  case  e*2<i<e+2«d.  Since  i  is  beyond 
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the  correction  capability  of  the  code,  NTo-0.  To  compute 
we  reason  as  follows: 

The  decoder  must  produce  a  nonzero  codeword  by  making  e  or 
fewer  changes.  Equivalently,  a  nonzero  codeword  must  result  in 
a  given  error  pattern  of  weight  i  by  making  e  or  fewer  changes. 
Now  a  codeword  has  exactly  d  ones  and  d-1  zeros.  If  we  change  j 
zeros  into  ones,  we  must  also  change  d+j-i  ones  into  zeros  in 
order  to  obtain  a  vector  of  weight  i.  There  are  ways  of 

changing  i  zeros  into  ones  and  for  each  of  these,  there  are 
^•d+j-i^  wa^s  °*  changing  d+j-i  ones  into  zeros.  Finally,  since 
the  total  number  of  changes,  given  by  j+d+j-i,  is  limited  to  e, 
j  may  vary  from  0  to[ — ^ — ],  where  the  symbol  [x]  denotes  the 
largest  integer  not  larger  than  x.  In  summary,  a  given  nonzero 
codeword  is  associated  with 


Z 

3-0 


(d-iw  d 
^  j  J  '•d+j-i 


) 


different  error  patterns  of  weight  i  and  since  there  are  n  nonzero 
codewords  we  have 


N 


Z 

3-0 


All  remaing  error  patterns  of  weight  i  are,  of  course,  detect¬ 
able  by  virtue  of  the  fact  that  they  cannot  be  turned  into  code¬ 
words  or  with  e  or  fewer  changes. 

Thus, 
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N.  =* 

L 


In  similar  fashion, we  derive  the  values  of  N  ,.W. ,  and  N.  for 
2e+2<Ze+2  and  3e+3<i<de+3 .  These  results  are  summarised  in 
Table  S . 3 . 

2.3  Performance  Over  the  Binary  Symmetric  Channel 

Consider  now  the  performance  of  these  codes  over  the  binary 
symmetric  channel  with  error  probability  p  and  under  the  decoding 
algorithm  given  above.  We-  evaluate  wo  parameters  of  code 
performance . 

1.  The  decoder  output  bit  error  rate,  defined  as 

3ER  »  ~£  {number  of  errors  in  decoder  output  sequence1- 

2.  The  rate  of  undetected  decoder  output  bit  errors, 
defined  as 

UBER  »  jE  {number  of  undetected  errors  in  decoder  out¬ 
put  sequence) 

Here  E{x)  is  the  expected  value  of  the  random  variable  x. 
For  the  situation  where  decoder  output  sequences  known 
to  be  in  error  are  deleted  we  also  define 

3.  The  rate  of  deleted  data  bits,  given  by 

DBDR  ■  -  Pr  {decoder  output  sequences  is  deleted) 

With  the  aid  of  Table  5.3  we  then  obtain 

3ER  -  §  Ce*l)  Ce!JL)?**Iqa'C*+1) 


ri-e-2, 

l — 5 — 1 


r  dn  z  d..)  pV*1 

i*e*2  j»Q  2  1 
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Decoding  Algorithm 


3e+2  1  2  J 

*  Z  dn  2  Ci^JCpPV 
i-2e+3  j-0  1  a  J  j 


d-1  .  ,d. i_n-i 


2e+2 
Z  i 
i*e+2 


V-n  }z.o  cli1,c4-i) 


I 

?  3 


3e+2 
S  i 
i»2e-*,5 


j-o 


i  n- 
p  q 


-le+3 

%  .  ,n,_i  n-i, 

+  £•  1UP  q  / 

i*3s+3  1 


Using  the 'fact  that  the  mean  of  the  binomial  distribution  is 
given  by  np,  this  reduces,  after  some  algebra,  to 


BER  » 


t  kJjpV'1 

i-1  x 


2e*2 

*  Z  Cd-i) 
i»e+2 


[±^! 
z 

3-0 


C  j  Hd.i,j)P  q 


C3.2) 


3e*2 

l  Cd-U 

i»2e-*-3 


[3e^i] 

Z 

j-0 


(  d-1  \  fd-.  i  n*i 
Ci-d*j)Cj}P  * 


The  undetected  bit  error  rate  is  also  easily  found  and  is  given 

by 


7*  - 


USER 


2e*2  1 2 J  .  ,  , 

-5  %  -  rC-l>r  « 

a  .  -  ,  C  j  Ud+j-i)p  ^ 

i*e+2  3*0  *'  J 


l  n-i 


+  d  Z 


ri-e-2, 
3e  +  ?  L — l J 


(5.3) 


t  t  c  *v  —  i n* i 

.  t  -  -n  Ci-d+j3  ^0?  * 
i*2e+3  j*0  J  * 


For  our  last  performance  parameter,  namely  the  bit  deletion 
rate,  we  obtain,  after  some  simplification, 


DBDR  *  |  (1  -  I  C?)p1qn"1} 
n  i*0  1 


rl-e- 2  ■* 

2e+2  L  2  J  ,  1  A  .  . 

■k  .  i  < 

x-e+2  3*0  J  u  ; 


(3.4) 


r5e+2-i-, 

3e+2  *■  2  ^  «  «  j 

•V  t  7  r  n  /d>  i  n*i 

k  1  L  Ci.d^i j H )P  q 

i* 2e+3  3*0  1  Q  ^  2 


If  we  assume  binary  coherent  PSK  communications  over  the  channel, 
the  relationship  between  the  channel  error  probability  p  and 
the  ratio  of  symbol  energy  E$  to  noise  power  spectral  density 
NT0  in  the  channel  is  given  by  the  well  known  formula 


1  -rff£s 


o  *  ▼  •  e 
•  » 


.  “  *»  .•‘wn’clw  '■ir.k  kh  «•*.!*<«.  cm*  w>iWU<Ha»>i 


Table  5.4  lists  this  relationship  for  several  values  of  ?.  In 
terns  of  the  energy  per  data  digit  we  also  have 

Eb  ,  n  Es 

With  the  aid  of  these  equations  we  may  now  calculate  3ER,  USER, 
and  D3DR  as  a  function  of  =^/^0-  ?h*  results  are  shown  in 

Tables  3.5  through  5.10,  for  n*7,  15,  31,  63,  127  and  225. 

Figure  S-4.‘is  a  plot  of  logBER  vs  (Eb/NQ)  in  decibels  for 
the  maximal  length  shift  register  codes  using  hard  decisions  and 
the  bounded  distance  decoding  algorithm.  The  characteristic  steep 
threshold  behavior  is  evident.  Thus  for  E^No  above  approximately 
10  d3  the  bit  error  rate  starts  to  fall  off  very  rapidly.  This 
threshold  corresponds  to  the  limiting  performance  capability  of 
these  codes  in  a  binary  symmetric  channel  and  is  consistent  with 
the  asymptotic  performance  bounds  derived  in  the  next  section. 

As  a  point  of  contrast,  it  should  be  observed  that  the  best 
codes  operating  with  soft  decisions  over  the  additive,  white 
Gaussian  noise  (ANGN)  channel  would  have  threshold  at  -1.6  d3 
Cthe  Shannon  limit  of  the  AWGN  channel) .  Indeed  these  very 
same  codes  when  soft  decision  decoded  with  a  bank  of  correlators 
in  fact,  approach  the  performance  specified  by  the  Shannon 
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3.23 
1.87 
0.83 
-0.04 
-0.82 
-1.62 
-2.34 
-3.03 
-3.76 
-4.42 
-S.27 
-3.99 
-6.84 
-7.74 
-8.69 
-9. 37 
-ID. 55 
-11.88 
-13.18 
-14.71 
-16.99 
-18.93 
-23.01 
-27.43 


Table  5.S  Bit  lirror  Hates  uiul  Hit  Deletion  Hates  for 
Maximal  Length  Sequence  Codes  Over  HSC 
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Table  5.6  Bit  lirror  Rates  ami  Bit  Delct ion . . . (etc) 
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Table  5.8  Hit  Hrror  Rates  and  Hit  Deletion  Rates  l?or 
Maximal  l.engtli  Sequence  Codes  Over  the  HSC 
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Table  5.9  Hit  lirror  Kates  ami  Bit  Deletion  Kates  for 
Maximal  Length  Sequence  Codes  Over  the  BSC 
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Table  5.10  Bit  Hrror  Rates  and  Bit  Deletion  Rates  for  Maximal 
Length  Sequence  Codes  Over  the  BSC 
n  ■  255 


log  I1HR 


limit.  This  is  discussed  in  section  2.4  and  the  corresponding 
curves  are  shown  in  Figure  5-9. 

On  the  basis  o£  these  tables  and  curves  we  make  the  fol¬ 
lowing  observations: 

1)  If  we  arbitrarily  define  an  acceptable  bit  error  rate 
as  BER<10*'5  and  an  acceptable  data  bit  deletion  rate  as  1%  of 
the  data  transmitted,  i.e.,  DBDR<0.01  k/n,  then  the  minimum  re¬ 
quired  value  of  E^/N  increases  with  code  length  from  7db  at 
n*7  to  11  db  at  n*223. 

2)  For  the  most  practical  cases  of  n»127  and  n*2S5  the  un¬ 
detected  bit  error  rate  U3ER  is  essentially  zero,  for  all  values 
of  E^/Nq.  Thus,  even  though  the  error  correcting  ability  of  the 
decoder  becomes  marginal  at  E^/N0  less  than  7db,  virtually  none 
of  the  residual  decoder  output  errors  go  undetected.  This  fea¬ 
ture  appears  to  be  the  one  of  major  importance  ir  the  applica¬ 
tion  of  the  maximal  length  sequence  codes  to  spread  spectrum 
systems. 
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2.4  Asymptotic  Performance  Bounds 


For  large  values  of  n  the  Plotkin  bound  on  minimum  distance 
for  a  binary  linear  block  code  takes  the  asymptotic  form  (1) 


k 

n 


< 


1-4 


d 

In 


Letting  e*2d-l  can  be  the  code's  error  correcting  capability, 
this  can  be  rewritten  as 

.  s  <  — 

The  implication  of  this  formula  is  that  the  error  correcting 
capability  of  a  bounded  distance  hard  decision  decoding  algorithm 
for  sufficiently  long  linear  (n,k)  block  codes  is  limited  to 
one  fourth  the  number  of  parity  digits  in  a  codeword. 

If  this  code  is  used  on  a  binary  symmetric  channel  with 
bit  error  rate 


P  > 


lfn-k. 


f 


the  code  is  unable  to  correct  even  the  average  number  of  errors 
np  in  a  codeword.  Thus,  in  an  absolute  sense  of  the  upper 
limit  on  the  correctable  channel  bit  error  rate  is  given  by 


P 


max 


< 


n-k 

in 


For  the  maximal  length  sequence  codes  n-k-n  as  n  becomes 
large  and  Pna:c-1/**  T^is  is  equivalent  to  Es/N'o»<:.25  db ,  which 
represents  the  minimum  channel  signal-noise  ratio  for  which  the 


maximal  length  sequence  codes  are  useful  under  a  hard  decision 
decoding  rule,  using  BER  as  the  criterion  of  performance. 

Let  us  now  derive  an  asymptotic  expression  for  BER  in  the 
limit  as  n  becomes  large,  assuming  p<l/4.  Without  difficulty 
we  may  show  that  the  double  sum  in  equation  (5.2)  is  then  neg¬ 
ligible,  so  that 


BER  *  p-i  I  i C^)piq 
n  i-1  1 


n-i 


Since 


n 


2  UfjpV1  *  np 

i-1  1 


this  can  be  rewritten  in  the  form 


BER  £  |  Z  i (f)piqn’i 
n  i-e+1  1 


If  we  then  let  i-1  d  r,  then 


n-1 


BER 


n- 1>,  r  n- r- 1 


P  I  rp*)p‘q 

r»e 


The  sum  is  now  the  "tail"  of  a  binomial  distribution  which 
we  can  overbound.  From  Peterson  and  Weldon  [1]  we  have 


":l  C":1)  ?V!-'  < 

A (n-1)  *  ■ 
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where  u  *  1-X 

provided  \>p.  [This  is  the  sane  as  saying  that  the  sun  starts  t 
the  right  of  the  mean  of  distribution.] 

For  the  codes  in  question, 


e  *  hr  i-  xCa’l3 

or 

n-3 

:  .  T~  1 
A  I=T  *  T 

as  n-*<®.  Hence  far  all  p<0.25  the  above  inequality  holds  asympto 
ically.  Then  substituting  we  have 


-  X(a-l) 
3ER  <  pC|) 


uCa-i) 


for  large  n,  X-i  and  upl¬ 
and 


BER  <_  p(9. 48  pq^) 


as  an  asymptotic  upper  bound  on  BER  for  all  p<l/4. 

Next  we  establish  approximations  to  BER,  LIBER,  and  DBDR  valid 
for  values  of  p<0.02.  In  this  case  the  first  terms  of  equations 
(5. 2),  (5. 3)  and  (S.4)  suffice  and  we  have 


BER 


(n-l)pe*iqa-e 


1 
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UBER  %  dCf)?e  +  2qn‘e'2 
DBDR  %  “Y(n*1)pe+1qI1*e'1 

As  an  example  o£  the  quality  of  these  approximations  we  evaluate 
equations  (3.6)  for  the  p»0.02  and  n»7.  Then 

BER  £  6p2q3  •  2.17  x  10"3 
UBER  £  16p3q4  -  1.18  x  10-4 
DBDR  %  9p2qS  »  3.23  x  10-3 

Note  the  excellent  agreement  with  the  exact  values  of  Table  5.3. 
2.3  Correlation  Decoding 

With  the  fundamental  limitations  of  maximal  length  sequence 
codes  with  hard  decision  decoding  well  established,  we  turn  our 
attention  to  their  performance  under  a  correlation  decoding 
regime. 

We  assume  that  each  of  the  n-1  maximal  length  sequence 
codewords  ui »u2 » •  •  •  un- 1  cont?os«^  of  n  elemental  antipodal 
signals  sQ(t),  s^(t),  where  0<t<T/n  and  T  is  the  duration  of 
the  codeword.  The  set  of  codewords  then  forms  a  so  called 
transorthogonal  or  regular  simplex  code  with  correlation 
coefficient  -l/(n-2)  between  unlike  codewords  (2). 

Let  the  codewords  be  transmitted  with  equal  probability 
and  received  in  white  Gaussian  noise  of  power  spectral  density 
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Ng  watts  per  hertz.  With,  coherent  reception  the  optimum  decoder 
is  then  easily  shown  to  take  the  font  indicated  in  Figure  5.3. 
Here  v(t]  is  the  received  signal  corresponding  to  the  trans¬ 
mission  of  a  codeword  and  the  maximisation  of  the  outputs  from 
the  n-1  integrators  takes  place  at  the  end  of  the  time  period  T. 

We  may  also  show  that  under  this  decoding  scheme  the  code 
bit  error  rate  is  minimum  and  given  by 

'"J 


3ER 


n-1 

Tn~ 
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1 

,/Z7 


? 

-x-/2 


k£.  (n-1) 

erfc(X” 


where  is  the  energy  per  information  bit  and  erfc(x)  is  defined 
by 

x  2 

erfc(xj  *  — i—  e*^  ^2dy 

•  03 

Figure  5.9  shows  BER  as  obtained  from  a  numerical  evaluation 
of  the  above  formula  ( 2 ),  as  a  function  of  H^/Nq,  for  a  number 
of  values  of  k.  If  we  again  take  10* J  as  the  maximum  accepta¬ 
ble  value  of  BER,  then  the  minimum  required  H^/h'o  decreases 
with  increasing  k  from  4  db  at  k*3  to  2  db  at  k*8.  Comparing 
these  values  with  Section  2.4  it  is  clear  that  the  correlation 
decoding  significantly  improves  upon  the  performance  of  these 
codes  over  the  Gaussian  channel.  The  complexity  of  implemen¬ 
tation,  however,  is  likely  to  be  much  greater  than  for  hard 
decision  schemes. 


Section  3.  Low  Rate  Convolutional  Codes 
3.1  Code  Definition  and  Properties 

A  binary  convolutional  code  of  rate  R  ■  1/n  and  constraint 
length  k  is  a  mapping  of  a  semi- infinite  sequence  of  data 
digits  into  a  semi- inf inite  sequence  of  channel  digits  such 
that  each  data  digit  corresponds  to  n  channel  digits  which  are 
functions  of  the  last  k  data  digits.  The  implementation  of 
such  a  coding  scheme  is  accomplished  in  terms  of  a  socalled  for¬ 
ward  shift  register  consisting  of  k  stages,  from  which  n  sub¬ 
sets  of  the  k  stages  are  connected  to  modulo  -  2  adders.  The 
outputs  from  these  adders  are  then  sampled  once  per  data  digit 
and  transmitted  over  the  channel. 

Letting  1  and  0  denote  connection  or  no  connection,  respec¬ 
tively,  of  a  given  stage  to  a  given  modulo  -  2  adder,  the  code 
is  completely  defined  by  a  set  of  n  binary  k-tuples,  the  so- 
called  code  generators  G.,G,...G  . 

As  an  example,  consider  the  case  k*3.  One  possible  form 
of  the  encoder  is  shown  in  Figure  5.  io.  Clearly,  the  two 
generators  are  ■ 

G1  •  111  ;  G2  -  101 

and  the  code  has  rate  R-l/2. 

In  general,  if  we  insist  on  each  generator  being  distinct, 

v 

then  n  is  limited  to  n<2  -1.  In  the  remainder  of  this  section 

]( 

we  consider  the  case  n«2  -1,  so  that  the  generators  are  all 
binary  nonoero  k-tuples.  The  rate  of  these  codes  is  then 
R«1/(2X-1)  and  the  channel  bandwidth  is  expanded  by  a  factor 


Figure  5.10  3inary  Convolutional  Code  of 

Rate  of  1/2  and  Constraint  Length  3 


of  : 


In  Table  3.11  we  have  listed  the  set  of  generators  as  the 

rows  of  an  array.  Without  great  difficulty  we  may  show  that 

except  for  the  absence  of  all  tero  k-tuple,  the  rows  of  the 

array  form  a  vector  space  of  dimension  n  over  the  binary  field. 

Each  column  and  each  modulo- 2  sum  of  two  or  more  columns  con- 
k- 1  k- 1 

tains  exactly  2  ones  and  2  -1  zeros.  In  fact,  the  columns 

can  be  thought  of  as  the  basis  vectors  for  a  vector  space  of 
dimension  k  which  is  exactly  the  set  of  codewords  of  the  (n,k) 
maximal  length  sequence  code  discussed  in  Section  2. 
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Table  5.11  Generators  of  n*  2^- 


1  Convolutional  Code 
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Under  the  usual  definition  for  convolutional  codes  [1]  , 
the  general  form  of  the  k  x  !<C2k-l)  generator  matrix  G  is  given 
by 


r 


c  - 


T 

c2 


ck 


T 

ck 


T 

where  is  the  transpose  of  the  ith  column  of  Table  5.11  and 
the  empty  places  in  G  contain  all  zeros . 


Since  the  first  digit  of  c^  is  always  one  and  the  first 
digits  of  c2, c3>''*c’£  ar®  always  zero,  the  code  is  clearly 
systematic.  The  parity  check  matrix  is  therefore  given  by 


where  cV  is  equal  to  with  the  first  digit  deleted  and  I  and  0 

k  k 

are  unit  matrices  and  zero  matrices  of  sice  2-2x2  -2. 


As  an  example,  consider  the  case  k*3.  Then  the  generators 
in  array  form  are 


C1  c2  c3 

1  0  0 

0  10 
0  0  1 

1  1  0 

10  1 
0  11 
111 


G  « 


10011010101011001  0*  111 
000000010011010101011 
000000000000001001101 


and  its  parity  check  matrix  takes  the  form 


From  the  standpoint  of  performance,  the  most  important  para 


meter  of  a  convolutional  code  is  its  free  distance  d-_  ,  de- 

fined  as  the  number  of  ones  in  that  semi-infinite  channel  sequence 

with  the  smallest  nonzero  number  of  ones.  We  now  derive  this 

quantity  for  the  codes  defined  above. 

Suppose  the  data  sequence  to  be  encoded  consists  of  a 

leading  one,  followed  by  an  arbitrary  semi- infinite  sequence  of 

ones  and  zeros.  Then  as  the  leading  one  enters  the  leftmost 

stage  of  the  encoding  register  and  traverses  to  the  last  stage, 

the  successive  channel  n-tuples  correspond  to  sums  of  columns 

of  the  array  in  Table  5.11  and  are  hence  codewords  of  a  maximal 

length  sequence  code.  The  total  number  of  ones  in  each  output 

k-1 

n-tuple  is  therefore  always  2  and  since  the  leading  one  tra¬ 
verses  k  stages,  the  output  sequence  corresponding  to  a  semi- 

k  “1 

infinite  input  sequence  must  have  at  least  k2  ones.  Conse¬ 
quently, 


d£re^-k  jk'1 


On  the  other  hand,  for  the  data  sequence 


100 . 0 . 

k-1 

the  output  sequence  has  exactly  k  2  ones  and  we  have 


dfree  "  k  2 


,k-l 


C3 


Next,  we  compare  d£  in  C5.")  with  a  well  known  upper 
bound  (3)  which  for  codes  of  rate  R»l/C2  -1)  we  may  state  as 


We  note  that  for  large  values  cf  k  this  bound  approaches 

v.  1 

k2‘  ,  so  that  the  codes  defined  above  have  maximum  free  dis¬ 

tance  for  large  constraint  lengths. 

In  Table  3 12  we  compare  equations  C3-7)  and  C3.3]  for 

v 

k<10.  Also  included  are  the  free  distances  of  good  R-1/C2*  -1) 
codes  constructed  from  optimum  R*l/2,  R»l/3  and  R*l/4  codes 
by  generator  replication  [3] . 

Another  important  parameter  of  performance  is  the  code’s 
minimum  distance,  defined  in  the  usual  way  as  the  minimum  dis¬ 
tance  among  k.  successive  blocks  of  encoder  output  digits  ccrres 
ponding  to  k  successive  encoder  input  digits  corresponding  to  k 
successive  encoder  input  digits  of  which  the  first  is  a  one. 
Since  each  block  of  2* -1  decoder  outputs  is  a  codeword  in  a 

maximal  length  sequence  code,  the  distance  between  such  blocks 
k- 1 

is  exactly  2‘  .  Thus  we  have  immediately 

d  -  -  k  Z*^'1  -  dr 

mm  free 

and  the  code  can  therefore  correct  all  patterns  of  up  to 

V-  7 

e  -  k  2*  -1 


v 

errors  among  any  successive  k C 2  - 1 )  channel  digits. 

Although  the  calculation  of  a  convolutional  code's  per¬ 
formance  over  a  given  channel  and  under  a  specific  decoding 
algorithm  is  generally  a  difficult  problem,  good  upper  bounds 
on  the  basic  bit  error  rates  can  readily  be  obtained  if  the  so- 
called  generator  function  of  the  encoder  state  transitions  is 
known.  We  now  turn  our  attention  to  a  derivation  of  this  ouan: 


100 


K 

max  dfree 

1,  1 
h.2*"  * 

constructable 

d£ree 

3 

18 

12 

18 

4 

50 

32 

50 

»• 

3 

124 

80 

124 

6 

288 

192 

283 

7 

6  55 

448 

63S 

8 

1457 

1024 

1402 

9 

3212 

2304 

3066 

10 

7014 

5120 

6905 

Table  5.12  Free  Distance  Comparisons  of 

v 

R*l/(2  -1)  Convolutional  Codes 
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Since  the  calculations  involved  are  somewhat  tedious,  we  restri 


ourselves  to  the  case  I<*3  for  purposes  of  illustration  and 
merely-  state  the  result  for  general  values  of  fc. 

Let  the  states  of  the  encoding  register  of  Figure  3.6 
correspond  to  the  four  combinations  00,01,10,11  of  the  contents 
of  the  last  two  register  stages  and  denote  these  states  by  x,  , 
x2>  Xj,  x^,  respectively.  Then  it  is  clear  that  under  a  tero 
input  the  states  undergo  the  transitions 


where  x£  is  the  tero  state  reached  from  a  nontero  state. 
Likewise,  under  a  one  input  we  have 


Using  standard  signal  flow  graph  techniques  we  now  reason 
as  follows: 

For  every  state  transition  the  output  from  the  encoder  is 
always  a  codeword  from  the  maximal  length  sequence  code  and  so 
has  weight  4,  except  for  the  transition  x^-x^.  The  input  to  th 
encoder  corresponding  to  a  given  state  transition  is,  of  course 
either  tero  or  cne.  Let  the  exponents  of  the  variables  D  and 
N  correspond  to  the  weights  of  the  output  and  input,  respective 
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5.2  Performance  Bounds 

With  the  aid  of  equations  (5.7)  and  (5.8)  we  may  easily 
obtain  upper  bounds  on  the  residual  bit  error  rate  when  the  above 
codes  are  used  for  signalling  over  the  additive  white  Gaussian 
noise  channel  and  decoded  via  a  maximum  liklihood  algorithm. 

The  pertinent  formulas  are  (4) :  r 

a.  For  hard  decision  decoding: 


N-l,  D 


where  p  is  the  channel  bit  error  rate  and  q  »  1-p. 
b.  For  soft  analog  decision  decoding: 


P,  <  erfc'  J  ET 


dr  RvT* 
e  free  N0 
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<  triz'Jl 


We  may  compare  this  result  with,  the  case  of  no  coding  by  noting 
that  for  the  latter 


pb  •  er£c  v2K7 


Clearly  then,  for  large  the  so  called  coding  gain  of  the  codes 
discussed  here  satifies  the  bound 


k- 1 

Coding  gain  >  dfre.R  -  Sfc—  -  | 


CS. 


which  is  a  monotonically  increasing  function  of  the  constraint 
length. 

Eb 

For  finite  the  denominator  of  (5.11),  which  for  reasonably 
o 

large  k  assumes  the  form 


a 


k-l 

I 

i-1 


i  Eb 


2 


must  of  course,  be  taken  into  account.  For  k*"  and  5^/No»3db, 
for  example,  its  value  is  approximately  0.2.  The  numerator 
in  this  case  equals  erfc'  /l4 . 1  %  10~4,  so  that 


<  3x10 
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Mote  that  this  bound  compares  favorably  with  the  performance 

of  the  standard  rate  1/2,  constraint  length  7  convolutional 

code  decoded  via  the  soft  decision  Viterbi  algorithm  [4] . 

As  our  last  special  case,  let  E^/N  -HD.  Then  the  denomi- 

•> 

nator  of  (3-11)  approaches  the  value  C2-k)".  For  k>3  this  is 
no  smaller  than  one  and  thus  the  codes  exhibit  a  coding  gain 
which  is  bounded  by  C5.12) 

3.3  Concatenated  Convolutional  Codes 

An  alternate  means  for  obtaining  a  low  rate  convolutional 

code  is  by  means  of  cascading  or  concatenating  a  single  moderate 

rate  code.  For  example,  if  the  output  of  a  rate  1/2  encoder  is 

fed  into  a  rate  1/2  encoder  the  result  is  a  rate  1/4  code.  Then 

v 

concatenating  k  times  can  yield  a  rate  1/2  code.  Ir  the  same 
code  is  used  in  each  stage  of  the  concatenation,  then  a  single 
Viterbi  Algorithm  decoder  can  be  used  for  decoding.  The  first 
inner  decoder  will  give  a  coding  gain  in  that  the  BER  at  its 
output  will  be  improved.  If  the  output  of  the  first  inner 
decoder  is  then  applied  to  the  input  of  the  second  inner  decoder, 
the  raw  bit  error  rate  should  again  be  improved  provided,  of 
course,  that  the  second  decoder  bit  input  errors  behave  as  if 
they  were  derived  from  a  binary  symmetric,  memoryless  channel. 
Unfortunately,  the  first  decoder  output  errors  are  known  to  be 
bursty.  However,  it  is  nevertheless  possible  to  achieve  the 
coding  gain  per  stage  by  interleaving  and  deinterleaving  .  The 
coding  gain  advantage  of  these  codes  have  been  explored  in 
this  study  but  it  is  felt  that  because  of  the  excessive  delays 


in  the  interleaving-deinterleaving  process  and  because  the 
"complete"  convolutional  codes  described  in  the  previous  section 
appeared  so  promising,  that  we  should  postpone  further  inquiry 

into  concatenated  convolutional  codes. 

\ 

\ 

3.4  Conclusions 

We  have  investigated  the  basic  properties  of  a  class  of  low 

\ 

rate  convolutional  codes  for  the  purpose  of  assessing  their 

usefulness  in  direct  sequence  spread \spectrum  applications. 

'  \ 

In  summary,  we  note  the  following  points: 

a.  Compared  to  the  generation  of  PN' sequences ,  the 
k 

2-1  tap  configurations  of  the  convolutional  encoders  imply 
a  somewhat  larger  complexity,  although  from  the  standpoint  of 
LSI  technology  this  seems  of  little  real  consequence. 

b.  Since  maximum  liklihood  decoding  of  convolutional  codes 
is  in  a  sense  self-synchronizing--although  at  the  loss  of  some 
data--the  need  for  PN  sequence  synchronization  at  the  receiver 
is  obviated.  This  appears  to  be  a  major  advantage  of  the  low 
rate  coding  approach.  Of  course,  bit  synchronization  is 
required  in  both  cases,  representing  the  fundamental  limit  to 
spread  spectrum  operation  at  low  signal-to-noise  ratios. 

c.  Due  to  their  extremely  long  branch  sequences  the  codes 
discussed  above  appear  to  be  of  particular  value  on  bursty 
channels,  elimating  or  at  least  reducing  the  need  for  interleaving. 
The  savings  in  storage  requirements,  decoding  delay  and  inter¬ 
leaver  array  synchronization- -to  name  just  a  fev--are  likely 

to  be  appreciable. 


d.  Unlike  standard  direct  sequence  methods,  low  rate  coding 
exhibits  a  coding  gain  beyond  the  normal  processing  gain  achieved 
through  correlation  decoding  of  the  PN  sequences.  For  high  and 
low  signal -to -noise  ratios  this  coding  gain  increases  linearly 
with  constraint  length  and  is  therefore  limited  only  by  the 
constraint  length. 

e.  Compared  to  the  alternative  configuration  of  direct 
sequence  spectrum  spreading  in  conjunction  with  rate  1/2,  con¬ 
straint  length  7  convolutional  codes  decoded  via  the  Viterbi 
algorithm,  it  is  clear  that  law  rate  coding  has  the  advantage 

of  requiring  only  one  level  of  implementation.  Another  advantage 
resides  in  the  fact  that  the  coding  gain  for  the  same  total 
bandwidth  expansion  can  be  much  larger  than  the  S  db  achieved  with 
rate  1/2  cades.  Finally,  we  note  the  possibility  of  better  per¬ 
formance  in  a  bursty  environment  where  Viterbi  decoding  of  rate 
1/2  codes  requires  substantial  interleaving  and  usually  results 
in  a  reduced  coding  gain. 
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CHAPTER  6 

Adaptive  Null -Steering  Arrays  Transient  BehaviOT 
Introduction 

This  Chapter  of  the  final  report  presents  a  novel  approach 
for  analysing  the  transient  behavioT  of  adaptive  array  for  both 
deterministic  and  random  jammers.  Using  this  technique  the  be¬ 
havior  of  adaptive  arrays  in  the  presence  of  step  and  periodic 
jammers  is  analyzed.  The  analysis  has  been  carried  out  for 
adaptive  nulling  antennas  which  have  stationary  or  fixed  main 
beams.  This  corresponds  to  the  tactical  scenario  of  a  line  of 
sight  link  in  a  jamming  environment. 

The  first  section  of  the  report  presents  a  review  of  the 
Widrow^  and  Applebaum^2 ’ algorithms.  These  two  algorithms 
represent  the  starting  point  of  most  modern  adaptive  theory. 

The  Widrow  or  least  mean  square  error  (LMS)  algorithm  was  de¬ 
veloped  by  Widrow  and  his  co-workers.  They  applied  methods 
used  in  adaptive  filters  to  phased  arrays.  About  the  same  time 
Applebaum  developed  a  similar  algorithm  based  on  maximizing  the 
signal  to  noise  ratio.  The  algorithm  is  sometimes  called  the 
MSN  algorithm.  Following  Widrow' s  initial  work  the  LMS  algo¬ 
rithm  was  developed  further  by  Griffiths ^  and  Frost. ^  They 
found  that  one  can  maintain  a  chosen  frequency  characteristic 
for  the  array  in  a  desired  direction  while  discriminating 
against  noises  coming  from  other  directions.  3rennan  and 
Reed^*'^  extended  the  domain  of  the  MSN  algorithm  to  include 
adaptivity  both  in  the  spatial  as  well  as  the  temporal  (Doppler 


1 

I 

4 


i 
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filter)  domain.  Recent  emphasis  in  the  field  has  been  on 
speeding  up  loop  convergence.  A  number  of  techniques  such  as 
orthogonalication,  ^  null  control  law  modification  have 
appeared  in  the  recent  literature  in  this  connection. 

The  second  section  of  the  report  is  devoted  to  a  determin¬ 
istic  analysis  of  an  analog  implementation  of  the  L>1S  algorithm. 
Here  it  is  assumed  that  deterministic  signals  are  incident  upon 
the  array  and  that  receiver  noise  can  be  neglected.  Three  par¬ 
ticular  examples  are  explored:  first,  the  case  of  a  step  plane 
wave  normally  incident  upon  the  array;  second,  a  step  jamming 
signal  incident  obliquely  on  the  array;  and  third,  a  periodic 
pulsed  jamming  signal  obliquely  incident  on  the  array.  For  all 
three  cases,  an  explicit  solution  is  obtainable. 

The  third  section  of  the  report  treats  the  wide  spectrum 
random  jammer  applied  to  an  array  utilising  the  tfidrow  algo¬ 
rithm.  The  approximate  equation  for  the  mean  is  found  exactly 
by  using  special  group  properties  of  the  governing  equations. 
This  exact  solution  leads  to  explicit  expressions  for  the  decay 
coefficients.  These  decay  coefficients  are  related  explicitly 
to  the  number  of  elements  in  the  array,  the  signal  amplitude 
and  the  jammer  noise  power.  The  case  of  a  step  jammer  is 
treated  in  detail  and  certain  conjectures  are  made  about  the 
behavior  of  the  array  when  a  periodic  jammer  is  applied. 

Adaptive  Array  Principles 

A  review  of  the  principles  of  operation  for  the  Widrow 
and  Applebaum  algorithms  will  now  be  presented.  Although  there 


exist  many  variants  o£  these  two  basic  structures,  an  under¬ 
standing  of  the  motivation  and  operation  of  these  algorithms 
is  fundamental. 

Let  us  consider  the  Widrow  algorithm  first.  Let  X^(t) , 

X,  (t) , . .  .XjjCt)  denote  the  quadrature  outputs  of  an  N/2  ele¬ 
ment  phased  array.  They  are  assumed  to  be  random  functions  of 
time.  Suppose  these  functions  are  multiplied  by  weighting 
factors  ,W. , . . .W^  and  added  to  form  the  output 

N 

SCt)  -  l  W.X.ft)  (6.1) 

1-1  *  * 

as  indicated  in  Figure  6.1. 

Next  S(t)  is  compared  to  a  reference  function  R(t)  which 
may  be  taken  as  the  "desired”  response)  with  the  object  of  mini¬ 
mizing  the  error  e(t),  or  more  precisely,  its  mean  squared,  viz, 

<| cCt) | 2>  (6.2) 

Thus  one  would  like  to  adjust  the  weights  ,W2 ,W_ , . . .W^  such 
that 

2  N  2 

<!e(t)  r>  -  <|RCt)  *  I  W.X.(t)r>  (6.3) 

1-1  4  1 

is  minimized.  Both  X^(t)  and  the  reference  function  R(t)  (a 
stochastic  process)  are  assumed  real.  The  case  of  complex 
functions  will  be  considered  later.  Thus  (the  symbol  £  will' 
stand  for  mean  squared  error) . 

£  =  <ie(t) |2>  =  <e(t)2> 


mo 


'«***«*.  ay- 


iT'jw 


,  N 

-  <R'(t)>  -  2<R(t)  I  W  X,(t)>  *  l  W  W  <X ,(t)X  ft)> 

l-l  1  i,a  m  m 

Note  that  the  weights  W£  are  not  functions  of  time  and  are  taken 
as  nonrandom  parameters.  Differentiating  e  with  respect  to 
and  setting  each  derivative  to  zero,  yields 


N 


°  •  !§_  *  ■  2  <RCt)XnCt)>  ♦  2  i  w£  <  xtCt)XnCt3> 


n 

n-1,2, . . .  ,N 
Now  define  column  matrices 


(6.4) 


'<R(t)X1(t)>- 

rWloptl 

31 

w 

w_ 

.<R(t)XN(t)>. 

—opt  - 

2opt 

• 

NOpt 

(6.5) 


where  W£0pt*  W£  that  satisfies  (6.4)  £«1 ,2 , . . . ,N. 
Similarly  define  the  matrix 


C  •  <X(t)X*(t)>  .  [<x1(t)XB(t)>] 


(6.6) 


with  X(t) 


Xi(t). 


Xjj(t) 


and  Xr  denoting  the  transpose  of  X.  The 


solution  of  (6.4)  then  becomes 


2cPt  "  c  £ 


(6.-) 
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These  are 


Equation  (6.7)  gives  the  weights  that  minimize  £. 
expressed  in  terns  of  the  covariance  natrix  C  that  depends  on 
the  statistics  (second  order)  of  the  N  stochastic  inputs  and 
the  column  matrix  (vector)  S  that  involves  the  cross  correlations 
between  the  inputs  and  the  desired  response  R(t) .  Thus,  to  find 
the  required  weights  one  must  know  the  covariance  matrix  and  the 
statistical  dependence  of  R(t)  on  X^t).  In  many  applications 
the  covariance  matrix  of  X^Ctps  would  not  be  known  and  the  ex¬ 
plicit  inversion  as  indicated  by  (6.7)  could  therefore  not  be 
carried  out.  An  adaptive  array  carries  this  inversion  process 
out  automatically  by  sensing  the  error  e(t)  and  feeding  it  back 
to  adjust  the  weights  until  5  is  minimized,  i.e.,  (6.7)  is 
satisfied.  A  possible  closed  loop  system  (adaptive  filter)  for 

achieving  this  is  shown  in  Figure  6.2.  Here  the  symbols  <g)  and 

A 

/ f\  denote  a  multiplier  and  integrator,  respectively. 

Consider  a  typical  loop.  Performing  the  indicated  opera¬ 
tions  yields 


(t) 

ZK  Xt(t)  e(t)  - 


N 


Since  s(t)  *  R(t)  -  £  W4(t)  X^t)  ,  one  finds 


2*1 


dw.ct) 


N 


-  -ZX  XACt)  I  \(t)  XB(t)  *  2X  R(t)  x4(t) 

m*l 

or,  in  matrix  form 

-  -  2X  X(t)Xt(t)  W(t)  -  2X  R(t)  X(z) 


(5.31 
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This  is  a  stochastic  differential  equation  for  the  vector  W(t) . 

It  can  be  solved  rigorously  only  if  the  detailed  statistics  of 
X(t)  and  R(t)  are  known.  However,  if  one  assumes  that  the  loop 
response  time  is  much  longer  than  the  decorrelation  time 
of  X(t)  (in  other  words,  the  X2(t)  vary  much  more  rapidly  than 
the  time  constants  of  the  adaptive  process  of  adjusting  the 
weights) ,  then  one  has 

<X(t)  XC(t)  W(t)>*<X(t)  XC(t)><W(t)>  ,  (5.9 

i.e.,  W(t)  remains  essentially  constant  during  the  averaging 
operation. 

Taking  averages  of  both  sides  of  (6.8)  yields 
d<W(t)> 

— 3—  -  -2K<X(t)  Xi:Ct)><W(t)>  +  2K<R(t)  X(t)>  (5. 1C 

If  the  convariance  matrix  C  *  <X(t)  Xt(t)>  is  positive  definite 
then  (6.10)  is  stable  and  in  the  steady  state  (i.e.,  t~«) 

d<W(t)> 

— ~ 0 

at 

Thus,  one  finds  that  as  t—» 

-2K<X(t)  Xt(t)><W(t)>  *  2K<R(t)  X(t)  >-*0 , 
or. 

lim  <W(t)>  -  W  , 

t—  * 

as  given  by  (6.7).  Thus,  subject  to  the  requirement  that  C  be 
positive  definite  and  that  K>0,  the  adaptive  filter  yields  (5.7) 
as  the  steady  state  solution.  The  settling  time  (i.e.,  the  time 


122 


r 


to  reach  steady  state)  is  determined  by  K  and  the  eigenvalues  of 
C.  It  will  be  convenient  to  examine  the  solution  of  (6.10)  with 
the  aid  of  Laplace  transforms.  Let  £(s)  be  the  Laplace  transform 
of  W(t) .  Then  (6.10)  is  equivalent  to 

-W(Q*)  *  s£(s)  +  2KC  tf(s)  •  —  £  (6.11) 

Let 

W' (s)  *  TW(s),  where  T  diagonalizes  C  and  TTt«I.  Then 
(6.11)  may  be  written 

-W(0  +  )  ♦  STW*  (s)  ♦  2KT1  A  TW(s)  » 

or 

-TW (0*)  ♦  SB'  (s)  ♦  21  A  tf’ (s)  «  2K  TS/s  (6.12) 
where  A  «  diag.  (X^,  Xj. . .  X^) . 

Solving  (6.12)  for  UM (s) , 

W'(S)  «  (s  +  2K  A)*1  TW(0+)  ♦  fi(s  «•  2K  A)*1  TS  (6.13) 

Inverting  the  transform  in  (6.13) 


-  TT 

ro 

i- 

*  0 

* .-2KX,t 

e  * 

• 

TW  (0+) 

••x  Ls 

X  ;  lt  ’ 

Az.  ; 

1 

*  *  i 

JL  - 

3  *— i 

JL  ’•  — 

To  find  W(t) ,  the  above  is  multiplied  by  Tr  from  the  left  to  yield 


w(t)-r 


o  1 

-  ZktX 


1  TW (0*)  -T1 


- 2KtX  „  :  , 

t  K  !  TS  «•  C  XS  (6.14) 
AS,.  i  “ 
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Evidently,  the  time  constant  of  the  system  is  given  by 


M 


(6.15) 


where  X  .  „  is  the  smallest  eigenvalue  of  C. 
mm  * 

The  preceding  discussion  followed  closely  the  adaptation 
process  described  by  Widrow. ^  The  implementation  is  sometimes 
referred  to  as  the  LMS  algorithm.  A  somewhat  different  phrasing 
and  implementation  of  the  adaptive  process  is  due  to  Applebaum. 

It  is  most  appropriate  for  application  to  array  antennas.  First, 
consider  an  N-element  array  with  outputs  at  the  array  ports  u^, 
u, ,...u^..  The  outputs  at  the  array  elements  are  complex  time 
functions  corresponding  to  the  envelope  in  the  analytic  signal 
representation  of  the  composite  signal  (including  modulation)  and 
noise.  When  no  noise  (generalized  interference  which  may  be 
ambient  noise  together  with  jamming)  is  present 


uk  -  a  Sk  (6.16) 

where  Sk  are  weighted  element  patterns  as  they  usually  appear 

i  k  *  D 

in  the  array  factor  (e.g.,  Sk  «  mk  eJ-  — k)  while  a  ■  o(t)  is 
the  information  bearing  signal.  When  only  noise  (interference) 
is  present  at  the  array  output  ports 

uk  s  nk(t)  (6.17) 

which  may  be  composed  of  the  ambient  noise  n  .  (t)  and  jammer 

O  f  & 

noise  nj  k(t) ,  vis. 

uk  ■  "o,k(t)  *  nJ,kft)  (6-l8) 
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Suppose  that  in  the  "quiescent"  environment  (rij  ^(t)  ■  0)  the 
weights  for  optimum  signal  reception  are  W,  ,  or 


then  upon  summation  the  desired  signal  is  W^_t,  i.e.,  t»s_  is 
the  vector  that  determines  the  correct  array  illumination  under 
quiescent  conditions.  When  noise  only  (interference)  is  present, 
the  power  output  after  array  output  summation  is 


<|ntW|2>  *  w*  <n*  n 5  w  •  W*  C  W 


(6. IS 


where  C  is  the  covariance  matrix  of  interference,  A  is  the 

* 

conjugate  to  transpose  of  A,  and  p  is  the  complex  conjugate  of  p, 
If  the  quiescent  noise  and  jammer  noise  are  uncorrelated 


C  *  Co  *  CJ 


The  signal  to  "noise"  ratio  can  now  be  defined  as 


2>- I'llii  2 


W  Of 


(6.20) 
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and  one  would  like  to  choose  W  so  that  is  maximum.  Writing 


|wrti2  -  w+tVt  -  w+(t*tr)w 

+  *  f 

W  (t  tZ)W 

I  "  <ifll ->  '  w+cw 

One  readily  finds  that  the  preceding  is  maximized  by 

W  -  W0TJt  -  C~lt*  (6.21) 

The  adaptive  implementation  that  realizes  this  optimum  is  shown 
in  Figure  6.4  (for  the  kth  element). 

Let  RC  ■  t  be  the  time  constant  of  the  integrator.  Then 


-uk  jx  +  zl 


1  dwk 

z  (Tat~  *  V  * 


Taking  averages  and  reverting  to  matrices 

-  cwct)  *  i  •  z  Ct  zr  *  ^Ct)) 

or 

dW  * 

Th=  +  W  -  -G  O  W+  Gt  (6.22) 

The  steady  state  solution  is  obtained  by  setting  the  time 
derivative  to  zero.  Hence 
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?ijur*  6-4.  Ispiestaiaticn  of  Applao aum  Algorithm 


-  123  - 


.  «*.  - -a.  _  -  ■«■ 


(6.23 


[G  +  I  £•]  W  ■  t* 

.  *.W  «  [C  +  |  T]"1  t* 

For  sufficiently  large  gain,  W  Wg  .  However,  when  G  is 
finite,  there  is  always  a  residual  error. 

Deterministic  Reponse  of  Adaptive  Arrays 

In  this  section,  the  adaptive  array  will  be  examined  from 
a  deterministic  point  of  view.  Receiver  noise  will  be  neglected 
and  any  signals  incident  upon  the  array  will  be  assumed  to  be 
deterministic.  Although  an  oversimplification  of  the  actual 
situation,  the  deterministic  behavior  of  a  device  many  times 
supplies  the  physical  insight  necessary  to  extract  information 
from  the  more  complex  stochastic  case. 

Our  discussion  will  be  limited  to  the  Widrow  algorithm 
shown  in  Figure  6.2.  There  an  M  port  array  having  N-M/2  in  phase 
and  N  quadrature  ports  is  considered.  The  transient  behavior 
of  the  array  will  be  examined  in  three  situations  of  interest. 

First,  the  behavior  of  the  array  will  be  examined  when  a 
step  plane  wave  is  normally  incident  upon  the  array.  Since  the 
Widrow  algorithm  that  has  been  chosen  contains  no  steering 
mechanism  its  optimum  direction  of  reception  is  normal  to  the 
array.  Thus  once  the  initial  transient  has  decayed,  the  array 
should  be  locked  on  the  normally  incident  plane  wave. 

Proceeding  with  the  analysis,  the  equation  fOT  the  weights 
(6. S)  is  given  by 
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dW(t) 

“at 


-  2KX(t)  XC  (t)  £(t)  *  2XR(t)X(t) 


2LCo)  *  Wq 


(6.24 


where  the  initial  weights  ^  have  been  arbitrarily  chosen. 
Assume  the  input  signal  is  of  the  form 


X(t)  -  X  u(t)  ,  X  *  A  £  (6.23) 

where  A  is  the  amplitude  of  the  plane  wave  normally  incident  on 
the  array,  a  is-  an  array  vector  given  by 


a 


2N  »  M 


(6.26) 


where  ljj  and  0^  are  column  vectors  of  l's  and  0's  respectively 
and  u(t)  is  the  Heavside  function.  The  array  vector  takes  the 
form  of  (6.26)  since  the  incident  plane  wave  is  assumed  to  excited 
only  in-phase  components.  Finally  the  reference  signal  is  chosen 
as  R(t)  »  A. 

Using  the  above  information  in  (6.24) 


dW(t)  . 

-33 —  *  -2XX  X  W(t)  +  2KAX  t  >  0  (6.27) 


W(o)  -  ^ 

tH 

This  is  simply  a  2N  1  order  system  of  ordinary  differential 
equations  (o.d.e.'s)  with  constant  coefficients .  A  solution 
can  be  written  as 

♦» 

W(t)  ■  2XA  f  exp [ 2XX  XC ( t ' - t) ] dt '  (5.23, 

0 

♦  exp[-2XX  XZtl^  ;  t  >  0 
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The  expression  for  W(t)  given  above  can  be  substantially  simpli¬ 
fied  by  making  use  of  the  following  identity: 


exp  a  X  X 1  ■  I  +  (exp  a  j  X 1 2  -1)  — 


X  X 


(6.29) 


2  t 

where  |X|  *  XX,  I  is  the  identity  matrix  and  a  is  a  constant. 

This  identity  is  obtained  in  Appendix  A. 

Using  (6.29)  in  (6.28)  with  | X | 2  »  A2N,  we  have 


W(t)  «  A(l-e‘2KNA“t)|  ♦  ^ 


?riJ,2r  X  X1- 
+  (e-2KNA  t.1)=fj=_ 


t  >  0 


(6.50) 


Since  an  explicit  expression  for  W(t)  has  been  obtained  the  out¬ 
put  of  the  adaptive  array,  y(t)  can  be  calculated 


y(t) 


y(t)  -  xrw(t) 

A[1.e-2KNA2t)  .  xt  -20A2t 


(6.31) 


In  the  steady  state,  the  desired  output  is  , 


lim  y(t)  ■  A 

? 

In  fact,  the  transient  decays  with  a  time  constant  t  *  1/2KNA*’. 
Thus,  the  larger  the  factor  KNA2  becomes  the  faster,  the  array 
locks  unto  the  desired  signal.  We  note  that  if  ■  X/N  then 
y(t)  *  A  for  t>0.  This  essentially  means  the  array  was  in  steady 
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state  at  t*G  and  no  transient  was  needed  to  natch  the  initial 
conditions  to  the  steady  state. 

Next  an  analysis  will  be  made  o£  a  step  jammer  obliquely 
incident  on  an  array.  We  will  assume  that  the  array  is  in  steady 
state  before  the  jamming  signal  intercepts  the  antenna,  let 


X(t)  •  X  ♦  XjU(t) 


(6.32) 


where  X  is  given  in  (6.2S)  and 


h  m  B£j  »  £j 


^JC 


Hjs 


(6.32a) 


and 


jaT  2jaT  NjaT  * 

j£jS  -  («  J . e  V  ,  aj  -  ^sin3J  (6.'') 


2JC  * 


Here  d  is  the  array  element  spacing,  X  is  the  wavelength,  3j 
is  the  angle  of  incidence  measured  from  the  normal  to  the  array 
and  3  is  amplitude  of  the  jamming  signal. 

Using  this  in  (6.21)  with  R(t)*A,  we  have 


dW(t)  ... 

—  -  -2XX.  X^W(t)  *  2m  (6.34) 

X 

W(o)  *  f  t>0 

where 

X  •  X  ^  Xj  (6.33) 

The  initial  setting  for  the  weights  has  been  chosen  so  that 
the  array  will  be  in  steady  state  (y*A)  at  t*0.  When  the 
methods  used  in  the  previous  section  are  applied  to  this  case, 
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we  f ind 


W(t) 


ACl-e 


- 2K  j  X j 


[I  «■  (e' 


•zKixrt 


1) 


X  XZ  X 

TTT2  i 


t>0 


A  direct  calculation  gives  the  output  o£  the  array  once  the 
weight  vector  is  known.  One  finds 


ytt)  •  X1  WCt) 


v(t)  .  A .  ^ 


.where 


J  *  2j  £  *  £*  £j 


The  decay  constant  t  is: 


Since 


T  •  [2KIXI2]'1 


!  X  l  2  •  |X  ♦  Xrl2  -  |  Ac  +  BoT|2 


12.  SJI 


«  (A2  *  B2)N  ♦  2ABJ] " 1 

we  have 

t  «  [ (A2  +  B2)N  ♦  2AB J] " 1 
For  large  B,  we  find 


Thus  we  see  from  (6.37)  and  (6.41)  that  although  the  transient 
has  a  higher  amplitude,  the  larger  the  jamming  signal  be- 


(6.36) 

(6.37) 

(6.38) 

(6.39) 

(6.40) 

(6.41) 
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comes,  the  decay  rate  also  becomes  much  more  rapid.  For 
large  N,  the  behavior  is  more  complex  since  J  depends  on  N’ 
and  can  be  negative. 

Before  preceding  on  to  the  periodic  case,  we  shall  point 
out  an  important  weakness  of  the  deterministic  approach  to  the 
problem.  We  have  seen  from  (6.37)  that  as  t  becomes  large  the 
output  of  the  array  approaches  the  required  steady  state,  how¬ 
ever,  an  examination  shows  that  no  null  has  developed  in  the 
array  in  the  direction  of  the  jamming  signal. 

This  anomalous  behavior  in  the  deterministic  case  leads  us 
to  the  conclusion  that  the  adaptive  processor  has  taken  advan¬ 
tage  of  the  special  properties  of  the  input  signals  to  the  array 
to  cancel  out  the  jammer  signal.  In  general,  one  would  suspect 
that  in  order  for  the  array  to  cancel  out  a  random  jamming 
signal,  a  null  would  have  to  be  formed  in  the  direction  of  the 
jamming  signal;  thus,  we  are  led  to  the  conclusion  that  cer¬ 
tain  characteristics  of  the  array  could  be  different  in  the 
random  case.  This  is  indeed  the  case  as  we  shall  show  later. 

Before  treating  the  random  case,  we  shall  consider  the 
deterministic  steady  state  response  of  the  array  to  a  periodic 
jammer.  Although  the  method  suffers  from  the  previously  men¬ 
tioned  weakness,  the  work  serves  to  indicate  methodology  and 
gives  insight  into  certain  aspects  of  the  arrays  behavior  in 
the  presence  of  periodic  signals. 

To  proceed,  we  take  X(t)  to  be 


where  S(t)  is  a  periodic  function  of  period  7  consisting  of 
pulses  of  unit  height  and  duration  t  as  shown  in  Figure  6.3. 
We  now  look  for  solutions  to  (6*24)  that  are  periodic, 

i.e. , 

W(t+T)  •  W(t) 


Our  technique  is  as  follows.  First,  assume  W(t)  is  known  at 
t»0,  i.e.,  W(o)  ■  W^.  Now  solve  {6.24)  in  the  region  of  0<t<t. 
Next,  use  the  continuity  of  W  at  t»t  to  establish  initial 
conditions  for  the  region  r<t<T.  Now  solve  (6.24)  in  this 
region  and  set  W(T)  *  W^.  This  gives  an  algebraic  system  of 
equations  for  W  .  Upon  execution  of  the  foregoing,  we  find 


W(t)  -  A(l-e 


-2K 


ii>S 


“1]  w, 
x  z 


C< 


0<  t<T 


and 


W(t)  »  A(l-e'2KNA  ♦  [I  +  (e 


-2KNA4(t-t) 


i<t<T 


The  equation  for  is  given  by 


Ho  •  Apo  I  *  ti-po  T-)tHo  *  P17T72 CA  -  i*So)l 

X  I 


(( 


where 


p  .  i.e" 2KNA" (T-t) 
o 


?i  .  i.e^KSXl  “T 


.43) 


(6.44) 


•45) 
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For  arbitrary  N,  the  system  of  equations  for  is  dif¬ 
ficult  to  solve  explicitly.  Since  only  y(t)  is  required,  the 
complete  solution  to  (6.45)  is  not  required  but  only  the  pro- 

A  A  ** 

jection  of  in  the  X  and  X  directions,  i.e.,  Xx  W  and 
•'  — o  —  —  —  — o 

Xt  W  .  Let  us  define 
—  —o 

Y*  -  Xr  W  (6.40 

o  —  — o 


■  r  »<,  (6-*7> 

where  Y*  are  scalar  quantities  which  determine  v(t).  This  can  be 
seen  from  the  fact  that 


y(t) 


X1  W(t)  ,  0<t<T 

X*  W(t)  ,  r<t<T 


(6.48) 


If  (6.43)  and  (6.44)  are  used  in  (6.48),  we  obtain 


y(t) 


1  +  (Y*  -  l)e'2K^  x  ,  0<t<T 


1  +  CY'  -  !)e+2KN’A^(T-t) ^  „<t<T 


(6.49) 


Equations  for  Y*  can  be  obtained  by  multiplying  (6.45)  from  the 
left  by  X*  and  Xt.  We  find 


?  X*  X 

“^”"5  =Yo  *  pCp, 


'*r  T  ~ 

X1  X  X1  X 


N  jX 


1—  -V  Yo  *  8] 


(6.50) 


a  A 

X*"  X 

Po  Y0  +  PCP.-D-T'--,"  y!  -  g, 

lit 


(6.51) 
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where 


t  r 

p  r  xr  x 


M  I  X.I 


P.  ?  X1  X 

Sr»  *  -M - 777 - 

|X|‘ 


X*  X 


The  solution  of  this  system  (6.50-51)  for  Y“  is  straight  for¬ 
ward  and  yields  Y*  *  1.  3y  using  this  result  in  (6.4.9),  we 
find  that 


y(t)  -  i  ,  o<t<i 


(6.32) 


The  physical  interpretation  of  this  result  is  that  the 
adaptive  array  has  had  an  infinite  amount  of  time  to  adapt  to 
the  presence  of  a  periodic  jammer,  and  as  in  the  case  of  the 
step  jammer,  it  adapts  perfectly*  This  is  an  important  re¬ 
sult.  It  indicates  that  in  order  to  defeat  an  adaptive  array, 
the  jamming  signal  must  be  constantly  changing  in  a  random 
fashion. 


Random  Response  of  Adaptive  Arrays 

Before  concluding  this  Chapter,  we  would  like  to  present 
an  interesting  analysis  of  the  approximate  mean  equation  given 
in  (6.10).  Although,  the  transient  character  of  this  equation 
has  been  treated  in  the  literature  and  reviewed  in  this  report 
(6.11-6.15),  many  of  the  symmetry  properties  of  the  equation 
were  not  adequately  utilised.  An  indication  of  this  is 
illustrated  by  (6.13).  Here  the  transient  properties  of  the 
array  are  given  in  terms  of  the  minimum  eigenvalue,  \  .  , 
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however,  no  information  about  the  behavior  of  X  .  with  resoect 

mm 

to  the  system  parameters  such  as  N ,  A,  etc.,  is  given. 

We  shall  rectify  the  situation  by  utilizing  the  group 
properties  of  the  system  of  equations  of  the  weights  W(t) .  It 
is  shown  in  Appendix  A  that  by  proper  parameterization  of  the 
problem,  the  2N  order  system  of  equations  can  be  replaced  by  a 
fifth  order  system  for  any  N.  This  system  is  solved  explicitly 
in  the  constant  coefficient  case,  and  a  generalization  of  the 
identity  used  in  the  deterministic  case  is  derived.  We  use 
this  identity  in  this  section  to  obtain  our  main  result. 

Although  the  method  we  employ  has  broad  application,  we  use 
it  just  to  analyze  the  case  of  a  broad  spectrum  step  jammer. 

The  analysis  proceeds  as  follows.  Consider  (6.10)  with 
R(t)  •  A: 


-  -2K  <X(t)  Xt(t)><W(t)>  ♦  2KA  <X(t)>  (6. S3) 

W(o)  -  W 
—  — o 

where 

X(t)  -A a  *  Nj(t)  £j  (6.54) 

Here  a  and  Oj  are  array  factors  as  defined  in  (6.26)  and 
(6.32a),  respectively  and  Nj(t)  is  a  zero  mean  stationary 
Gaussian  processes  with  spectral  density  Nqj. 

Using  this  information,  we  have 

<X(t)  X  ( t )  > " A  £  £r  +  N'oJ  £j  £j  ,b.o5) 


139 


and 


<X(t) >  *  A £  (6 . 5o, 

Putting  (5.53)  and  (6.56)  into  (6.33)  an  equation  with  con¬ 
stant  coefficients  results.  The  solution  to  this  equation 
is 

<W(t)>  *  2KA  I  gQCt’  "c)  dt,  ^  +  (6.37) 

^0 


where 

Q  •  2XtA2  a  »'  .  MoJ  ^  a]]  (6.53) 

The  expression  in  (6.37)  can  be  substantially-  simplified  by- 
using-  the  identity  given  in  Appendix  A.  We  have 

eQt  -  f0(2)  ♦  f^(z)eX+t  ♦  f_(t)eX  1  (6.59) 


where 


(a+b)N  ±/(a-b)Z^NZ  »  labJZ 


(6.60) 


with 

a  -  ZKA2  ,  b  ■  2KNoJ  (6.61) 

The  vector  z_  appearing  in  (6.59)  consists  of  five  components. 

Each  component  is  a  matrix  s  These  matrices  form  a 

complete  set  of  generators  for  a  group  under  matrix  multi¬ 
plication.  It  is  the  realization  that  expQt  is  also  a  member 
of  this  group  that  leads  to  the  identity  (6.39).  The  explicit 
formula  for  z^,  i*0,...i  and  the  functions  fQ  and  f _  are  given 
in  Appendix  A. 

The  important  property  of  (6.39)  to  notice  is  that  exp(Qt) 
has  been  represented  in  terms  of  three  scalar  functions  of 


-  U0  - 


time  whose  exponential  decay  rates  are  given  explicitly  in 
terms  of  the  parameters  of  the  system.  Using  (6.39)  in 
(6.37)  and  integrating,  we  obtain: 


1  -  e” ^  ^ 

<W(t)>  «  2KA[  ■■  -  f+(2).o 

X 

1  -  e  *  ^  ^ 

♦  i - S -  f. (t)-£] 

X 

where  we  have  used  the  fact  that  f0(£)-£  *  0. 

Now  that  we  have  the  exact  solution  to  the  approximate 
mean  equation,  we  will  examine  its  behavior  in  the  steady 
state.  Let 


(6.62) 


<vto>  -  lia  <W(t)> 


By  letting  t  approach  infinity  in  (6.62),  we  find 

f+(£)  £.(£) 

<*»>  -  2Ka[ — ♦  — r— ]  £ 

X  X 


(6.63) 


(6.64) 


Using  the  expressions  for  fQ  and  f+  given  in  Appendix  A  to 
simplify  (6.64),  we  obtain 


<w„> 


(No  -  j£j) 

V”  J2) 


(6. 63) 


Since  the  output  of  the  adaptive  array  is  given  by  y(t)  * 

Xx(t)  W(t) ,  we  have  for  t-<-»: 

<yae>  ■  lia  <y(t)>  *  A  £r  <WflB>  (6.66) 

A(N  S*  £  -  J  £l  £j) 
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m 


(6.67) 


<y»>  *  a 


Thus  we  see  Chat  the  j  asm  lag  signal  has  been  completely 
cancelled  out  in  the  steady  state.  One  may  recall  that  this  is 
exactly  the  result  obtained  in  the  deterministic  case. 

Next,  we  would  like  to  determine  if  the  cancellation  of 
the  jamming  signal  was  due  to  a  null  in  the  phased  array  in 

i 

the  direction  of  the  jammer.  To  do  this,  we  define  the  array 
amplitude  directivity  factor  F(3,t)  by 


where 


and 


F(8,t)  -  oj  <W(t)> 


(6.68) 


(6.69) 


Sec 


£es 


(e 


jo  -Zja 


TNja, 


(6.70) 


a 


2*d 

-r 


sin8 


We  note  that  the  directivity  factor  is  a  function  of  time 
since  as  the  weights  change,  the  radiation  pattern  of  the 
antenna  changes. 

We  would  now  like  to  compare  F(o)  to  ?(3j)  for  t— .  Let 

F„(8)  -  lia  F(3,t)  -  cl  <Wso> 

t— 


and  then,  using  (6.63)  we  have 

F (o)  ■  1 


(6.71) 


and 


F(6j)  -  0  . 


Thus  we  conclude,  that  an  infinitely  deep  null  has  been  de¬ 
veloped  in  the  direction  of  the  jamming  signal.  Thus  we 
see  that  in  the  random  case  the  array  factor  is  responsible 
for  nulling  out  the  jamming  signal  rather  then  the  signal  can¬ 
cellation  that  occurred  in  the  deterministic  case.  The  array 
behavior  in  the  random  case  requires  just  one  set  of  steady 
state  weights  to  cancel  out  any  jamming  signal  entering  at 
3  *  8j  while  the  array  behavior  in  the  deterministic  case  re¬ 
quires  a  different  set  of  weights  for  each  different  jamming 
waveform  at  S  ■  6j. 

Before  concluding  this  section  we  would  like  to  discuss 
the  settling  time  of  the  array  in  the  random  case.  An  examin¬ 
ation  of  (6.62)  shows  that  there  are  two  decay  factors  A  +  and 
X"  given  in  (6.60).  This  is  in  contrast  to  the  step  jammer 
considered  in  the  deterministic  case.  In  that  case  only  one 
decay  factor  appeared. 

The  time  constant  of  the  array  to  reach  steady  state  is 
given  by 


T 


sJ 


1 

min(A+,x”) 


(6.72) 


where 


min (A  ,  A*} 


X*,A’>A+ 

*  ♦  m 

A  y  A.  >4 
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For  large  jamming  amplitudes  b>>a  or  Nqj>>A~,  expression 
(6.60)  can  be  substantially  simplified.  We  find 

X+  *  bN 

and 

A."  «  J  C N2  -  J2)  >  0 

This  shows  that  the  \+>  is  the  counterpart  to  the  decay  fac¬ 
tor  found  in  the  deterministic  case  (6.41)  while  X*  represents 
a  new  effect.  When  b  becomes  large  it  is  clear  that  for  fixed 
M,  A  becomes  the  slowest  decay  constant  and  thus  determines 
the  transient  time  of  the  array. 

It  is  generally  believed  that  the  larger  the  jamming 
signal,  the  faster  the  array  will  lock  in.  It  is  clear  that 
the  presence  of  X  will  tend  to  mask  the  above  effect  for  very 
large  b  since  the  array  settling  time  will  depend  only  on  X* 
and  this  parameter  is  independent  of  b .  Thus  there  appears  to 
be  a  saturation  effect.  As  b  is  increased,  the  transient  time 
decreases  until  approximately  A+  >  A*  and  then  the  transient 
time  remains  relatively  fixed. 

We  point  out  that  the  above  results  are  of  a  preliminary 

i 

nature.  Only  the  approximate  equation  for  the  mean  has  been 
solved, not  the  exact  mean  equation.  In  addition,  no  check  of 
the  variance  has  been  made  which  is  essential.  Work  is  con¬ 
tinuing  on  these  topics. 

The  application  of  the  techniques  used  in  this  section, 
to  other  types  of  jammers  is  certainly  possible.  In  particular, 
the  investigation  ot  a  periodic  wide  spectrum  jammer  could  be 
treated  by  these  methods.  Although  the  problem  has  not  been 


attacked  formally,  it  is  our  conjecture  that  the  steady  state 
array  factor  would  have  a  perfect  null  in  the  direction  of  the 
jammer.  This  would  indicate  that  to  effectively  jam  an  array 
the  use  of  aperiodic  modulation  of  the  wide  spectrum  jamming 
signal  would  be  more  effective  than  a  periodic  modulation. 


Conclusion 


We  have  shown  that  the  methods  are  at  hand  for  assessing 
the  vulnerability  of  adaptive  null  steering  antenna  arrays  to 
broadside  pulse,  step  CW,  periodic  and  random  jammers.  A  re¬ 
markable  result  is  that  for  periodic  jammers  and  for  noise  jam¬ 
mers  the  steady  state  null  can  be  perfect  (assuming  a  perfect 
reference).  This  leads  us  to  the  conclusion  that  null  steering 
antennas  would  be  most  vulnerable  to  blinking,  non  periodic, 
wide  spectrum  jammers.  For  this  situation  it  is  important  to 
understand  the  parameters  that  determine  the  array  time  constants. 
This  we  have  arrived  at  explicitly  for  both  deterministic  and 
random  signals.  A  further  interesting  aspect  of  the  analysis 
is  that  by  exploiting  the  symmetries  in  the  array  equations  we 
are  able  to  show  that  the  array  settling  time  constants  are  de¬ 
termined  by  three  scalar  functions  of  time  whose  exponential 
settling  time  are  given  explicitly  in  terms  of  the  parameters 
of  the  system.  This  is  true  for  any  sice  array.  This  result 

i 

will  be  of  enormous  benefit  in  the  design  and  implementation  of 
adaptive  arrays.  Adaptive  arrays  are  the  only  ECCM  technique 


which  allows  substantial  anti- jam  performance  without  bandwidth 


expans  ion. 
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Appendix  to  Chapter  6. 


In  this  Appendix  we  will  derive  an  identity  that  will  be 
useful  for  simplifying  expressions  of  the  form  exp{aoor  +  bcj£^)t) 
where  o  and  Oj  have  been  defined  in  (6.26)  and  (6.32a),  respectively 
and  a  and  b  are  arbitrary  constants. 

We  start  by  defining  the  matrices  zi ,  i*0,l,...,4.  They  are 


“0 


I 


zi  •  £  lz  "4  £j  £j 

"2  *  £  £j 


(A6-1) 


Next  we  recognize  that  exp{azQ  +  bz4)t}  is  the  fundamental  solu¬ 
tion  to 


-  [azQ  ♦  bz4]T(t)  (A6-2) 

T(0)  •  I 

Following  this  we  expand  T  in  terms  of  the  z..  We  have 

4 

T(t)  -  I  a.(t)  z.  (A6- 3) 

i-0  1  1 


Using  (A6-3)  in  (A6-2) ,  we  obtain  the  following  system  of 
equations  for  <a(t): 


where 


(A6-4) 
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A  * 


0  0 
a  aN 
0  0 
0  bJ 
b  0 


0  0 
0  aj  0 
aN  0  aJ 
0  bN  0 
bJ  0  bN 


The  eigenvalues  of  A,  i.e.,  and  J  *  £C  £j* 


(A6-5) 


det|A-  I  A |  ■  0  (A6-6) 

are:  A  ■  0  of  multiplicity  1,  and  X"  each  of  multiplicity  2 
where 


X£  «  Ca^bl  N  t  /(a-b)  “  NZ  *  4  abJ“  ;  Q 

One  should  note  that  although  A*  has  a  multiplicity  of  two 
no  growing  terms  in  t  result  since  <ach  of  the  two  eigenvalues 
for  A*  have  independent  subspaces.  The  same  is  true  for  x". 

The  matrix  A  can  be  diagonalized  by  the  similarity  trans¬ 
formation  T  as  follows: 


AT  »  TA 


CA6-3) 


where 

A  -  diagfO,  A+,  A+,  A^*)  (A6-9) 

and  T  is  a  3  x  S  consisting  of  the  eigenvectors  of  A.  If  we 
define 


a(t)  »  T  eft) 


-  U3  - 


(A6-10 


then  (A6-4)  becomes 


de 

at 


Ae 


CA6-11) 


eCt»  -  T'1  Oq 

The  desired  identity  can  now  be  obtained:  first,  solve  (A6-11) 
for  e(t),  find  a(t)  from  (A6-10)then  use  CA6-3). 

The  result  is 


where 


'  'D2llt  -  £0C=)  •  f.Cj).***  *  t.Cs).*** 


i.  9  (2o  *  21 »  ‘  "  ' 


and 


f  (2)  -  [(J2-N2)  z0  ♦  NSj  *  Jz2  -  Jz3  +  Nz4]/CJ'i-N4') 


—  \T*»  ^ 


f^Cz.)  *  [aJZj.  +  (X"  -  aN)  z^j/aJN 
£.(2)  *  {-  jj^CaJ  zx  CX’  -aN)  z.) 


?  -1 
1  - 


♦  aJ  z,  +  (X*  -aN)  z4)  •  [aNCJ*  -  N  )] 

Xt  are  given  by  (A6-10)  and  J  ■  1*  £j  • 

The  formula  simplifies  substantially  when  b*0.  We  have 


(A6-12) 


(A6-13) 


CA6-14) 

(A6-15) 


ae  3 

6  — “  " 1 * 


e^-l  o  aZ 
T - H 


CA6-16) 
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CHAPTER  7 


Modulation  and  Spreading  Techniques 

In  conjunction  with  the  dynamic  allocation  of  channel 
capacity  in  Chapter  3  of  this  report  the  use  of  an  appropriate 
combination  of  modulation/ spreading  technique  will  ensure  that 
the  system  will  accrue  the  maximum  advantage  with  respect  to  a 
jammer.  This  latter  topic,  namely  which  modulation/spreading 
performs  the  best  when  used  over  a  channel  with  a  fixed  bandwidth 
constraint,  was  studied  in  depth  in  a  companion  study  performed 
by  S  Consulting  Services  entitled  "Adaptive  Techniques  in  Multi¬ 
channel  Transmission."  References  [7.1]  and  [7.2]  document  both 
the  analyses  and  the  results.  Because  these  results  were  documented 
in  the  above  references  quite  extensively,  they  will  not  be 
repeated  here.  Rather  an  overview  of  the  key  conclusions  will 
be  presented,  specifically  with  respect  to  theiT  relevance  to 
an  overall  system  employing  the  dynamic  allocation  of  channel 
capacity  algorithm. 

Two  basic  combinations  of  modulation/spreading  were  considered, 
the  first  being  a  coherent  direct  sequence  CDS)  system  using 
either  BPSX,  QPSX,  or  16-ary  QASX,  and  the  second  being  a  frequency 
hopped  C?H)  system  employing  noncoherent  FSK.  In  particular,  in 
the  latter  system,  the  MARX  and  SPACE  channels  were  allowed  to 
overlap  (i.e.,  they  were  not  necessarily  orthogonal)  in  order  to 
provide  mere  slots  in  a  given  bandwidth  over  which  to  hop. 
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Consider  the  DS  systems.  These  systems  were  analysed  in 
the  presence  of  both  tone  jamming  and  broadband  noise  jamming. 

It  was’  found  that  for  tone  jamming  the  QPSK  modulation  technique 
performed  the  best,  while  with  noise  jamming  BPSK  and  QPSX  per¬ 
formed  identically,  both  of  them  outperforming  16-ary  QASJC. 

Physically,  what  is  happening  is  the  following.  Since  a 
higher  level  alphabet  (e.g.,  16-ary  QASX)  has  a  signal  constel¬ 
lation  which  in  general  has  a  denser  location  of  signal  vectors 
(under  a  constraint  of  either  fined  average  power  or  fined  peak 
power) ,  the  higher  order  system  pays  a  definite  penalty  with 
respect  to  thermal  noise.  On  the  other  hand,  for  a  fined  infor¬ 
mation  rate  the  bandwidth  of  a  higher  order  system  can  be  reduced; 
hence  it  has  the  advantage  of  a  higher  processing  gain  in  a 
fined  overall  bandwidth.  Thus  there  are  two  apposing  forces  at 
work,  and  for  the  systems  considered  here,  the  former  one  was  the 
dominant  one.  That  is,  the  greater  vulnerability  of  the  16-ary 
system  to  thermal  noise  outweighed  the  decrease  vulnerability  of 
the  system  to  the  jammer. 

Figures  7.1  -  7.4  show  plots  of  probability  of  symbol  efror 
for  the  QPSX  system  under  different  jamming  threats.  Figures  7.1 
and  7.2  correspond  to  a  tone  jammer,  and  Figures  7.3  and  7.4 
correspond  to  a  noise  jammer.  One  qualification  has  to  be  men¬ 
tioned  with  respect  to  Figures  7.1  and  7.2.  The  tone  jammer 
for  those  figures  was  located  at  the  carrier  frequency  of  the 
QPSK  signal,  and  an  entire  period  of  the  spreading  code  was 
contained  in  one  data  symbol.  As  described  in  [7.2],  this 


combination,  while  allowing  prec.ie  results  to  be  calculated, 
leads  to  an  overly  optimist'-  answer  with  respect  to  the  situation 


of  say  a  tone  jammer  locat  roximately  1/7  He  away  from 

the  carrier  frequency  where  s  the  symbol  duration.  In  other 
words,  for  this  latter  situation  the  ratios  of  rammer  power-to- 
signal  power  (J/S)  for  which  the  results  of  Figures  7.1  and  7.2 
apply  would  be  much  smaller  than  what  is  indicated  on  the  figures. 

With  respect  to  the  FSK-FH  systems,  it  was  found  that  some 
type  of  error- correction  coding  was  needed  in  order  for  the  system 
to  perform  satisfactorily  in  the  presence  of  a  partial-band 
jammer.  Otherwise  the  jammer  could  put  all  its  poweT  in  just  a 
single  slot  and  cause  an  error  rate  of  approximately  1/N,  where 
N  is  the  number  of  slots  over  which  the  signal  can  hop,  indepen¬ 
dently  of  the  input  signal- to-noise  ratio-  Because,  of  this,  the 
system  was  evaluated  with  a  (7,4)  Hamming  cude  and  a  (23,  12) 

Golav  code.  As  is  well  known  ([7.3]),  the  former  code  corrects 
all  single  errors  and  the  latter  code  corrects  all  combinations 
of  three  or  fewer  errors.  Results  of  system  performance  for  an 
orthogonal  FSK  system  are  shown  in  Figures  7.5  and  7.6.  When  the 
Golay  encoded  system  was  evaluated  for  both  a  251  and  a  50%  overlap 
between  the  MARK  and  SPACE  channels  (and  hence  a  corresponding 


increase  in  slots  over  which  to  hop)  it  was  found  there  was  very 
little  difference  in  performance  from  what  is  shown  in  Figures 
7.5  and  7.6  ([7.2]). 

It  is  now  possible  to  put  these  results  in  perspective  with 


the  results  of  Chapter  3.  The  question  that  was  left  unanswered 
at  the  end  of  Chapter  3  was  how  much  does  one  gain  in  the  way  of 
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system  improvement  by,  say,  reducing  the  number  channels  by 
505  and  correspondingly  increasing  the  processing  gain  by  305. 

From  the  figures  it  is  clear  that,  depending  upon  the  processing 
gainE|J/NQ,  and  J/S,  one  could  pick  up  anywhere  from  0  d3  to 
tens  of  d3’s  in  improvement. 

For  example,  from  Figure  7.3,  corresponding  to  J/S  *  13  d3, 
if  the  probability  of  error  of  the  system  is  10 ’4,  going  from 
311  channels  to  1000  channels  results  in  an  improvement  of  about 
0.1  dB.  From  Figure  7.4  however,  corresponding  to  J/S  *  2  3  d3 , 
the  improvement  at  the  same  error  rate  is  essentially  unbounded, 
since  the  curve  for  M  a  311  has  started  flattening  out  at  around 
2.6  x  10 .  Therefore,  while  the  actual  improvement  in  system 
performance  that  one  obtains  by  dynamically  varying  the  number  of 
channels  in  use  depends  upon  the  specific  operating  conditions 
of  the  channel,  it  is  seen  that  for  any  system  that  is  being 
highly  stressed,  the  improvement  can  be  enormous. 

In  summary  then,  using  either  coherent  QPSK  in  a  DS  spread 
spectrum  system,  or  noncoherent  FSK  in  a  FH  system  (in  conjunction 
with  an  error-correcting  code)  can,  when  combined  with  a  dynamic 
channel  capacity  allocation  algorithm,  result  in  significant 
improvement  of  performance  of  a  system  being  severely  jammed. 
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CHAPTER  8 


Conclusions  And  Recommendations  Pot  Further  Work 

In  this  document  we  reported  on  the  analysis  of  a  variety 
of  techniques  for  reducing  ECM  vulnerability  of  tactical  communi¬ 
cation.  The  specific  techniques  chosen  for  study  included: 

1)  The  use  of  traffic  control  and  channel  capacity  allocation 
to  permit  certain  classes  of  users  to  maintain  communication  under 
severe  jamming  conditions. 

2)  A  technique  using  extreme  value  theory  for  the  purpose  of 
monitoring  the  bit  error  rate  which  can  then  be  used  in  conjunction 
with  Cl)  to  provide  a  completely  adaptive  and  dynamic  channel 
allocation  scheme  that  will  respond  directly  to  jamming  threats. 

3)  The  search  for  more  bandwidth  efficient  digital  modulation 
schemes  which  would  provide  a  net  gain  in  jamming  protection  when 
combined  with  spread  spectrum  coding. 

4)  The  possible  use  of  a  more  protected  frame  synchronization 
scheme  in  multichannel  digital  communications. 

3)  A  means  for  exploiting  the  possibilities  of  coding  gain 
by  the  use  of  low  rate  error  correcting  codes  as  an  alternative 
to  conventional  spread  spectrum. 

6)  Techniques  for  adaptive  null- steering  antennas  to  provide 
significant  AJ  protection  independently  of  any  spread  spectrum 
gain. 

7)  Various  frequency  hopping  techniques  and  novel  ways  to 
remove  some  of  the  disadvantages  of  this  method  of  achieving 
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spread  spectrum  for  jamming  protection. 

The  choice  of  topics  was  based  on  the  special  constraints 
imposed  on  systems  design  for  tactical  communications  in  general 
and  on  multichannel  radio  communications  in  particular.  For 
example,  the  communication  nodes  can  be  expected  to  be  small, 
stationary  and  sometimes  located  in  isolated  positions;  available 
bandwidth  for  spreading  is  not  expected  to  be  large  ( if  it  is 
available  at  all);  the  jammers  are  not  expected  to  be  highly 
sophisticated  even  though  they  can  be  expected  to  know  the  weak¬ 
nesses  and  vulnerabilities  of  the  communications  systems;  and  it 
is  assumed  that  it  is  more  desirable  to  reduce  the  grade  of 
service  for  communications  under  a  jamming  threat  in  some  propor¬ 
tion  to  that  threat  than  it  is  to  allow  the  system  to  degrade 
totally  below  acceptable  levels. 

The  analysis  carried  out  to  date  shows  significant  promise 
for  achieving  the  goal  of  tactical  ECCM  within  the  constraints 
imposed.  In  particular,  we  can  draw  conclusions,  with  some 
confidence,  about  the  effectiveness  of  the  techniques  studies  in 
this  final  report.  They  are: 

1)  It  is  possible  to  effect  a  rational  and  reasonable  re¬ 
allocation  of  communication  channel  capacity  assignment  which 
meets  the  threats,  imposed  by  a  jammer,  of  possibly  disrupting  all 
communications.  A  method  for  reallocation  of  channel  capacity 
based  on  offered  traffic  and  priority  classes  of  users  has  been 
proposed  and  has  been  illustrated  with  examples.  These  examples 
show  that  a)  pooled  channel  capacity  requires  fewer  total  channels 
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for  a  given  traffic  and  grade  of  service;  b)  under  jamming  con¬ 
ditions  it  is  possible  to  reassign  the  availability  of  channels 
based  on  priority  so  that  all  users  get  protection  and  high  prior¬ 
ity  users  retain  their  grade  of  service  (probability  of  being 
blocked)  while  lower  priority  users  have  this  grade  and  service 
somewhat  degraded.  Furthermore,  it  is  possible  to  monitor  bit 
error  rate  by  using  the  methods  of  extreme  value  theory  to  obtain 
rapid  indication  of  the  threat  of  being  jammed  so  that  appropriate 
reallocation  of  channel  capacity  can  be  effected.  The  combination 
of  channel  bit  error  rate  monitoring  and  the  method  of  realloca¬ 
tion  based  on  priority  users  can  result  in  a  dynamic  and  adaptive 
scheme  for  making  most  effective  use  of  the  available  bandwidth 
based  directly  on  real  time  traffic  needs,  priorities  and  jamming 
conditions. 

2)  Detailed  calculations  have  been  made  which  now  ^llow  us 
to  assess  the  -tradeoffs  of  using  bandwidth  efficient  (high  bits 
per  second  per  Hertz)  signalling  alphabets  to  achieve  greater 
protection  gain  against  the  greater  vulnerability  to  interference 
and  jamming  that  these  schemes  invariably  have.  Detailed  equations 
were  derived  for  FSK,  and  M-ary  QASK  including  the  effects  introduced 
by  Gray-coding  the  constellation.  For  straight  noise  jamming  it 
appears  that  a)  reduced  MARX- SPACE  frequencies  below  the  orthogo¬ 
nality  condition  in  FSX  is  counter-productive;  b)  QPSK  appears  to 
be  the  best  tradeoff  when  J/S  is  high  (40dB) .  The  probability  of 
bit  error  equations  can  be  used  to  investigate  for  the  particular 
vulnerabilities  present  in  multichannel  systems  where  some  bits  in 
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a  frame  may  be  more  important  than  others  (e.g.  frame  synch'. 

Also  the  bit  error  rate  rate  results  can  be  used  to  assess 
the  value  of  dynamic  allocation  of  capacity  as  described  above, 
and  it  is  demonstrated  chat  gross  benefits  can  accrue  in  perfor¬ 
mance  by  dynamic  allocation  when  high  jamming  conditions  are 
present.  Detailed  curves  based  on  the  computations  are  presented 
and  can  be  used  for  other  investigative  and  design  purposes  as 
well. 

5)  Low  rate  error  correcting  codes  are  a  viable  and  often 
a  very  effective  alternative  to  spread  spectrum  for  jamming  pro¬ 
tection.  '<tfe  have  demonstrated  several  particular  classes  of 
codes  one  of  which  is  related  to  the  linear  PN  sequences  which 
are  a)  easy  to  encode;  b)  relatively  easy  to  decode;  c)  for 
hard  decision  decoding  exhibit  good  bit  error  rates  near 
Ejj/Nq^IO;  d]  perform  better  than  DS  spread  spectrum  with  the 
same  spreading  factor;  e)  provide  the  added  bonus  of  detecting 
virtually  all  errors.  Numerical  calculations  of  performance 
have  been  carried  out  to  verify  these  conclusions.  Asymptotic 
analytic  bounds  have  been  found  which  support  the  numerical 
results. 

For  soft  decision  decoding,  these  codes  are  effective  in 
providing  coding  gain  for  low  H^/Nq.  The  coding  gain  can  be 
10  dB  or  higher  but  would  require  banks  of  matched  filters. 
Another  class  of  low  rate  convolutional  codes  have  been  found 
which  exhibit  excellent  performance  and  which  show  premise  of 
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relatively  easy  decodability .  There  are  distinct  advantages  of 
using  low  rate  codes  either  in  place  of  direct  sequence  spread 
spectrum  or  as  an  adjunct  to  any  of  the  spread  spectrum  schemes 
being  used  now. 

4)  It  is  now  possible  to  analyse  the  complete  steady- state 
and  transient  behavior  of  a  linear  mean  square  (LMS)  algorithm 
adaptive  null  steering  antenna  array.  It  should  therefore  be 
possible  to  a)  determine  the  dynamics  of  the  array  when  it  is 
in  the  process  of  being  jammed;  b)  analyse  the  vulnerability  of 
the  antenna  array  to  "blinking"  or  periodic  jammers;  c)  give 
insight  to  algorithm  design  which  is  simple,  cost  effective  and 
achieves  expected  performance;  c)  provide  a  starting  point  for 
the  understanding  of  array  behavior  under  noise  and  barrage 
jammers,  CW  tone  and  burst  tone  jammers. 

Properly  designed  adaptive  antenna  arrays  offer  great 
promise  for  achieving  a  high  degree  of  jamming  protection  without 
expending  additional  bandwidth  [which  may  not  be  available  in  the 
tactical  environment.  Furthermore,  in  principle,  antenna  arrays 
can  be  "front-ended"  to  existing  equipment  without  major  modifi¬ 
cation  in  modulation  format,  channelisation,  etc. 

5)  Frequency  hopping  systems  offer  good  promise  as  an  ECCM 
technique.  Many  ideas  associated  with  frequency  hopping  need  to 
be  developed  to  exploit  the  full  potential  of  FH.  Alternate  means 
of  correlation  at  the  receiver  can  be  used  to  change  the  character 
of  the  "sliding  correlation  function"  and  permit  sidelobeless 


outputs.  The  effects  of  partial  band  jamming  on  FH  can  be  deva¬ 
stating  if  FH  modulation  is  approached  naively  since  it  allows 
a  jammer  to  optimize  a  strategy  for  maximum  communications 
damage.  These  problems  can,  however,  be  overcome  by  diversity » 
coding  and  other  techniques. 

Recommendations  for  Further  Work 
A.  Investigative 

•  Refine  algorithms  for  adaptive  allocation  of  capacity 
under  jamming  threats.  This  will  also  require  the  development 
of  an  effective  nonparametric  technique  for  sensing  jamming  con¬ 
ditions  and  estimating  3ER  "on-line". 

•  Explore  the  properties  of  the  complete  Clow  rate)  con¬ 
volutional  codes  with  a  specific  view  of  defining  a  relatively 
simple  algorithm  for  soft  decision  decoding  of  these  codes. 

•  Investigate  the  periodic  response  of  adaptive  array 
antenna  algorithms  in  order  to  uncover  means  for  avoiding  the 
vulnerability  to  blinking  jammers.  Since  the  analysis  available 
appears  to  be  able  to  handle  multiple  jamming  signal  approaching 
from  different  angles,  this  should  be  given  particular  attention. 

•  Analyze  further  the  behavior  of  adaptive  arrays  in  noise 
and  random  jamming  signals.  This  analysis  has  never  been  carried 
out  until  the  work  reported  here.  The  payoff  in  understanding  th 
response  of  the  array  to  stochastic  signals  can  be  great  since 
real  jammers  are  best  modeled  by  random  signals. 


B.  Developmental 


3ased  on  the  direct  results  of  our  research  and  also  on  the 
contact  with  other  problems  that  we  have  had  in  the  course  of  our 
investigations,  we  recommend  a  number  of  development  projects 
which  would  hopefully  exploit  the  theoretical  ideas  discovered 
and  take  advantage  of  the  opportunities  that  a  successful  VHSI 
and  VLSI  endeavor  will  produce  by  the  mid  1980's. 

•  Implementation  of  all  the  signal  processing  functions 
of  direct  sequence  spread  spectrum  receiver  on  a  chip.  This 
would  include  correlators,  matched  filters  for  acquisition, 
logic,  clocks,  sequence  generators,  data  detectors  and  alarms. 

•  Implementation  of  a  digitally  controlled  frequency 
synthesizer  in  VLSI.  This  would  be  a  challanging  task  since  fre¬ 
quency  synthesis  functions  are  predominantly  analog  such  as 
filters,  mixers,  etc.  The  new  analog  integrated  circuit  concepts 
such  as  capacitor  switching  may  be  invoked  here.  Success  in  this 
venture  would  make  frequency  hopping  systems  much  more  palatable. 

•  Develop  an  integrated  circuit  matched  filter  bank  for  the 
Reed-Muller  codes  (PN  Sequence  Codes)  and  other  block  codes  shown 
to  be  effective  over  Gaussian  Channels.  This  would  require  the 
equivalent  of  about  1000  tapped  delay  line  systems  and  could 
conceivably  be  done  in  VLSI.  The  achievement  of  integrated, 
matched  filter  banks  could  result  in  substantial  coding  gain  as 
shown  in  this  report.  Additionally,  such  a  development  could  be 


used  for  rapid  acquisition  and  other  signal  processing  objectives 
in  EC  CM. 

•  Develop  the  front  end  and  signal  processing  functions  of 
an  adaptive  array.  At  present  adaptive  arrays  are  expensive 

since  they  require  multiply  many  front  ends  and  the  signal  processing 
that  goes  with  each*  One  far  out  possibility  is  to  produce  a  VH5I 
''sampled-data''  adaptive  array. 

•  Develop- a  convolutional  decoder  in  a  very  small  package. 

This  might  be  a  maximal  likelihood  (Viterbi)  decoder  of  moderate 
constraint  length  (3-4-)  on  a  chip  or  a  hard  decision  decoder  based 
on  a  possible  algebraic  algorithm  for  the  complete  Clow  rate)  codes 
described  in  this  report. 

•  Develop  "Smart"  multiplexors  for  digital  multichannel 
communications  which  can  switch  traffic,  assign  priorities, 
operate  at  variable  data  rates  and  variable  channel  allocations. 

This  would  require  the  multiplexor  to  have  a  microcomputer; 
direct  access  memory  and  programming  capability.  One  of  many 
applications  of  such  a  smart  multiplexor  would  be  to  implement  a 
priority  based  adaptive  channel  allocation  scheme  as  described  in 
this  report. 
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Appendix  A.  Spread  Spectrum  Background  -  Frequency  Hopping  Systems 
A.l.  Introduction 

The  ECCM  technology  has  produced,  over  the  past  20  years  a 
plethora  of  ideas  for  jam-resistant  communications.  Some  indi¬ 
cation  of  the  diversity  of  the  proposals  is  provided  by  examina¬ 
tion  of  the  list  below: 


1. 

Direct  Sequence  Pseudo-Noise 

CPN) 

2. 

Frequency  Hopping  (FH) 

5 . 

Frequency- Time  Dodging 

4. 

Time  Dodging  Burst 

3 . 

Chirp  Techniques 

6 . 

Interleaved  Error  Correcting 

Coding 

7. 

Moving  Notch  Filters 

8. 

Adaptive  Signal  Cancellers 

9. 

Receiver  Blanking  Techniques 

10. 

Adaptive  Antenna  Nulling 

11. 

Weak  Signal  Amplifiers 

The  first  six  items  on  this  list  involve  some  or  a  large 
degree  of  spectral  spreading  so  that  the  amount  of  bandwidth 
required  per  bit  per  second  of  raw  data  is  larger  than  would 
normally  be  necessary.  The  other  techniques,  while  not  utilizing 
the  bandwidth  as  a  resource  for  achieving  jamming  protection 
require  a  priori  knowledge  of  the  nature  of  the  jammer  or,  in 
the  case  of  adaptive  null  steering  antennas,  an  array  receiver. 

In  this  appendix,  a  review  of  the  analysis  and  major  design 
issues  for  frequency  hopping  will  be  presented. 
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On  the  premise  that,  in  any  event,  certain  resources  must 
be  expended  to  achieve  anti- jamming  capability,  it  is  important 
to  first  address  the  question  of  how  one  measures  EC CM  effective 
ness.  Perhaps  the  most  commonly  used  measure  of  protection 
against  jamming  is  a  quantity  known  as  the  processing  gain. 

In  the  following  pages,  after  defining  Processing  Gain  and  its 
concomitant  system  measure,  Jamming  Margin ,  we  will  show  how 
singularly  misleading  this  measure  can  be  in  general  and,  in 
particular,  when  assessing  the  performance  of  frequency  hopping 
systems . 

A.  2.  Processing  Gain 

The  notion  of  processing  gain  (P.G.)  derives  from  the  fact 
that  a  signal  to  noise  (jammer)  ratio  improvement  can  be 
obtained  by  a  receiver  which  possesses  a  replica  of  the  wideband 
modulation  carrier  by  correlation  or  matched  filter  processing. 
This  is  shown  in  Figure  A. la  where,  for  illustration  purposes, 
the  broad  band  process  in  a  pseudo -noise  (?N)  sequence.  The 
pseudo-noise  at  the  transmitter  causes  the  transmission  to  be 
spread  and  a  bandwidth  B^j.  which  is  large  compared  to  the  data 
bandwidth,  3^.  Correlating  the  received  signal  with  a  replica 
(time  synchronised)  of  the  transmitted  spreading  signal  causes 
the  resulting  signal  to  ’’coherently  collapse”  in  bandwidth  to 
that  of  the  data  spectrum. 

Subsequent  post-correlation  filtering  then  permits  all 
of  the  data  spectrum  to  be  recovered  while  accepting  only  a 
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fraction  3^c/3^  of  the  interfering  signal, 
is  one  of  many  wavs  of  viewing  the  achievement  of  the  pro¬ 
cessing  gain  and  is  easily  formalised  at  baseband  as  follows: 
Transmitted  signal  «  s (t)  *  XCt)  d(t) 

Where  XCt)  *  pseudo  noise 

d(t)  ■  data 

Receiver  correlation  output  R(t)  : 

rCt.r)  -  [XCt)dCt)  •  XCt-t)]  *  h(t) 

Where  *  h(t)  denotes  convolution  with  the  impulse  response 
of  the  receiver  filter  and  ?  is  the  local  time  offset  which 
is  adjusted  to  cero  during  the  process  of  acquisition  and  main¬ 
tained  near  cero  by  tracking  loops.  The  receiver  output  is 
then 

r(t)  *  d(t)  *  h(t) 

or,  in  terms  of  the  Fourier  transform, 

RCfO  -  D(f)  •  H(f) 

so  that  if  D(f)  «  0  for  |f|>3jj  it  is  possible  to  choose 

' 1  ;  1*1  <  3n  • 


HC  £) 


>  3t 


and  thus 

r(t)  »  d(t) '  [data  recovered] 

An  interfering  signal,  not  matched  to  x(t)  however  will  yield: 

rLCt)  -  N(t)  •  X(t)  *  hCt) 

Where,  if  N(t)  ,  the  interfering  signal  is  a  narrowband  spot 
frequency  jammer,  say,  then  the  mean  squared  power  out  of- the 
correlation  receiver  becomes 
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and  if  $x Cf) ,  the  power  spectrum  transmitted  bandwidth  we 


have 
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and  finally 
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Where  *  input  interference  power  which  we  observe  to  have 
been  reduced  by  the  processing  gain,  i.e. 

B, 


P.G.  * 
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While  the  processing  gain  concept  and  the  analysis 
leading  to  it  is  valid,  the  number  has  only  a  second  moment 
interpretation-  That  is,  the  processing  gain  can  only  tell  us 
how  much  advantage  we  gain  in  average  interference  power 
that  finally  is  allowed  to  pass  through  the  correlation 
receiver.  3y  itself;  the  processing  gain  cannot  be  used  as  a 
measure  of  performance  in  digital  communications.  However, 
with  some  care,  and  under  certain  circumstances  it  is  possible 
to  obtain  meaningful  anti- jamming  performance  measures  from 
the  processing  gain.  For  example,  for  long  binary  ?N  sequences 
which,  when  filtered  behave  like  Gaussian  noise,  the  probability 
of  error  can  be  estimated  by: 

P(e)  »  f (y* P.G) 

Where  -  -  , 

s_  ,  signal  power 

'  7  j  aiming  power 
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and  f(>)  depends  on  the  modem 


e  .  2 .  for  noncoherent  FSX 


fCx) 


T  SXp-X 


CA.2) 


tor  coherent  binary  PSK 

f(x)  »  j  erf:  /x 


CA.S) 


for  binary  DPSK 

f(x)  ■  j  exp-4  (A.  4) 

as  shown  in  Figure  A. lb. 

Thus,  to  achieve  a  bit  error  rate  (3.E.R.)  on  P(s)  of  10  3 
we  require 

x  -  y  •  PG  *  9 . 6  d3 

now  if  ?G  *  30  d3  (1000:1  bandspread)  then  we  could  tolerate  J/3 
of  9.6  -  30  *  20.4  dB.  In  other  words,  the  jammer  would 
require  more  than  30  d3  or  100  times  more  power  at  the  receiver 
before  the  performance  of  B.E.R.  *  10’^  would  be  degraded'. 

Indeed,  except  for  the  ommission  of  implementation  losses  and 
the  approximation  of  the  filtered  ?N  sequence  as  Gaussian  noise, 
this  calculation  can  be  used  to  compute  a  direct  performance 
measure,  i.e.  B.E.R. 

The  protection  ratio  afforded  by  any  spread  spectrum  system 
is  theoretically  determined  by  B^/Rq  where  3-^  is  the  radio 
frequency  bandwidth  and  Rq  is  the  data  rate  (bps)  (and  is  also 
the  nominal  bandwidth  Bp  that  would  be  required  to  transmit 
binary  data  without  spreading}.  The  protection  ratio  is  simply 
the  theoretical  advantage  on  processing  gain  that  a  desired 
signal  (that  is  matched  in  frequency,  code  and  time)  has  ever 
a  jamming  signal  (wideband  noise  over  the  band,  or  a  spot, 


176 


im  ok  u  in 


single  or  multiple  frequency)  that  is  not  properly  informed. 

That  is,  for  the  same  received  signal  power  the  desired  signal 

will  emerge,  after  processing,  with  3  j^/Rq  more  power,  than  the 

jamming  signal.  Or,  to  put  it  another  way,  for  the  jammer  to 

affect  a  response  equal  to  that  of  the  desired  signal,  he  would 

have  to  provide  3^p/R^  more  power  than  the  desired  signal.  Thus 

if  G  *  then  (J/S)  out  *  G’^(J/S)  in, where  (J/S)out  is  the 
r  K0  P 

output  jamming  to  signal  power. 

The  theoretical  processing  gain,  however,  is  not,  by  itself, 
a  measure  of  how  well  the  system  is  capable  of  performing  in  a 
jamming  environment.  For  this  purpose  we  usually  introduce 
the  ideal  of  jamming  margin  (M. ) .  This  is  the  residual  advantage 
that  the  system  has  after  we  subtract  Cin  dB) : 

1.  The  minimum  CS/N) Qut  necessary  to  achieve  the  transmission 
performance  required.  (Typically,  this  would  be  specified  in 
terms  of  Hb/NQl  the  energy  (joules)  per  information  bit  divided 

by  the  '’noise”  spectral  density  in  watts/Hc  required  after  the 
processing  gain  to  achieve  a  specified  probability  of  error). 

2.  The  implementation  loss  (L^)  due  to  imperfect  timing, 
correlation,  etc. 

3.  The  loss  in  protection  (Lj)  that  would  be  incurred 

by  a  more  intelligent  jammer  chat  takes  advantage  of  knowledge 
of  the  system  which  allows  him  a  greater  degree  of  penetration 
than  that  obtainable  by  brute  force  Of  or  noise  jamming. 

Thus  in  decibels , 
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1  GP  •  _ 
*  out 

S00  KHZ ,  R  - 


*  LI  *  LJ 
15  Kbps, 


For  example,  if  3 


then 


Typically,  Lj 


1.5 


r  .  SOxlU5 
Gp  "16x10*5 

-  2  dB 


■  9,72  dB 


*r  *  Cw)  •  10  dB  CFot  BER  »  10*3,  BPSK.no  coding) 
l  o  w  out 

say,  Lj  *  2  dB  due  to  partial  band  j arming 
Then 


Mj  -  9.72  -  2  -  10  -  2  -  -4.28  dB 
This  means  that  the  system  will  operate  within  the  prescribed 
performances  (BER  <  10  *3)  up  to  the  point  when  the  jammer 
power  is  as  much  as  4.28  lower  than  that  of  the  desired  signal 
at  the  receiver  input.  If  the  jammer  power  is  larger  than 
this,  than  although  the  processing  gain  still  appears  to  favor 
the  desired  signal,  the  performance  will  drop  below  the  require¬ 
ments  (BER<10'6). 

Several  additional  points  should  be  made: 

a.  In  a  complicated  signal  structure  involving  FH,  PN, 

TH,  Gp  may  not  be  B^p/R^  but  it  certainly  will  not  be  greater 
(against  a  noise  jammer).  Gp ,  in  fact,  implies  a  reference 
jammer  which  is  usually  taken  as  stationary  Gaussian  noise 

all  of  whose  power  is  uniformly  spread  across  the  bandwidth  in 
question.  This  defines  N  in  watts/Hi  once  Bpp  is  specified. 

b.  The  loss  due  to  an  intelligent  jammer  depends  strongly 
on  the  perceived  threat.  That  is,  it  depends  on  how  much  value 
an  enemy  is  willing  to  place  on  taking  advantage  of  his  knowledge 
of  the  system  (it  is  always  a  tacit  understanding  that  an  enemy 
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knows  the  system  design  and  the  only  thing  denied  to  him  is 
the  code  of  the  day) .  If  the  enemy  choses  not  to  take  advantage 
of  this  knowledge  then  his  use  of  brute  force  noise  or  Of  will 
cost  him  Lj  d3  (or,  we  can  pick  up  the  Lj  d3  to  our  jamming 
margin  or  reduce  our  requirements,  etc.) 

C 

c.  The  Cjf)out  can  be  reduced  (thereby  increasing  Mj 
for  specified  Gp  or  vice  versa)  by  the  use  of  coding.  Typically, 
we  can  pick  up  anywhere  from  2-5  d3  with  coding.  That  is, 
to  achieve  a  specified  probability  of  error  it  may  require 
2  to  5  d3  less  out  and  this  may  be  used  to  increase  Mj 
or  reduce  G  required  for  a  given  MT.  It  should  be  noted  that 

p  j 

5>r  a  frequency  happing  system  it  is  usually  essential  that  some 

coding  be  used  because  the  protection  ratio  provided  by  that 

method,  being  effectively  an  evasion  tactic,  may  not  be  effective 

even  against  a  spot  jammer  since  all  the  jammer  must  do  is  to 

degrade  the  system  performance  beyond  the  acceptable  range. 

Thus  if  a  frequency  hopping  system  uses  100  frequencies,  a 

spot  jammer  sitting  on  one  of  them  will  cause  an  average  error 

„2 

rate  of  10  .  Even  though  it  appears  that  the  protection  gain 

achievable  is  20  dB,  the  system  is  totally  vulnerable.  Ceding 
forces  the  jammer  to  occupy  many  frequency  slots  and  since  he 
doesn't  know  the  pattern,  wiping  out  just  a  few  of  the  code 
symbols  still  leaves  the  full  power  of  the  code  to  recover  the 
data.  Since  most  block  codes  (including  Reed-5olomon  M-ary 
symbol  codes)  are  effective  provided  not  too  many  symbols  are 
caused  to  be  in  error  in  a  given  block,  it  is  usually  also 
necessary  to  interleave  the  data  before  transmission  and 
deinterleave  it  at  the  receiver  so  that  bursts  of  errors 
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caused  by  a  high  power  pulse  jammer  is  distributed  among  many 
blocks  and  the  code  retains  the  effectiveness . 


d.  The  jamming  margin  provided  above  is  based  on  spread 
spectrum  signal  structure  alone.  Additional  and  significant 
CIO  *  20  dB  or  more)  anti- jamming  protection  may  be  obtained 
from  adaptive  null  steering  antennas  which  discriminates 
against  a  jammer  because  of  the  physical  separation  of  the 
desired  transmitter  and  the  jammer.  This  technique  does  not, 
in  principle,  require  bandspreading  but  it  is  most  naturally 
implemented  in  conjunction  with  PN  spread  spectrum  which 
provides  a  necessary  reference  signal. 

e.  Finally,  the  actual  jamming  margin,  Mj  required  is  an 
operational  specification  that  is,  it  is  determined  by  the 
perceived  jamming  threat  scenario.  Thus,  if  a  requirement  of 
say,  l2dB  jamming  margin  is  based  on  the  assumptions  of  say, 
2KW  jammers  deployed  in  a  certain  way  and  with  100  Watt 
friendly  transmitters ,  then  all  other  assumptions  remaining 
the  same,  3  dB  of  jamming  margin  can  be  gained  by  increasing 
the  power  to  200  watts  per  transmitter.  The  total  electro¬ 
magnetic  "clutter"  caused  by  this  move  would,  of  course,  in¬ 
crease  by  3  dB  but  if  this  could  be  tolerated  by  systems 
occupying  the  same  band  then  one  could  presummably  obtain  the 
same  performance  against  the  specified  jamming  threat  with 


3d3  less  processing  gain  Cot  one  half  the  required  B^p 


if  we 


choose  to  trade  on  that) . 


In  this  appendix  we  hope  to  introduce  some  of  the  properties 


and  problems  associated  with  frequency  hopping  and  especially 


as  they  relate  to  the  possible  use  of  such  systems  to  minimite 
jamming  vulnerability  or  to  maximize  jamming  margin  as  defined 
above.  The  objective  is  to  record  some  analysis  which  is  not 
generally  available  or  known  and  to  point  to  the  possible 
pitfalls  that  await  naive  application  of  formulas. 

A.  3.  Frequency  Hopping 

Frequency  hopping  is  an  evasion  technique  by  which  the 
” instantaneous"  frequency  of  a  carrier  is  changed  pseudo- 
randomly  every  T  sac.  1/T  is  the  hopping  rate  and  the  pattern 
is  illustrated  in  Figure  A.  1c  where  the  frequencies  are  chosen 
from  the  set  { f^ }  corresponding  to  pseudo- random  code  i.  3ecause 
the  spectral  occupancy  of  a  finite  burst  of  carrier  is  approxi¬ 
mately  1/T,  it  is  necessary  to  space  the  separation  between 
frequencies  by  the  order  of  1/T.  Thus  if  there  are  N  frequencies 
in  the  set  over  which  happing  takes  place,  then  B^p,  the  total 
range  is  N/T  as  shown  in  Figure  A. 2  where  the  time -frequency 
occupancy  is  shown. 

A  typical  implementation  of  the  Frequency  Hopping  (FH) 
transmitter  is  shown  in  figure  A. 3.  A  digitally  controlled 
Frequency  Synthesizer  (to  be  discussed  later)  is  driven  by 
a  PN  generator  to  generate  a  pseudo-random  hopping  pattern.  A 
common  way  to  impart  data  is  to  allow  the  least  significant 
bit  CLS3)  of  the  frequency  synthesizer  control  input  to  cause 
frequency  shift  keying  of  the  hopped  signal  in  order  to  generate 
mark  and  3pace  symbols.  The  FH  transmitter  could  either  be 
coherent  and  preserve  relative  phase  in  transitions  or,  as  is 
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more,  commonly  the  practice,  the  transmitter  phase  is  not  pre¬ 
served  and  we  have  non-coherent  FH.  The  FH  receiver  is  more 
complicated  since  provision  must  be  made  for  acquisition  and 
tracking.  A  simplified  block  diagram  is  shown  in  figure  A. 4. 

The  process  of  demodulation  after  code  acquisition  has  taken 
place  is  illustrated  by  examining  Figure  A.S  where  the  local 
synthesizer  effectively  "dehops”  the  signal  to  a  standard 
FSK  signal  for  demodulating  the  data. 

In  order  to  fully  appreciate  the  properties  of  FH  for 
purposes  of  acquisition  and  spectral  analysis,  we  shall 
derive  the  auto  correlation  function  and  the  power  spectral 
density  of  a  coherently  hopped  signal,  the  hopping  pattern  of  which 
is  assured  to  be  random. 

A. 4.  Autocorrelation  Function  of  Coherent  Frequency  Hopping 


Consider  a  wave  form  generated  by  selecting  at  random  one 
of  2N  sine  functions  sin  (nAwt)  n  »  tK,  ±(K  *  l)...e(K  +  N  -  1) 
every  T  seconds.  Assume  further  that  K>>l/Tw  so  that  there  are 
many  r.f.  cycles  in  any  one  interval,  T. 

The  total  waveform  can  be  represented  by 


s(t) 


K+N-l 

l 


n*K 


(t)  sin  nAwt 


(A. 5) 


where 


is  a  random  pulse  train  of  1,  o, 


and  -1  so  that 


and 


av.<  P^n^>  ■  o 
av.  PT^nJ  2>-  0 


'(A- 6) 

(A.:) 
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and 


Af  ■  ZirLv  is  the  frequency  spacing. 


A  typical  s (t)  is  shown.  in  Figure  A. 6.  This  signal  may 
be  generated  by  the  scheme  shown  in  Figure  A. 7. 

The  auto  correlation  function  $  (t)  of  s  (t)  taking  into 
account  the  fact  that  K>>1  is: 

x  K+N-l 

*sCt)  '  1  n-K  =0S  ni"T 

(n) 

Now  is  the  same  for  all  n  and  is  readily  seen 

?T  (?) 

to  be 

*pt(t)-  ^  (1  -  if1:  ;  |T|  <  i 
-  o  ;  |  t  |  >  T 

Hence, 

,  It  I  K**-1 

♦  S(T)  ■  (1  -  -Hjr*-)  I  cos  n  Awt  ;  |t  i  <  T 

n*K 

or 


0  s  (?) 


,  .  ,  r  K+N-l  K-l 

•to  (1  *  V*)  7  cos  nAwr  -  l  cos  nAwx 

*  L  n-0  n-0  J 


The  finite  trigonometric  sum  has  a  closed  form  expression*  so 
that 


os  (?) 


3  C1  ' 


cos  -*2'* 1  Awt  sin  ^to--Awt 


sm 


Awt 


JC- 1  K 

cos  to—  Awt  sin  t  Awt 


sm 


Awt 


*See  R.P.  Adams  Smithson,  Mathematical  Formulae,  p.  81. 


(A.  8) 


(A.  9) 


(A. 10) 
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Figurs  A.? 


Figurs  A. 3 
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UNCLASSIFIED 


-  1 

TK 


[sin  Aw.  .  sin  2^1  .VT] 


sin 


iWT 


$s  (") 


1 

itt 


Awt 
sxn  -y 


sin  t  Awt  ^  cos-^yli-  awt;|t|  <  T  (A- 11) 


<  T 


Since  (2K+N-1)  >>  N  the  bracketed  term  is  recogniced  as  the 
envelope  of  the  autocorrelation  function.  Furthermore,  when 
N  >>  1  then 


1 

ITT 


Cl- 


sm 


is  in  the  envelope  of  the  envelope. 

Equation  (A. 11)  can  be  normalized  by  dividing  by  <j>s(0)  ■ 
1/2,  and  setting  t/T  ■  0.  Then  we  get 


*SC6)  , 

S7HTT  s  u 


3  0 


5in  f  —  cos  ifclfcl  4wTS 


sm 


AwT9 


(A. 12) 


;  |8|  >  1 
I  ©  1  >  1 


Some  interesting  results  are  easily  obtained  form  the 
bracketed  function.  For  example  the  first  significant  sidelobe 
occurs  when  j  AwTS  *  Then  if  N  is  large. 


1 


5A~t 

ITT 


(a. is: 
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♦  »  0.212  (N  very  large) 

3-(?) 

whe re {—=-7  v:is  the  bracketed  envelope  ractor  or  the  ncmalited 
$sLo) 

autocorrelation  function  (A. 13). 

In  order  to  avoid  a  large  spurious  peak,  it  is  necessary 

that 

sin  j  dwTt  3*  0  0  _<  3  <  1 

It  is  sufficient  that 

7  AwT  <  7 
or, 

dW  <  ijr  or  df  <  -»rijr  :,A.Id) 

There  are  2N  lobes  and  an  equal  number  of  nulls. 


A.S .  Spectral  Density 


The  power  spectral  density  is  obtained  directly  from  the 
Fourier  transform  of  (6)  The  Fourier  transform  of  a  typical 
term  is 


F{jL  (i  .  Ijl)  cos  awj}  -  \  (Wp  )f-aAf)  *  W^Cf+ndf)) 

where  W0  is  the  Fourier  transform  of  (t) 

*T  *  T 

Now, 


-  F{J  Cl  -  -UU 


T  , sinrrfT,  2 

ts  c“rrrr”J 


!  c\ 


<  T 


<  a 


Hence  the  two-sided  spectral  density  of  the  frequency 
hopped  signal  is: 


130 


Ws(f) 


T  ,sinrTCf-nlf).2  .sinrTf  fmhf) .  2 

IoN  '•r  (f-nix3T  ""  L-rr(f+ndf)  T  ' 

u*K 


A  typical  spectrum  for  the  case  t£  •  y.  N«4  is  shown  in 
Figure  A. 8. 

Notice  that 


CA.15 


J 


/  W  ( £)d£  -  7  -  0,(0) 

-CD 

which  is  the  mean-squared  power  of  the  signal. 


Example 

For  the  case  N  *  16  and  the  frequency  separation  d£ 

is  — ,  the  envelope  normalized  autocorrelation  function 
T 

{®16  CS)  }/<P16  (o)  is  plotted  in  Figure  A.  9.  The  dotted  curve 
is  the  envelope  of  the  sidelobe  peaks.  The  first  significant 
sidelobe  has  a  normalized  amplitude  of  0.212  which  was  shown 
to  be  the  maximum  possible.  The  reason  for  the  large  central 
peak  is  due  to  the  fact  that  Af  does  not  satisfy  condition 
CA.14).  Physically,  the  peak  occurs  because  the  frequencies 
chosen  all  have  a  period  which  is  an  exact  even  submultiple 
of  the  quantizing  period  T.  Therfore  at  exactly  a  shift  of 
t  ■  T/2  or  9  ■  1/2,  all  the  frequencies  are  in  phase  again. 

This  central  peak  would  not  occur  if  if  ■  1/T. 

The  curve  of  Figure  A. 9  can  be  applied  directly  to  a 
practical  frequency  hopping  system.  For  example,  it  applies 
to  a  system  using  16  frequencies  separated  by  100  KC  and  hopping 
at  a  SO  KC  rate. 

The  large  central  peak  cf  Figure  A. 9  would  be  extremely  undesirable 
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for  acquisition  since  a  false  lock  condition  would  occur  with 
high  probability.  This  coupled  with  the  fact  that  the  spectrum 
is  not  smooth  (see  figure  A. 8)  neccesitates  that  we  choose 
Af  *  i-  which  insures  1)  no  spurious  peak  2)  smooth  spectrum 
3^  orthogonality  of  the  frequency  bursts  for  good  discrimination. 
For  Af  »  y,  the  autocorrelation  function  is  shown  in  Figure  A. 10. 

A. 6.  Alternate  Derivation  of  the  Autocorrelation  Function 


The  derivivation  given  above  is  appealing  because  of  the 
fact  that  it  follows  the  physical  means  by  which  the  signal  is 
created.  There  are,  however  certain  rough  spots  so  that  the 
alternate  more  rigorous  derivation  is  warranted. 

The  autocorrelation  of  stochastic  process  s(t)  is  defined 


as 


4>s  *  E[s(t1)  x  s(t2)] 


PCs^s 


where  E  is  the  statistical  expectation  and  p(s 
probability  density  of  s  Ct^)  and  s(t2)  and 
*$s(ti,t2}  «  E[sCtx)]  x  E[s(t2)]  -  o 


2)s1s2dslds2 
1  s2)  is  the  joint 


since 

E[s(tx)]  -  S[s(t2)]  •  0 

If  tx  and  t2  are  in  the  same  interval,  then  p($x,s2)  ■  0 
except  when  s(t^)  and  s(t2)  are  on  the  same  curve  in  which  case 
p(s^,s?)  is  the  probability  of  the  curve.  If  there  are  2N 
curves 


(ti,t2) 


2N 

l  (tn)s„Ct,) 

n»l 


Where  ?  is  the  probabilit 
n 

of  the  curve. 
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9  91 


if  all  the  curves  are  eauiorobable  P„  •  — .  Now  if 

n  2N 

s  »  sin  Awt,  n  *  =K,  =(K+1),...,  =  {K  ♦  N  -  1) 
n 

,  *C|*N-lD 

$sCti»t?)  *  I  (sin  nAwt^  sin  nAw  Awt,} 

n-iK 

l  i(r+N-l) 

*s  ^1  ’  t2'*  *  477  {cos  nAw  (t2'tl^  '  cos  nAw  ^t2*tl^ 

and 

1  *CK+N-1) 

$ j  C”  » " i )  *  l  {cos  nAwr  -  cos  nAw  (2t1  ♦  r) ; 

n*-K 

i  KiN_ 1 

yv  l  {cos  nAw  r  -  cos  nAw  (2t,  +  f) } 

‘  n-K 

If  the  process  is  stationary  t.  is  distributed  uniformly 

* 

over  [Q»  f] .  Therefore,  we  wish  to  averate  *s  (T.t^)  over  t^. 

This  is  done  by  integrating  over  the  shaded  areas  of  Figure 
A.  11  where  $s(t1,t2)  and  $s  (t.t^  possess  non- zero  value. 

Hence  for  t  >  0 


4>s(t)  -  y  ‘  $sCt,  tx)  dt:  ;  0  <  T  <  T 


,  K+N- 1 


2  (T  -t)  cos  nAw-r  +  sin  nAw  (2tx  *t)  ]lq 

n*K  o 


K+N-l 


*  (?)  “  78  I  {(I*  y)  cosAwT  ■  J nTSw  sin  nAwT} 
n-K 

if  K  >  >  ,  the  second  term  may  be  dropped  and 

.  _  K+N-l 

cs(t)*  (i  -  *0  L  cos  nAw:  ;  9  <  t  <  T 

n»K 


(A. 16) 


(A.1-) 


T  >  T 
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Since  $  (t)  is  an  even  function  (A. 17)  iS  identical  with  (A. 10) 

A  7.  Correlation  Function  for  Non-coherent  FH 

For  reasons  having  to  do  with  rapid  acquisition,  it  is 
common  to  have  the  receiver  in  a  FH  system  operate  incoherently. 
That  is,  the  correlation  operation  for  dehopping  is  followed 
by  a  square  law  or  an  envelope  detector  so  that  the  phase  of  the 
r.f.  signal  is  not  relevant  to  the  detection  process.  The 
consequence  of  this  is  that  the  fine  structure  of  the  corre¬ 
lation  function  is  obliterated.  A  simple  calculation  can  show 
this  phenomenon.  The  diagram  of  figure  A. 12  shows  the  basic 
receiver  structure  where  the  local  FH  synthesizer  is  s ' Ct  *  t) 
and  is  a  “-shifted  version  of  the  received  signal  s(t)  and 
has  random  phase  -relative  to  the  received  signal. 

Thus 

s(t)  «  i  p^(t)  sin  n  dwt  (as  in  A. 5) 
n 

s'  Ct)-  l  p^(t)  sin  (ndwt  ♦  ?  ) 
n  n 

¥  are  uniform  (0,2ir);  i.i.d. 
then  the  (cross  correlation  function  of  s  and  s'  is 

* 

r__,  (t)  »  s(t)  s'  (t+T)dt 

^period 

*  S  T  I  cos  (ndwr+'?n) 

n 

for  J  r |  <  T  and  appears  at  point  A.  Next,  consider  that  the 
non-coherent  detector  performs  a  squaring  and  averaging  so 
that  the  output  (at  B)  is 


»„,(-)  -  H[r2(T)3  -  [  |(l-ifL]  *j  I  (y  *  f  cos(2nAwT*2*n)] 

The  overbar  indicates  expectation  or  average  and  is  clearly 
equal  to  cero  so  that  finally, 

iss<t)  -  4  Cl  -  41)2;  1*1  <  T 

as  shown.  The  lack  of  fine  structure  allows  for  much  cruder 

timing  accuracy  than  would  otherwise  need  to  be  the  case  for 

acquisition  or  tracking.  The  non- coherent  autocorrelation  function 

resembles  that  of  a  PN  direct  sequence  with  one  exception.  The 

reciprocal  of  the  width  of  the  autocorrelation  function  which 

is  one  half  the  hopping  rate  is  not  proportional  BRp,  the 

bandspread.  Furthermore,  the  processing  gain  G^,  determined  earlier 

does  not  depend  on  the  hopping  rate  and  thus  for  a  given  Bj^. 

and  data  rate,  we  can  achieve  any  processing  gain  desired 

even  with  arbitrarily  low  hooping  rate’.  This  is  easy  to  see 

B 

since  the  processing  gain  G  ■  5^ —  ■  ^73  *  N  when  the  data 

P  aDATA 

is  modulated  at  1  symbol/hop  (we  shall  see  the  evil  conse¬ 
quences  of  this  later] .  Thus  G_  is  equal  only  to  the  number 

P 

of  frequencies  usedl  This  is  of  course  appealing  since 
it  implies  that  if, for  example,  we  use  N  *  1000  frequencies 
and  the  Jammer  doesn't  have  code  tracking  ability,  his  single 
spot  frequency  jammer  would  cause  interference  only  1  in  1000 
thus  giving  the  desired  signal  a  30dS  advantage.  The  folly 
of  this  reasoning  is  that  without  coding,  that  single  spot 
jammer  can  induce  an  error  rate  of  1  in  1000  which  may  be 
intolerable.  As  we  will  see  later,  this  necessitates  the  use 
of  coding  for  ?H  systems.  Without  coding  however,  the  achieve- 
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ment  of  high  processing  gain  is  an  illusion  as  far  as  jamming 
protection  is  concerned. 

A. 3  Phase  Modification  far  Improved  Correlation  in  Coherent 

Frequency  Hopping 

In  section  A. 4  we  derived -the  autocorrelation  function  for 
coherent  frequency  hopping.  The  autocorrelation  function  is 
characterised  By  many  sidelobes  which  can  present  great  difficul¬ 
ties  during  sequential  acquisition  and  also  during  tracking  when 
"tau"  modulation  or  delay  locking  is  used.  The  problem  of  side- 
lobes  is  avoided  By  non-coherent  correlation  as  was  shown  in 
A. 7,  But  with  potential  loss  in  performance,  especially  in  pro¬ 
cessing  gain  against  a  partial  hand  jammer.  It  is  shown  here  that 
a  modification  of  the  phase  structure  of  the  local  frequency 
synthesizer  relative  to  that  of  the  transmitted  Cor  received) 
signal  can  Be  used  to  eliminate  sidelobe  effects  even  when 
coherent  correlation  is  used.  Generalizations  are  possible. 

As  in  A. 4  we  write  the  input  frequency  hopped  waveform  as 

'  K+N-l 

S-Ct)  ■  l  sin  n  dwt  ■  [A.  3] 

1  a-K  1 

The  local  reference  frequency  hopping  signal  is 
5C*N-1 

S4Ct)  -  l  ?T°°  sinEniCc+T]^]  [A.  13] 

n*  SC 

When,  as  before  is  a  random  pulse  train  of  1,  0  and  -1 

and  are  a  selection  of  deterministic  phase  functions  which  are 
to  be  determined  to  remove  the  sidelobes  in  the  cross  correlation 
function 
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K+r-i  ,  . 

0-(t)  *  7  (l--^ri)cos  (nAwt+9„)  [A.18- 

n-K  1  n 

The  resulting  function,  of  course,  reduces  to  the  previous 

results  if  6  »0  identically  for  all  n.  In  the  interest  of  writing 
n 

results  relative  to  the  center  of  the  band  of  hopped  frequencies, 
it  is  convenient  to  make  the  substitution  n-M±r  where  r  is  the 
harmonic  number  around  the  center  frequency  MAw.  Then  for  N  even, 

0_(-O  ■  (1--U-4  {  cos  [MAwt  ♦  +  el 

s  1  r»-N/2  2  r 

expanding  the  cosine  and  letting  9r«-9_,r  we  obtain* 

0c(t)  *  (1— Wot  f  cos9_cos  Cr-^)  Awx ]  [2cos  MAWt)  [A.  19] 

5  1  r-1  r 


Similarly  for  N  odd 

N-l 

|r  {  T“ 

0  (r)  »  (1 — =~)[l  +  2  l  cos8_cosrAwr] CcosMAwr)  [A. 20] 
5  A  1  * 


The  last  factor  in  each  case,  cos  MAwt,  represents  the  rf 
center  frequency.  The  other  factors  are  the  envelope  of  the 
correlation  function.  When  9r»0  identically  for  all  r  we  obtain 
the  result  in  section  A. 4.  However,  if  we  relate  cos  Sr  to  the 
coefficients  of  the  binomial  expansion  we  can  adjust  the  0s(t) 
function  to  avoid  sidelobes.  Thus,  let 


cos 


•noting  that 

N72 

l  sin[(r-l/2)AWT+S_]  »  0 
r»-n/2 
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and  substituting  in  A. 19,  we  obtain  for  N  even 
and 


M/2 


M/2 


A  COS [ {T'l/l) dWT*6_]  -  2  l  cos  9„ 
r— M/2  r  r-1 

Tie  envelope  of  the  correlation  function  is 


cos(r-l/2). 


Env  0sCt) 


cos(r-l/2)dvT 


[A 


This  series-  may  be  evaluated  in  closed  form  by  examining 
:he  binomial  expansion 

N-l 

•iZiWT/  L 

y 

T 


M-l  awt  iAwr/2  -ihwr/2 
cosN  1  -  C- - ^ - ) 


M-l 


21-M  ‘  ^N-lj  eiXAwr/2  .  #1C»-1  -K) Awt/2 


K-0 


21-M  tN-lj  fti CN-1-2K) Awt/ 2 


K*0 


since  M  is  even  this  can  be  written 


-  i  C.N-1-2K)  iwr/2 


v  /  ■> .  i 

,1-N  -  1  ,N-1 


*  C5-1-I3  « 


K-0 


where  we  Lave  changed  the  index  of  summation  in  the  second  series 


from  K  to  N-l-K.  Now  since 


N-l 

CN*_1_K 3  ve  Lave 


,  22-* 


*2  CNvX)  cos  CN-1-2K)  Avt/2 
K*Q 


letting  r  ■  N/2-K 


N/2-1  v  , 

l  ^N/2-r^  cos  C.r-1/2)  A.WT 

K-Q  4  r 


[A. 22] 


equation  A. 22  contains  the  summation  in  the  envelope  of 
the  correlation  function,  and  we  can  finally  write  for  N  even 

Env  0S(/O  -  Cl- '^3  •  2N*2  •  cos1*'1  dwr/2  [A.  23] 

and  for  N  odd, 


\j  - 1  i .  v 

cos'  Awt/2  *  2 

and  we  would  choose 


N/2-1  V  1 

Cl+2  I  CN.i  3  COS  KhWT 3 
K-l  ^-i-K 


Csee  eauation  A. 203 


cos  9r  ’  'N-J’/tr1 

In  either  case  the  envelope  of  the  correlation  function  given 
by  A. 23  in  the  range  jt!<T  and  cero  elsewhere  is  strictly  a 
monotonic  function  in  this  range  with  no  sidelobes.  A  normal iced 
plot  of  this  function  for  N»16  and  T"0.9/if  is  shown  in  figure 

-  203  - 

\  - - - 


A. 13.  Tie  required  relative  phases  of  the  hopping  frequencies 
is  shown  in  figure  A. 14. 

The  implementation  of  such  a  correlation  requires  that  the 
receiver  frequency  synthesizer  adjust  the  relative  phase  of  the 

4.  e. 

selection  frequencies  such  that  the  r  frequency  from  the  center 
of  the  band  has  a  relative  phase 

*r  *  cas'1(V/2-r)  [A-24] 

for  N  even  and 

Sr  *  corVl/Vr)  CA.  25] 

for  S  add. 

Actually,  the  binomial  phase  law[A.24]  and  [A. 25]  is  similar 
to  the  technique  used  in  phased  antenna  arrays  to  eliminate 
sidelobes.  Extensions  to  this  technique  are  possible.  If 
the  amplitude  of  the  individual  frequencies  that  are  used  for 
local  correlation  are  also  amenable  to  variation,  a  possible 
advantage  here  is  the  prospects  of  obtaining  very  narrow  auto¬ 
correlation  functions  for  frequency  hopping  which  can  then  be 
used  for  accurate  timing  and/or  ranging. 

A. 9  Partial  3and  Jamming 

It  was  established  previously  that  although  the  processing 
gain  of  a  frequency  hopping  system  is  N,  the  number  of 

frequencies  used,  that  the  probability  of  error  or  bit  error 
rate  C3E£)  is  1/N  if  we  transmit  at  the  rate  of  1  bit  per  hop  and 
the  jammer  transmits  on  one  frequency.  This  is  an  extreme  case 
of  partial  band  jamming  since  the  jammer  needs  only  to  concen- 
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A  U/'C  (Degrees) 


trace  all  Lis  power  at  one  frequency.  A  more  extreme  case 
ensues  if  tLe  frequency  Lopper  uses  non-conerent  FSK  signalling 
with,  deterministic  or  adjacent  mark  and  space  frequencies.  Under 
these  circumstances  the  jammer  need  only  detect  the  instantaneous 
frequency  and  transmit  the  other  frequency.  Thus  if  a  mark 
is  sensed  then  a  space  is  transmitted  and  vice  versa.  Indeed  this 
regimen  will  degrade  the  BER  to  a  totally  useless  1/2.  This 
illustration  points  out  the  need  at  least  to  pseudo-randomice 
the  mark-space  frequencies.  This  is  especially  the  case  under 
slow  hopping  where  a  potential  repeater  jammer  has  ample  time 
to  measure  the  received  frequency  and  take  appropriate  action. 

In  order  to  fully  appreciate  the  effects  of  partial  band 
jamming  we  will  consider  several  modulation  schemes  incorporated 
in  F.H.  systems  and  the  effects  of  a  band  of  tone  jammers. 

To  begin,  suppose  we  have  uncoded  binary  coherent  phase  shift 


keving  CPSK)  and  the  signal  is 

iFb 

sCt)  «  i  1-^  C0S(wQt) 

Then  a  jamming  frequency  will  cause  a  bit  error  if  it  is: 

1.  on  frequency  wQ 

2.  of  phase  and  magnitude  to  cause  a  polarity  change 
in  the  output  of  the  demodulation. 

The  jamming  signal,  if  censored  on  wQ  is: 


[A. 26] 


j  (t)  cos^0t*e^ 


where  E.  is  the  total  jamming  energy/bit  spread  among  K  carriers, 
and  5  is  uniformly  distributed  bewteen  0  and  2r.  The  probability 


of  a  bit  error  is  therefore  equal  to 

?(e)  •  P (jammer  produces  error  I  jammer  on  frequency) 


•  ?  i 'J  a  Timer  on  frequency] 


PCs] 


cos  3  <  0) 


:< 

5T 


N 

I 


[A. 27] 


Where  it  is  assumed  that  the  jammer  phase  is  uniformly  distributed 
(0,2t)  and  where  J  ^  E,/N  is  the  jammer  density  since  it  measures 

V  w 

the  jammer  power  at  each  frequency  if  the  total  jammer  power  were 
equally  spread  amongst  each  of  the  .V  hopping  frequencies.  From 
A. 27,  it  is  apparent  that  their  exists  a  K/N  which  maximiced  PCs). 
Treating  K/N  as  a  continuous  variable  for  simplicity  allows  us  to 
find  the  maximum.  We  can  differentiate  and  set  the  result  equal 
to  zero  when  the  maximum  falls  inside  the  boundaries  1/N<1. 

Whether  the  maximum  falls  in  this  range  depends  on  Hb/J0  as 
depicted  in  figure  A.  15.  When  E^/ J"a  is  very  small  the  PCs) 
versus  K/K  cause  is  dominated  by  K/ttN'  and  it  is  seen  that  the 
maximum  within  the  boundary  occurs  at  K/N*l.  The  actual  maximum 
value  of  for  which  this  condition  holds  is  «  0.61 

C* *1 . 9dB] .  3eyond  that  range  the  peak  PCs)  occurs  inside  the 
boundaries  until  E^/N^O . SIN.  The  resulting  maximum  is  then 
given  by 
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1/rr  cos'1  VEb/J0  ;  Hb/J0  <.  0.64;  X*»N 

0.64  <  E./J  <  0.64N;  ICC^OJ 
D  0  — 

1/ir  cos'1  ^Eh/XJQ  ;  Eb/ JQ  >  0.64M;  K*«N 

Where  X*  is  the  optimum  number  of  jamming  frequencies  to  do  the 
most  damage  i.e.  maximize  PCs).  To  illustrate,  suncose  N-1000 
frequencies,  then  for  Eb/JQ  <  540  C*2S.ld3)  there  exists  a  strategy 
for  the  jammer  which  does  not  require  him  to  spread  his  power  over 
the  entire  band  of  1000  frequencies.  If,  for  example,  Eb/J  »10 
C10d3)  then  such  a  strategy  can  induce  a  PCs)  *  0.013  which  is 
unacceptable.  Recall  that  for  Eb/NQ-10  in  Gaussian  noise,  PCs)<10' 
All  of  the  above  discussion  was  for  uncoded,  coherent  binary 
phase  shift  keying  with  random  frequency  hopping  for  ECCM,  and 
discrete,  partial  band  tone  jammers.  Coherent  3PSK  was  used  to 
illustrate  and  highlight  the  mechanism  in  a  way  which  is  trans¬ 
parent  and  easy  to  comprehend.  The  same  analysis  can  be  carried 
out  for  non-coherent  nonbinary,  F5X,  DPSX,  etc.  and  for  partial 
band  noise  jamming  as  well  with  the  same  dismal  results. 

Since  most  frequency  hopping  systems  do  employ  non-ccherent 
detection,  a  simple  calculation  with  partial  band  noise  jamming 
is  useful  in  highlighting  the  phenomenon  at  work.  For  non-co - 
herent  FSK  in  full  band  Gaussian  noise  with  power  spectral  density 
N0  watts/hi,  the  minimum  probability  of  error  receiver  achieves 

PCs)  -  \  exp  C-=b/:N0) 


max  5Cs)  £44- 

0<X/N<1  bb/Jo 
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Now  if  a  two  level1  additive  noise  has  power  density 

N’0/p  in  a  fraction  0<p<l  of  the  hand 

and 

Q  in  a  fraction  1-p  of  the  band 

Then  a  random  frequency  hopper  will  encounter  this  noise  density 
with  probability  p  and  will  be  (virtually)  noiseless  or  error-free 
with  probability  1-p.  Thus  for  this  partial  band  jammer 
PCs)  *  p/2  exp  -  (p  Bb/2NQ) 

Cp  here  plays  the  same  role  as  K/N  did  in  the  partial  band  tone 
jammer  above].  It  is  apparent  that  the  p  which  maximics  PCs) 
is  p  -  2N  /H.<1  or 

m*xP(el  Z  7  ST7K"  '  I$7 

0<P<1  b  0  b  0 

With, equality  when  Eb/NQ  *  2  in  which  case  full  band  (p«l) 
jamming  is  required  to  induce  the  given  probability  of  error. 

As  an  example,  again  take  Eb/NQ»  13. 7dB  which,  recall  results  in 
P(e)  <  10  3  for  full  band  Gaussian  noise  for  noncoherent  FSK. 

A  partial  band  jammer  occupying  a  fraction  p»0.2  of  the  band 
(optimal)  will  cause  the  error  rate  to  be  0.016  which  is  totally 
unacceptable.  A  different  way  to  look  at  this  issue  which 
is  even  more  dramatic  is  to  ash  what  Eb/NQ  is,  in  fact,  required 
with  the  worst  case  partial  band  jammer  to  achieve  1Q~3  error  rate 


'It  can  be  shown  that  there  is  no  worse  distribution  than 


two  level 


The  answer  is  iS.68dB!  The  use  of  multiple  level  signalling 
actually  makes  matters  worse  since  if  we  employed  M-ary  rSx, 
the  results  for  PCs)  would  be  multiplied  roughly  by 

The  basic  behavior  in  either  case  fcoherent  PS X,  non-co - 
herent  FSK]  is  the  same.  The  error  rate  reduces  only  inversely 
with  H|3/No  rather  than  exponentially  as  is  the  usual  situation  with 
wideband  noise.  This  phenomenon  is  characteristic  not  only  of 
partial  band  jamming  but  of  Rayleigh  fading  as  well  since  the 
effects  of  the  two  are  very  much  related.  With  this  inverse 
behavior  signalling  can  he  virtually  distrupted  with  quite  little 
power  expenditure  on  the  part  of  a  jammer.  The  only  solution  is 
to  devise  a  modification  of  the  data  encoding  on  the  frequency 
hopper  which  would  force  the  jammer  to  spread  and  so  that  he  will 
achieve  no  advantage  from  partial  band  jamming..  Two  techniques 
are  available.  One  is  diversity.  I.e.  send  data  in  a  multiple 
such  as  on  different  frequencies.  In  frequency  hopping  this  is 
readily  achieved  by  having  a  bit  sent  over  the  duration  of  many 
frequency  hopping  "chips''.  This  is  known  as  bit -splitting  and 
it  has  the  effect  of  preventing  only  a  few  spot  jammers  from 
obliterating  any  bit  but  it  reduces  the  data  rate  accordingly. 

The  second  method  is  coding.  This  too  requires  effective  bit¬ 
splitting  but  if  the  codes  chosen  are  sufficiently  powerful  then 
it  may.be  possible  to  achieve  the  desired  performance  without 
sacrificing  the  data  rate.  Indeed  diversity  itself  is  a  repe¬ 
tition  code  which  is  not  optimum  in  this  sense. 


A. 10  Coding  for  Frequency  Hoppins 

Because  partial  band  jamming  can  be  so  effective  against 
a  frequency  hopper,  it  is  natural  to  inquire  whether  channel 
coding  can  help  to  mitigate  the  hazard.  It  will  be  shown  that 
even  simple  codes  such  as  repetition  codes  will  work,  but  more 
importantly,  that  with  coding  it  is  possible  to  achieve  not 
only  the  full  practical  processing  gain  expected  of  frequency 
hopping  but  that  a  coding  gain  can  be  realized  on  top  of  everything. 

To  illustrate  the  ideas  succinctly  we  will  use  a  simple 
odd- repetition  binary  code  to  show  how  partial  band  jamming 
strategies  can  be  made  ineffectual.  Consider  a  "slow  hopping" 
transmitter  that  takes  3  hops  to  present  one  information  bit. 

Then  a  noncoherent  dehopping  receiver  can  use  a  majority  vote 
decision  rule  to  decide  a  bit.  That  is,  the  receiver  announces 
a  binary  "1"  if  two  or  more  "l"s  appear  in  a  triplet  of  three 
hops.  Since  the  three  hops  during  one  bit  are  pseudo  random,  a 
single  spot  jammer  with  power  equal  to  that  of  the  received  signal 
cannot  cause  an  information  bit  error  since  this  jammer  will  at  most 
wipe  out  or  cause  an  error  in  one  out  of  three  chips  making  up  a 
bit  and  majority  logic  decision  will  always  be  correct,  providing 
that  the  frequency  triplet  is  distinct.  It,  in  fact,  now  requires 
at  least  two  spot  jammers  to  cause  a  possible  error,  (if  the  two 
spot  jamming  frequencies  happen  to  coincide  with  2  of  the  3  hop 
frequencies  and  to  cancel  them,  due  to  fortuitous  phase  relationships . ) 

For  example,  if  there  were  N«1900  random  frequencies  and 
one  jammer,  the  probability  of  a  bit  error  is  then  Pg  »  ^  p'(l-p) 

♦  p^  =  3  x  10*6  where  p  «  10*^  is  the  probability  of  a  spot 
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match. 


This  is  to  be  comcared  with  P_  *  10  J  if  we  had  no 

-  Q 

repetition  code.  This  result,  which  gives  close  to  three  orders 
of  magnitude  improvement  is  not  surprising  since  the  data  rate 
has  been  reduced  by  i.  Therefore,  in  terms  of  processing  gain: 

B  3 

P.G.  «  B— ■RF"  »  -§£  and  if  Rr  *  l  R;  P.G. '  -  3  P.G.  This  is  a 
SEATA_  A  J 

4.7  dB  increase.  If  there  is  bandwidth  to  "bum"  then  repetition 
codes  do  offer  protection  for  a  frequency  hopper  against  partial- 
band  jamming..  A  design  example  will  illustrate. 

Design  Example 

Let  R  *  Binary  data  rate 

]3rF  *  Radio  frequency  Cone- sided)  bandwidth 
N  *  Number  of  frequency  hop  channels 

T  »  Hop  rate 

M  *  Number  of  hops  per  bit 

J  »  Number  of  spot  jammers 

It  is  assumed  that  each  of  the  spot  jammers  is  of  sufficient 
power  to  cause  a  chip  error  if  the  frequencies  collide. 

Then  Pfi  =  where  p  *  J/N  and  [x]  ■  smallest 

integer  >,  x.  Then  since  orthogonality  requires  that  N  *  BRFT  and 


R 


1 

STf 


combining  terms  we  get: 


(A. 27) 
(A. 23) 
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We  usually  specify  p  ,  R,  and  B  and  we  need  to  solve  for  M 
necessary  to  deny  the  jammer  a  partial  band  advantage.  Instead 
we  will  illustrate  the  behavior  with  M  of  Pfi  with  the  B^j.  as 
a  parameter.  This  is  illustrated  in  figure  A. 16  foT  B^p  * 
lOMHr  and  40  MHz.  The  minimum  at  10  MHz  occurs  at  about  M  *  3 
while  at  40  MHz  is  minimized  at  about  M  •  11.  The  information 
rate  is  fixed  at  1KBPS  and  the  number  of  jammers  is  assumed  to 
be  100.  With  even  greater  bandwidth  a  lower  Pft  is  achievable 
but  the  data  rate  remains  the  same.  The  problem  is  that  this 
repetition  code  is  very  wasteful  of  the  bandwidth  since  similar 
reductions  can  be  obtained  by  the  use  of  a  more  sophisticated 
code. 


To  achieve  a  large  distance  between  codewords  with  more  bits/ 
chips,  (more  than  1/M  as  in  repetition  codes)  it  is  advantageous 
to  use  a  multialphabet  code.  Multialphabet  codes  are  natural 
for  frequency  hopping  since  it  is  straightforward  to  identify 
each  of  the  alphabet  symbols  by  a  particular  frequency.  Fo-  q  ary 
alphabet  with  q  ■  ps  (a  prime  raised  to  any  integer  power)  the 
class  of  Reed-Solomon  (R-S)  codes  can  be  constructed  with  (n,k)  « 
(q-1,  q-d).  n  is  the  number  of  channel  digits,  k  is  the  numbeT 
of  information  digits,  q  is  the  alphabet  size,  and  d  is  the  minimum 
distance,  such  that  the  codewords  have  the  greatest  possible 
minimum  distance  for  a  specified  Cn,  k) .  In  this  sense  the  R-S 
codes  are  optimal.  For  example  if  q  *  2J  *  8;  d*6;  (n,  k)  » 

(7,  2)  then  a  codeword  consists  of  7  octal  digits,  two  of  them 
being  information  digits.  Thus  a  seven  digit  word  contains  6 
bits  and  there  are  therefore  2^  *  8^  ■  64  codewords.  Since  the 
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distance  is  6,  a  jammer  must  cause  at  least  four  errors  to  result 
in  an  undetected  symbol  error.  A  repetition  code  would  have 
required  9  chips  for  rhe  same  performance.  This  is  a  savings  of 
2/7  or  28%.  \ 

The  R-S  codes  have  Vhe  further  property  of  1)  being  almost 
orthogonal,  2)  being  easil\  encodable ,  and  3)  having  a  sufficient 
algebraic  structure  to  allow'  for  a  variety  of  decoding  algorithms. 
The  R-S  codes  are  constructed  using  a  finite  2n  or  GF(2n) 
arithmetic.  For  example  for  the  arithmetic  of  2 3  ■  8  elements 
the  multiplication  and  addition  table  is  shown  in  figure  A. 17. 

The  encoder  is  a  mapping  of  octal  doublets  (6  bits)  into  octal 

's 

septuplets.  The  seven  elements  of  the  octal  codewoTk  C^,  C2, 

Cj,  C4,  Cs>  Cg,  C7  are  obtained  from  the  two  data  elements  a,  6 
by  the  equations 


C1  •  c  +  8 
C->  *  a  +  48 


C-  -  a  +  28 
C4  «  a  +  53 


C6  -  o  +'  78 
C-  *  a+33 


where  a,  8  are  elements  of  GFC2)  and  satisfy  the  addition  and  \ 
multiplication  table  of  figure  A. 17.  The  implementation  using 
binary  registers  is  shown  in  figure  A. 18.  Reed-Solomon  and  other 
nonbinary  codes  are  ideally  suited  to  being  implemented  in  a 
frequency  hopping  system  since  the  alphabet  symbols  can  be 
represented  by  different  frequencies.  For  example,  in  figure 
A. IS,  the  eight  digits  are  represented  by  eight  offset  frequencies 
separated  by  2 A  Hertz.  Thus  the  operation  of  the  encoder  in  this 
case  and  shown  in  figure  A. IS  proceeds  from  capturing  a  block 
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Figure  A. 18 
R-S  ENCODER 
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of  6  data  bits  (2  octal  digits)  given  as  a,  3  *  101,  100  and 
yielding  seven  octal  digits  1703264.  These  later  digits  in  tun 
induce  a  frequency  hop  shown  as  a  typical  word  in  figure  A. 19. 

The  numbers  and  the  transmission  format  are  those  used  in  TATS 
(Tactical  Transmission  System) .  A  block  diagram  of  the.  transmitter 
in  figure  A. 20  shows  not  only  the  R-S  freauency  hop  code 
but  the  major  frequency  hopping  which  can  be  invoked  when  A/J 
conditions  required  it. 

Figure  A. 21  shows  the  receiver-decoder.  A  bank  of  matched 
filters  is  used  to  detect  the  R-S  symbols  (incoherently).  The 
maximum  output  is  quantized  (for  soft  decisions)  and  is  used  to 
establish  a  received  word.  The  decoding  can  be  accomplished 
algebraically  using  the  Berlekamp  algorithm  if  hard  decisions 
are  used.  The  decoder  may  require  the  ability  to  do  arithmetic 
over  (GF(8)  and  to  interactively  solve  simultaneous  equations 
over  that  field. 

A.  11  Digitally  Controlled  Frequency  Synthesis 

The  fundamental  component  in  a  frequency  hopping  system  is 
a  synthesiser  which  has  the  following  properties: 

•  Fast  hopping  rate  (of  the  order  of  lusec  or  less) 

•  Low  spurious  response 

•  Fast  settling  time 

•  Digital  control 

All  of  these  requirements  can  be  met  with  a  direct  frequency 
synthesis  technique  which  can  be  modularised  into  a  sequence  of 
(frequency)  divide  and  (base  frequency)  add  operations.  The 
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basic  idea  is  illustrated  in  figure  A. 22. 


Here  we  see  two  stages 


of  a  possible  S  stage  system.  The  X  base  frequencies  and  the 
K—  subharmonic  of  the  center  frequency,  f  may  all  be  derived 
from  a  single  master  oscillator.  At  each  stage  there  is  a 
mixer-band  pass  filter  which  serves  to  derive  the  sum  frequency 
of  H  input  followed  by  a  divide  by  X  frequency  divider.  The 
divider  is  omitted  in  the  last  stage.  A  gate  selects  which  one 
of  the  X  base  frequencies  is  chosen  at  each  stage.  Thus  the 
control  bit  at  each  stage  consists  of  log,  K  binary  digits. 

An  example  here  illustrates  the  operation  of  the  synthesizer. 
We  wish  to  construct  a  selection  of  one  of  eight  frequencies, 
spaced  by  1  (unit  frequency)  and  having  a  band  center  of  10  (fc  " 
10).  Figure  A.  23  defines  the  problem  and  shows  that  the  tw»* 
base  frequencies  f^  »  3  and  *  ?  are  required  if  3  stages  are 
used  (23  *  8).  The  details  are  further  illustrated  in  figure  A. 24- 
where  we  fallow  the  outputs  of  stages  1,  2,  and  3  respectively  to 
finally  achieve  the  abjective.  Notice  that  the  control  in  this 
instance  is  3  bits,  one  per  stage.  Had  we  employed  say  10 
stages  with  one  bit  each  we  would  have  achieved  Z1  •  1024 
frequencies.  Alternatively  we  could  have  2  bits  per  stage  and  5 
stages  for  the  same  result.  In  that  case  four  base  frequencies 
would  be  necessary  but  faster  settling  time  would  be  achieved 
since  the  signal  would  have  to  propagate  through  half  as  many 
filters  and  frequency  dividers. 

We  summarise  by  listing  the  advantages  and  disadvantages 
of  digitally  controlled  synthesisers  currently  available. 
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Figure  A. 22 

DIGITAL  CONTROLLED  FREQ  SYNTHESIS  (ITERATED 
DIVIDE  AND  ADD) 
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DIRECT  SYNTHESIS  ADVANTAGES 

•  Fast  switching  (O.Lu  sec)  and  settling  tine 

.  -7 

•  Extremely  fine  resolution  (10  "  He  or  better) 

•  Good  spurious  rejection  (60  -  100  dB) 

•  Low  phase  jitter 

•  Wide  bandwidth  (>10^  MHz) 

•  Digital  control 

•  Modular 

DISADVANTAGES 

•  Limited  percentage  3W 

•  Large  number  parts 

•  Cost 

A.  12  Frequency  Hopping  Acquisition  and  Tracking 

One  of  the  principal  reasons  for  using  frequency  hopping 
as  a  means  of  spread  spectrum  is  due  to  the  fact  that  very  rapid 
acquisition  is  achievable  without  sacrificing  processing  gain 
during  the  acquisition  process.  Perhaps  the  simplest  acquisition 
process  to  examine  is  the  sliding  sequential  search  similar  to 
that  used  in  PN  sequence  acquisition.  In  frequency  hopping, 
the  autocorrelation  function  width  (in  seconds)  depends  only  on 
the  hopping  rate  (see  figure  A. 12)  as  it  does  in  direct  sequence 
signals.  Unlike  direct  sequence  signals  however,  the  processing 
gain  is  not  determined  by  the  hopping  rate  but  rather  by  the 
number  of  frequencies  used  in  the  pseudorandom  hep  pattern. 

This  allows  us  the  flexibility  of  having  a  slow  hop  with  broad 
autocorrelation  function  and  consequent  broad  timing  tolerance  for 
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rapid  acquisition,  while  still  permitting  full  and  arbitrary 
processing  gain  to  be  achieved  by  using  a  sufficient  number  of 
frequencies.  The  acquisition  time  is  given  by 


T  -  k  S9Z  v  I  ckiPS  „  -r 

Tacq  E  cnip  x  T  sec  x  ot 
where  k  -  #  bits /chip  searched 
R  ■  #  bits/sec  (data  rate) 

T  *  duration  of  a  chip  (sec) 

St  *  timing  uncertainly  (sec) 

The  expression  for  acquisition  time  assumes  that  either  we 
have  previously  acquired  and  have  lost  signal  so  that  reacquisition 
is  attempted  with  an  accumulated  timing  uncertainly  of  ot  seconds 
or  else  that  some  prearranged  preamble  begins  the  acquisition 
sequence  but  an  uncertainty  of  St  seconds  is  introduced  in  the 
process.  The  acquisition  time  is  simply  the  time  it  takes  to 
search  the  St  seconds.  In  terms  of  chips  this  is  simply  ^  chips 

in  order  to  integrate  over  the  data  rate  so  that  full  correlation 

v 

processing  gain  is  achieved,  we  must  spend  ^  seconds  for  each 
chip  that  is  searched.  As  an  example  for  fast  hopping,  suppose 
k  ■  3,  R  «  1K3PS  and  T  ■  ly  sec  (hop  rate  of  1  MHz) .  Suppose 
that  clock  stability  is  0.1  ppm  (10’7)  and  communication  is 
disrupted  for  10,000  seconds  (<3  hrs)  then  5  t 


10  ‘  x  104  »  1  ms. 


Then  the  clock  drift  would  have  accumulated  at  most 


1ms 


1000 


chips.  Since  we  must  spend  3  ms/chip  the  acquisition  time  is  at 
most  3  seconds.  For  the  same  example,  if  we  needed  a  processing 

4 

gain  of  -10  dS  we  woudl  use  10  frequencies  of  a  total  bandwidth 
of  10  MHz.  If  we  needed  SO  dB  we  would  use  10s  frequencies,  etc. 
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•  Constant  envelope  r.f.  signal  makes  linear  power  ampliers 
unnecessary 

•  Near  *  far  ratio  tolerable  »  10  ^ 

(*  20  dB  loss  for  line  of  sight  radio  link) 

•  Most  jammers  equivalent  to  noise 

•  Timing  accuracy  needed  <0.1 

•  Frequency  accuracy  needed  <0.1  R  (R  is  data  bank  rate) 

•  Slow  acquisition  unless  parallel  processing  on  matched 
filters  are  used 

•  Tac^  *  3  x  (P.G.)  x  <St  (5t  is  timing  uncertainty) 

•  Effective  against  "look-thru”  repeater  provided  codes  are  long 

•  High  spurious  responses  when  partial  period  correlation 
takes  place  or  when  data  rate  is  integrally  related  to  code 
repetition  rate 


-  Code  imbalance  and/or  equipment  misalignment  can  obliterate 
anti- jam  protection 

Frequency  Hotming 

•  Almost  arbitrary  band  spreading  is  possible.  Can  cover 
entire  technology  band.  Ex.  ■  1.5  GHZ  at  10  GHZ  with 

1500  frequencies  spaced  7XHZ  apart. 


•  Processing  gain  is  independent  of  hopping  rate 

•  Effectively  an  interference  avoidance  scheme  rather  than  an 

interference  averaging  system 

•  No  code  imbalance  or  equipment  alignment  problems 

•  Constant  envelope  signal 
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The  acquisition  time  would  still  be  3  seconds.  By  contrast,  if 
a  direct  sequence  system  were  used  that  required  40  dB  of  processing 
gain  with  the  same  information  rate  (lk3PS)  then  a  PN  chip  rate 
of  10 7  Hz  would  be  needed.  Using  a  serial  search  over  1  ms 
requires  that  we  search  through  1  ms/O.lus  »  10,000  PN  chips. 

At  3  ms/chip  this  would  require  30  seconds.  For  slow  frequency 
hopping,  the  result  is  even  more  dramatic.  For  example  if  the 
hop  rate  was  once  per  millisecond  then  the  acquisition  time  for 
the  frequency  hopping  system' would  be  3  ms!  Slow  hopping  has 
the  disadvantage  of  being  vulnerable  to  repeater  jamming. 

Frequency  hopping  acquisition  also  lends  itself  to  even  further 
improvement  in  acquisition  time  (as  does  PN  direct  sequence 
spread  spectrum)  by  the  use  of  multiple  processors  and  matched 
filters'. 

A. 13  Comparison  of  PN  Diiect  Sequence  and  Frequency  Hopping 

We  conclude  this  appendix  on  frequency  hopping  by  a  brief 
point-by-point  comparison  of  the  benefits  of  frequency  hopping 
and  PN  direct  sequence  spread  spectrum  modulation. 

Pseudo-Noise 

•  3inary  shift  register  generators 

•  Current-day  limit  B^  ■  100  -  400  MHz 

•  Effective  against  multipath  greater  than  1  chip 

•  Suitable  for  accurate  ranging  and  timing 

•  Selective  addressing  possible  for  multiple  access 
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*  Partial  band  jammer  may  be  effective  unless  coding  is 
employed 

•  Lends  itself  well  to  nonbinary  coding  which  is  effective 
against  burst  errors 


•  Selective  addressing  possible 


*  Timing  accuracy  required  is  <0.1  T  (Tc  is  the  chip 
inferred,  seconds)  '  <" 


•  Frequency  accuracy  required  is  <0.1  R 


-  Near  *  far  ratio  tolerable  »  10 


•  Slow  hop  rates  easy  to  achieve  and  permit  very  rapid 
acquisition 

•  Tac^  =  10  R*1  C»illiseconds) 

•  Not  suitable  for  accurate  ranging  and  timing 

•  Fast  hopping  requires  an  expensive  digitally  controlled 
frequency  synthesiser 


•  May  be  used  as  a  preamble  for  fast  acquisition  of  PN  system 
or  in  hybrid  configuration 


